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Abstract 

In this PhD thesis, we analyze and generalize the renormalization group approach 
to the resummation of large logarithms in the perturbative expansion due to soft 
and collinear multiparton emissions. In particular, we present a generalization of this 
approach to prompt photon production. It is interesting to sec that also with the 
more intricate two-scale kinematics that characterizes prompt photon production in 
the soft limit, it remains true that resummation simply follows from general kine- 
matic properties of the phase space. Also, this approach does not require a separate 
treatment of individual colour structures when more than one colour structure con- 
tributes to fixed order results. However, the resummation formulae obtained here 
turn out to be less predictive than previous results: this depends on the fact that 
here neither specific factorization properties of the cross section in the soft limit is 
assumed, nor that soft emission satisfies eikonal-like relations. We also derive resuma- 
tion formulae to all logarithmic accuracy and valid for all values of rapidity for the 
prompt photon production and the Drell-Yan rapidity distributions. We show that 
for the fixed-target experiment E866/NuSea, the NLL resummation corrections are 
comparable to NLO fixed-order corrections and are crucial to obtain agreement with 
the data. Finally we outline also possible future applications of the renormalization 
group approach. 
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Introduction 



The current theory that describes the strong subnuclear interactions is the non-abehan 
gauge theory of the SU(3) local symmetry group. This theory, called Quantum 
Chromo-Dynamics (QCD), has been tested to great accuracy in many experiments. 
Its perturbative regime will be crucial at the imminent high-energy proton-proton 
collider LHC The main target of LHC, is to find signals of the Standard Model Higgs 
boson and of new physics, i.e. supersymmetry, new interactions predicted by great 
unified theories, extra-dimensions, new gauge bosons and mini-black holes. In order 
to accomplish all this, an excellent understanding of QCD is necessary, both because 
LHC is a proton collider and because QCD backgrounds must be described accurately. 

In perturbative QCD (pQCD), it is well known that, when one approaches to 
the boundary of the phase space, the cross section receives logarithmical-enhanced 
contributions of soft and collinear origin at all orders. These large tems at order 
0{a'^) of the strong coupling constant in pQCD have, in general, the form 



log'"(l 



, m < 2n — 1, 



where 2; is a parameter that becomes close to 1 near to the phase space boundary 
and "+" denotes the plus distribution. These terms become important in the limit 
2; — s> 1 spoiling the validity of the perturbative fixed-order QCD expansion. Hence, 
they should be resummed to all-orders of pQCD in order to get accurate predictions 
of the cross sections of QCD processes. An example of the importance that can have 
resummation is given by the Higgs boson production at LHC near the threshold of its 
production. In this case z is given by the ratio of the searched Higgs mass over the 
center-of-mass energy of the partons in the colliding hadrons and the large logarithms 
arise from soft-gluon emissions. 

These large terms have been resummed a long time ago for the classes of inclusive 
hadronic processes of the type of deep-inelastic scattering (DIS) and Drell-Yan (DY) 
[H, [21 E] . Threshold resummation of inclusive processes can affect significantly cross 
sections and the extraction of parton densities [U [5] . For the case of small transverse 
momentum distributions in Drell-Yan processes, it has been shown that resummation 
is necessary to reproduce the correct behavior of the cross section [6]. 

These results has been obtained using the eikonal approaximation of Ref. [2] which 
generalizes to QCD the Sudakov exponentiation of soft photons emissions in electro- 
dynamics [7] or assuming suitable factorization properties of the QCD cross section 
[1]. More recently two other approaches to resummation have been proposed. The 
first is the renormalization group approach of Ref. [8] and the other is the effective 



V 



field theoretic (EFT) approach of Refs.jHl ^U\. There the resummation of the large 
logarithms for full inclusive deep-inelastic and Drell-Yan processes is obtained. 

In this thesis we will concentrate mostly on the study of the renormalization group 
approach of Ref. and on its applications and generalizations. This approach has the 
advantage of being valid to all logarithmic orders, and self-contained, in that it does 
not require any factorization beyond the standard factorization of coUinear singular- 
ities. It relies on an essentially kinematical analysis of the phase space for the given 
process in the soft limit, which is used to establish the result that the dependence on 
the resummation variable only appears through a given fixed dimensionful combina- 
tion. This provides a second dimensionful variable, along with the hard scale of the 
process, which can be resummed using standard renormalization group techniques. 
Beyond the leading log level, the resummed result found within this approach turns 
out to be somewhat less predictive than the result obtained with the other methods. 
In the other approaches references resummed results at a certain logarithmic accu- 
racy is fully determined by a ceratin fixed order computation, whereas a higher fixed 
order computation is needed to determine all coefficients in the resummed formula 
of Ref. [H] . The more predictive result is recovered within this approach if the depen- 
dence of the perturbative coefficients on the two dimensionful variables factorizes, i.e. 
if the two-scale factorization mentioned above holds. 

We shall show the generalization of the renormalization group approach to the 
resummation of the inclusive transverse momentum spectrum of prompt photons 
produced in hadronic collisions in the region where the transverse momentum is close 
to its maximal value. Prompt photon production is a less inclusive process than 
Drell-Yan or deep-inelastic scattering, and it is especially interesting from the point 
of view of the renormalization group approach, because the large logs which must be 
resummed turn out to depend on two independent dimensionful variables, on top of 
the hard scale of the process: hence, prompt photon production is characterized by 
three scales. The possibility that the general factorization Ref. [TT] might extend to 
prompt photon production was discussed in Ref. [12j, based on previous generaliza- 
tions |JL3j of factorization, and used to derive the corresponding resummed results. 
Resummation formulae for prompt photon production in the approach of Ref. [2] were 
also proposed in Ref. [H] , and some arguments which might support such resummation 
were presented in Ref. [15] . Our treatment will provide a full proof of resummation to 
all logarithmic orders. Our resummation formula does not require the factorization 
proposed in Refs. [121 [E] , and it is accordingly less predictive. Because of the presence 
of two scales, it is also weaker than the result of Ref. p] for DIS and Drell-Yan pro- 
duction. Increasingly more predictive results are recovered if increasingly restrictive 
forms of factorization hold. 

Moreover, we shall prove an all- logarithmic orders resummation formula for differ- 
ential rapidity of Drell-Yan and prompt photon production processes. The differential 
rapidity Drell-Yan cross section is used for the extraction of the ratio d/u of parton 
densities. The accurate knowledge of these functions is needed to study Higgs boson 
production and the asymmetry W^. The resummation of Drell-Yan rapidity distribu- 
tions was first considered in 1992 [16]. At that time, it was suggested a resummation 
formula for the case of zero rapidity. Very recently, thanks to the analysis of the full 



NLO calculation of the Drell-Yan rapidity distribution, it has been shown [T7j, that 
the result given in [16] is valid at next-to-leading logarithmic accuracy (NLL) for all 
rapidities. In this thesis, we shall give a simple proof of an all-order resummation 
formula valid for all values of rapidity. To do this, we will use the technique of the 
double Fourier- Mellin moments developed in [18]. In particular, we will show that 
the resummation can be reduced to that of the rapidity-integrated process, which is 
given in terms of a dimensionless universal function for both DY and and 
production. Then, we implement numerically the resummation formula and give pre- 
dictions of the full rapidity-dependent NLL Drell-Yan cross section for the case of the 
fixed-target E866/NuSea experiment. We find that resummation at the NLL level is 
necessary and that its agreement with the experimental data is better than the NNLO 
calculation of Ref . [I9] . In this case, we find also that the NLL resummation reduced 
the cross section instead of enhancing it for the parameter choices of this experiment. 
Threshold corrections to Higgs, and production rapidity distributions at high 
energy hadron colliders have also been studied in [201 EI] ■ 

Finally, we shall discuss the application of the renormalization group approach 
to the resummation another class of large logarithms that arise in Drell-Yan pro- 
cesses for small tranverse momentum distribution of the produces lepton pairs. These 
logarithmic-enhanced terms at order have, in general, the form 



log'-lgi) 



m < 2n — 1, 



where q± is the transverse momentum of the produced Drell-Yan pair. Also in this 
case the resummed results using the renormalization group approach are less pre- 
dictive than results obtained with the approach of Ref. as it is shown in Ref. [ffj. 
Furthermore the conditions that reduce our results to those of Ref.[n] in terms of 
factorization properties is still an interesting open question. 

This thesis is organized as follows. In Chapter [1], we review the basics concepts 
of pQCD. In particular the construction of the QCD Lagrangian, the asymptotic 
freedom of strong interactions, the structure of the cross section when initial state 
hadrons are present and the evolution equation of the parton densities. In Chapter 
[21 we discuss the importance of resummation at hadron colliders. Then, we describe 
how in the various approaches the large logarithms are exponentiated and resummed. 
They are the renormalization group approach, the eikonal approximation approach, 
the approach of non standard factorization properties and the effective field theoretic 
approach. In Chapter [31 we show in the detail the renormalization group approach 
to the resummation of all inclusive deep-inelastic and Drell-Yan processes and its 
generalization to the prompt photon process at large transverse photon momentum is 
shown in Chapter |H In Chapter [5l we prove the all-logarithmic orders resummation 
formula for the Drell-Yan and prompt photon production processes. We show also 
the impact of the NLL resummation for the DY process at the E866/NuSea experi- 
ment and discuss the numerical results. Then in Chapter [HI we turn to discuss the 
generalization of the renormalization group approach to resummation in the case of 
the small transverse momentum differential cross section for the Drell-Yan process. 
Finally, in the last Chapter, we summarize and determine the predictive power of the 



resummation formulae obtained with the different approaches, namely, the fixed-order 
computation needed to determine completely a resummation formula for an arbitrary 
logarithmic accuracy. 



Chapter 1 

General aspects of perturbative 
QCD 

1.1 Quarks, Gluons and QCD 

The quarks (the constituents of hadrons) are Dirac fermions. In the Standard Model 
(SM), as far as the electroweak interactions are concerned, the properties of quarks 
and leptons are similar. Indeed, as for the leptons, the six quark flavors are grouped 
into three SUl{2) left-handed doublets 




and six SU (2) ^ singlets, which are the right-handed parts of each flavor. Both, quarks 
and leptons, interact in a similar way with the electroweak gauge bosons (7, and 
Z°) of the group SUl{2) x Uy{1)- The main difference is that each quark flavor 
eigenstate is a unitary mixing of the quark mass eigenstate, while, according to the 
SM, this is not the case for charged leptons and massless neutrinos. However, in this 
decade it has been proven that neutrinos have mass and that there is also mixing in 
the neutrino sector thanks to observations of their oscillations. 

The peculiarity of quarks is that they have a speciflc property, the color charge, 
which is absent for leptons. Indeed, a quark of a given flavor has three different color 
states with equal masses and electroweak charges. The interaction of the quarks 
is mediated by the eight gauge bosons (gluons g) of the color group SUc{5). So, 
the quarks belong to the fundamental representation of SUc{S) and the gluons to 
the adjoint one. The gauge theory of this non-abelian group is called Quantum 
Chromodynamics (QCD) and is the current theory of strong interactions. 

Specifically, a gauged SUc{^) transformation of a quark {qa{x) with a — 1,2, 3) is 
given by 

qa(x) q'^{x) = Uab{x)qb{x) (2) 

Qa{x) q'^{x) ^qbUlix), (3) 

where the 3x3 matrix Uik{x) is the fundamental representation of the SUc{'i) group 



1 



that acts on an internal space defined at each space-time coordinate x. It satisfies 



UU^ = U^U = 1, det{U) = 1. (4) 

In this section the sum over all the repeated indices is implicit and the sum over 
spinor indices is omitted for brevity. The usual exponential representation of the 
gauge transformation matrix in terms of the basis of matrices (the generators of 
SUci'i)) of the corresponding algebra smc(3) is: 

U{x) = e-^"(")-^ = e-"^")-*', (5) 

where = {ai{x), . . . , as{x)) are the eight arbitrary parameters of the gauge trans- 
formation, A are the eight elements of the basis of the algebra suc{2) (or equivalently 
the eight generators of the group) and t are the eight color operators (in analogy to 
the spin operators of the group SU{2)). It is clear that, in order to respect Eq.(jl]), 
the generators of the group must be hermitian and traceless. The normalization of 
the color operators depends on the representation r 

tiitftf) = T,5^^. (6) 

The form chosen by Gell-Mann for the suc(3) basis in the fundamental representation: 




(7) 



With these definitions, the color matrices satisfy the following relations 

[t^,t^] = tf^B^t^ (8) 
tr(t^t^) = Tf6^^, T^ = 1 (9) 

where Tp is the normalization of the color matrices in the fundamental representation 
and f^^'^ are called the structure constants of the algebra smc(3) which are totally 
antisymmetric in {A,B,C}. The independent non- vanishing structure constants are 
given by: 

r' = h r = \. f''' = -\^ f'' = l f'' = l (10) 
f'' = l, f'' = -l^ f'' = ^, f'' = ^. (11) 



Furthermore, the structure constants provide the adjoint representation of the suc(S) 
algebra (the one which has the same dimension of the algebra). Indeed, if we define 
the adjoint representation as = —if^^'^', we can verify explicitly that this is the 
adjoint representatio because 

[T^,T^] = if^BCjC (12) 
tr(T^T^) = TaS^^, Ta = 3 (13) 

The Casimir operator Cr (the one which commutes with all elements of the algebra), 
for a certain representation r, is constructed as 

tX = C4d.xd.. (14) 

Now, since the contraction of this last equation is equal to the contraction of Eq.(l6]), 
we get 

drCr = 8Tr. (15) 

In particular, we find the Casimir operators are given by 

Cf = 1, Ca = S (16) 

for the fundamental and the adjoint representation respectively. 

We shall now show that the symmetry with respect to the gauge transformations 
Eqs. ( 121131) (together with the Lorentz invariance), can be used as guiding principle to 
construct the QCD Lagrangian. We start from the usual Dirac free Lagrangian for 
each quark mass and color eigenstate 

Cd{x) = tpfaix) {i^ - mf)ijfa{x), (17) 

where / is the fiavor index and a is the color quark index. This term is not gauge 
invariant. In fact under the gauge transformations Eqs.( l2|5l) . the Dirac free lagrangian 
Eq.( fT7Il transform as 

Coix) ^ Cd{x) + ijfb [zUl{x)d^Uac{x)^ l^i^fcix). (18) 

To restore gauge invariance, we introduce a gauge field matrix Afj_ab{x) made up of 
eight gauge fields A^{x) in this way: 

A,a,{x) = t^A^. (19) 

We assign to this matrix field the following interaction Lagrangian 

Ci{x) = gsipfaA^ab^'^i^fb- (20) 

Here Qs is the gauge dimensionless coupling analogous to the electric charge in QED. 
In order to cancel the symmetry breaking term of Eq. llTSl) . the gauge transformation 
of the field matrix A^ab{x) has to be 

A^abix) Uac{x)A^,d{x)U\,{x) - -d^Uac{x)Ul{x) . (21) 

9s 



Hence, with the introduction of the field matrix A^^ab, the sum Cd + is now gauge 
invariant. To complete the Lagrangian one has to add the pure gauge invariant term, 
which is 



CG{x) = --ii{G^,{x)G^\x)), (22) 



1 
2 

where the gluon field-strength tensor is given by 

Note that the third term of the gluon field-strength gives rise to the self interactions 
of gluons and that its origin stands in the fact that the gauge group SUc{'i) is non- 
abelian. The final form of the QCD classical Lagrangian is obtained by adding the 
three pieces introduced above and using Eq.(l9]): 

^3cd(^) = -\g%{x)G''^''{x) + i>fa{x){iD,a,r - mf5ab)i'fb{x). (24) 

where 

G^, = d,At - d^Al + ^ J^^^AjA^ (25) 

and 

Df,ab = Sabdf, - igsAf^ab (26) 

is the covariant derivative in the sense that D^abi'fbi.x) transforms as ipfa{x) under 
the gauge group. 

The quantization of the classical gauge invariant theory of QCD Eq. fl2^ . can be 
done with the Fadeev- Popov procedure. This procedure takes care of the fact that the 
equation of motion of the gluon field A^ can not be inverted and this prohibits to find 
the propagator. However, this is a consequence of gauge invariance which implies that 
the physical massless gluon has only two polarizations/spin states whereas the field 
A^ has four components. To make things work, an additional constraint on the gluon 
field is introduced, the so called gauge fixing condition which uses gauge invariance to 
define properly the gluon propagator. In QCD, however, this constraint is not linear 
and one should add specially designed fictitious particles (the so called Fadeev-Popov 
ghosts) which are Lorentz scalar anti-commuting fields and appear only in the loops. 
After this procedure (usually performed in the functional formalism), we have that 
the QCD Lagrangian from which we can calculate directly the Feynman rules is given 
by (see e.g. [22j page 514): 

^QCDix) = C^icDix) - ^{d''d''A^{x)A^{x)) - -c^d^Dfc^, (27) 

where A is the gauge fixing parameter, is the complex colored scalar ghost field 
and D^^ is the covariant derivative in the adjoint representation: 

Df = 5^^d, - ^ J^^^AJ. (28) 



The quark, gluon and ghost propagators and vertices for QCD are collected in Figure 
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Figure 1.1: Feynman rules for QCD in a covariant gauge for gluons (curly lines), quarks (solid 

lines) and ghosts (dotted lines). Here A, B, C, D arc the color indexes in the adjoint representation, 
a, b, c in the fundamental one, a, /3, 7, 6 are the gluon polarization indexes and i,j are the spinorial 
indexes. 



1.2 Asymptotic freedom and perturbative QCD 



Feynman diagrams in QCD are obtained employing the vertices and propagators as 
building blocks. However, the use of diagrams makes sense only if the perturbative 
expansion in Qs is meaningful. To respect this condition, the coupling 

= |j. (29) 

the QCD analog of the electromagnetic coupling a = e'^/An, has to be sufficiently 
small. We shall now show that the method of perturbation theory in QCD are useful 
at high energy. Indeed, the coupling constant is large at low energy and becomes 
smaller at high energy (asymptotic freedom). 

The simplest way to introduce the running coupling, is to consider a dimensionless 
physical observable R which depends on a single energy scale \/Q^- This is the case, 
for example of the ratio of the annihilation cross section of electron-positron into 
hadron with the annihilation into muons where = S the center-of-mass energy. 
We assume that this scale is much bigger than the quark masses that can be therefore 
neglected. Now, dimensional analysis should implies that a dimensionless observable 
is independent of Q^. However, higher order corrections produce divergences and 
so the perturbation series requires renormalization to remove ultraviolet divergences 
that in (i = 4 — 2e dimensions are regularized as 1/e poles. This poles can be re- 
moved defining a renormalized coupling constant at a certain renormalization scale 
nl- Consequently, in the finite e — limit, R depends in general on the ratio 
and the renormalized coupling as depends on /i^; we call the renormalized coupling 
at the scale /x^, a ^2. Since the renormalization scale is an arbitrary parameter intro- 
duced only to define the theory at the quantum level, we conclude that R has to be 
//^-independent. Formally this independence is expressed as follows: 



2 N ^ 



R = 0, (30) 



where 



P{a,2) = fil—t^ (31) 



Eq.( l30i) is a first order differential equation with the initial condition (at = fi"^) 
R{l,a^2). This means that if we find a solution of Eq. llHUj) . it is the only possible 
solution. This solution is easily found defining a function ^^(Q^) such that 

asifJ^l) = «M? (32) 



and that 



H-«!ur'^. (33) 



In fact differentiating this equation, we find that 

da^m /3(a.(Q^)) 



(35) 



and that R{l,as(Q^)) is the desired solution of Eq. (|30l) . because 

d da^m d , d , d ^2 dasm d 

~ A^rTTT = ~Q = -Q 0^9 o — 77=m- (36) 



This shows that all of the scale dependence in R enters through the running of the 
coupling constant ^^(Q^). To find explicitly this function, we need to know the j3- 
function so that we can solve Eq. flM|) . The /3-function can be calculate perturbatively 
from the counterterms of the renormalization precedure and a knowledge to order a"^^ 
requires a n-loop computation. The perturbative expansion of the /^-function is given 
by: 



00 



/?(a,) = -a,5^/5„(^)""'\ (37) 



An 

n=0 



At the moment, the QCD /3-function is known to order [23j: where in the MS 
scheme, 

Po = 11 -^iV/, Pi = l02-^Nf, (38) 

2857 5033 ,, 325 ^,2 

P2 = Nf + N^f (39) 

^ 2 18 54 ^ ^ ' 

/ 149753 / 1078361 6508 

+ 3564^3 



Ps = 1^— + 3564^3)-(^^ + ^e3)iV/ 



/ 50065 6472^ \ , . 1093 

+ + 6 N^f + N^f, (40) 

V 162 81 ^7 ^ 729 ^ ^ ' 

with Nf the numbers of flavors and ^ is the Riemann zeta-function (,^3 = 1.202056903 . . . 
The two loop solution of Eq. flM|) is given by: 



as[Q ) 



l + (/3o/47r)a,(/i2)log^ 



2 

Mr 

2 



a,(/i2)log(l + (/3o/47r)a,(/i2)logg 
l + (/3o/47r)a,(/x2)log§ 



+ 0(a^Mog^%). (41) 

For simplicity, in many cases, we will use another parametrization of the coeffi- 
cients f3n, which is obtained with the substitution: 

Pn = &n(47r)"+^ (42) 

With this parametrization, the perturbative expansion of the /5-function Eq. flH7|) be- 
comes 

00 

P{as) = -J2^nas+\ (43) 

n=0 

From Eq.f l41l) . we see that as{Q'^) is a monotonically decreasing function of Q^, be- 
cause the coefficients Po and Pi are positive (with Nf < 6). The running of asiQ"^) 
has been measured with great accuracy (see Figure fL2l) . The fact that at high energy. 




Fi gure 1.2: The running of the strong couphng constant. The asymptotic freedom is confirmed 
by the experiments 



the running coupling becomes small is a peculiarity of non-abelian gauge theories and 
is called asymptotic freedom. 

Hence, perturbative QCD can be applied when the relevant scale of a certain pro- 
cess is high enough such that the running coupling becomes small. A typical example 
is given by the annihilation of two high-energy electrons into hadrons. Perturbative 
QCD can also be applied to processes in which hadrons are present also in the initial 
state thanks to the factorization theorem. According to this theorem ^11], the cross 
section for the production of some final state with high invariant mass (the scale 
at which the running coupling constant is small) with two incoming hadrons is given 
by: 




For processes with a single incoming hadron the factorization theorem is simpler. For 
example for the deep inelastic scattering (DIS) of a lepton that exchanges a high 
square momentum Q"^ with the hadron, the cross section takes the form: 



a{P,Q') = y^[ dxF^{x,fi^)aa{xP,as{fi'),Q\fi'). (45) 

a 



In Eqs. fl44|45l) . Pi is the momentum of the incoming hadron Hi. A beam of hadrons 
of type Hi is equivalent to a beam of constituents (or partons) which are quarks 



or gluons. These constituents carry a longitudinal momentum XiPi characterized by 
the parton densities F^^{xi, fi"^). That is, given a hadron Hi with momentum Pi, 
the probability density to find in Hi the parton a with momentum XiPi is given by 
F^*{xi, fi"^). Furthermore, these functions are universal in the sense that they are 
process independent. The parton densities depends also on the so called factorization 
scale /x^. This scale is introduced to separate off the non-perturbative part of the cross 
section (the parton densities) from the perturbative one This is exactly the cross 
section where the incoming particles are the partons a (and b) and can be calculated as 
a perturbative expansion in agin^)- The parton densities have a mild dependence on 
the scale /i^ determined by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equations (see section [L4l) . Here, we have chosen the renormalization scale /if equal to 
the factorization scale //^ for simplicity. Anyway, in order to reintroduce the scale /i,^, 
we have only to rewrite ^^(/i^) in terms of fif. (see Eq.pTl)) and expand it consistently 
with the order of the calculation. The /i^ dependence in the parton densities is 
compensated by the fi"^ dependence in the partonic cross section a. However, with a 
fixed-order computaction of the partonic coefficient function at order the hadronic 
cross section will still depend on /i^ with a dependence which should be of order aj"*"^. 
Hence, this dependence can be used to estimate the theoretical error of a fixed-order 
computation. A simple discussion about the dependence of the hadronic cross section 
on the factorization and renormalization scale is given in Ref . . 



1.3 NLO DY and DIS cross sections 

We consider for simplicity the classical Drell-Yan (DY) hadronic process for the pro- 
duction of a dimuon pair through a virtual photon 7* (see Figure [USD: 

ifi(Pi) + H2{P2) = 1*{Q) + X{K), (46) 

where Hi and H2 are the colliding hadrons with momementum Pi and P2 respectively, 
Q is the momentum of the virtual photon and X is any number of additional hadrons 
with total momentum K. For the process of Eq.( H6i) . we define 



X ^ ^, (47) 

where S = {Pi + P2)^ is the usual Mandelstam invariant, which can be viewed as 
the hadronic center-of-mass energy. It is clear that Eq. lHTl) represents the fraction of 
energy that the hadrons transfer to the photon and, hence, < x < 1. 

According to the factorized expression of the QCD cross section Eq. (l44l) . the LO 

differential cross section is given by: 

'^^ 'x, Q^) = y c?xirfx2 [gj(xi,/i^)gi(x2,/i^) + gi(xi)gi(a;2,/i^))] (48) 



dQ 



d&i nv / N „o ,- / N nv / „o N Allot 



a^'{Q\x)QU{xiX2-x), a^'{Q\x) = —-x (49) 



Figure 1.3: Drell-Yan pair production. Here Q = M. 



where the functions qi{xj, fi'^){qi{xj, ^'^))aie the parton densities of a quark (or an 
anti-quark) of flavor i in the hadron j = 1, 2 at the scale /i^, a is the fine-structure 
constant and Qq. is the fraction of electronic charge of the quark gj. Now, if we define 
the dimensionaless cross section a{x, Q"^) as, 



a^''{x,Q')^^^^{x,Q'), (50) 



and use the identity. 



"1 

5{xiX2 — x)= I dz5{l — z)5{xiX2Z — x), (51) 



then Eqs.( H8fl9!) become: 

a^^(x, g^) = 

= dxidx2dz [qi{xi)qi{x2) + qi{xi)qi{x2)] Ql^Cqq{z) 6 {x 1X2 z - x 

i Jo 

r dxi r dx2 [^^(^3,^)^.(^3,2) + qi{xi)qi{x2))] Q^Cq 

■ Jx -^1 Jx/x\ -^2 



(52) 



(53) 



Cqq{z)=6{l-z), (54) 

where Cqq{z) is the LO Drell-Yan coefficient function. From Eq. fl52l) . we see that the 
new variable z that we have introduced is in general given by 



X 



X1X2 



(55) 



This means that at the partonic level, z can be viewed as the fraction of energy that 
the colliding partons transfer to the virtual photon. At LO it is clear that z = 1 
as can be explicitely seen from Eq.( |5^ . because there is no emission but the virtual 
photon. 



Beyond the LO the extra radiated partons in the final state can carry away some 
energy (so z < 1) and the gluon channel contributes. The NLO coefficient func- 
tions Cab{z) {a,b = q,g) recieves contributions that have infrared and ultraviolet. 
Infrared singularities cancel out (see e.g. |25J). The ultraviolet ones are reabsobed 
by renormalization of the bare parameters of the QCD Lagrangian, thus defining a 
renormalized strong coupling constant a^2 at an arbitraty renormalization scale /i^ 
(see section [L2l) . Collinear divergences are cut off by infrared physics. They can be 
absorbed multiplicatively in redeffnition of the parton densities [26], thus reabsorbing 
all dependence on soft physics in the parton distributions. The parton densities at a 
certain scale are determined by a reference process and their scale dependence is de- 
termined by the DGLAP equations (see section [L^ . However, there is an ambiguity 
on how to define the reference process, related to the fact that collinear divergences 
can always be factorized together with finite terms. The choice of these finite terms 
defines a factorization scheme. The most common factorization scheme is the MS 
scheme in which the collinear divergence (which is in = 4 — 2e dimensions a single 
pole 1/e) is factorized together the finite terms —■je + log47r, where ■je = 0.5772... 
is the Euler gamma. Now, in order to avoid the perturbative expansion to receive 
large contributions, the factorization and the renormalization scales are expected to 
be chosen of the same order of the scale of the process Q^. Here, for simplicity, we 
choose the factorization scale /i^ equal to the renormahzation scale /x^. We report the 
NLO Drell-Yan cross section (see e.g. [271 128] ): 



^ Jx -^l Jx/xx •''2 



X I [qi{xi, n^)qi{x2, /i^) + (1 ^ 2)] C^q yz, — , a,(/i^) 
+ [9{xi, /i^) (gi(x2, /i^) + gi(x2, /i^)) + (1 ^ 2)] Cqg I z, — , ^^(/i^ 



(56) 



where, in the MS scheme, 



Cqq [ ^^OisifJ-"^] 
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and 
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Figure 1.4: Deep inelastic electron-proton scattering 



where the "+" distribution is defined as follows: 

\zf{z)[g{z)]+^ f dz[f{z)~f{l)]g{z). (59) 
Jo 

Also for the case of the deep-inelastic scattering (DIS), we consider the simplest 
process in which a high energy electron scatters from a hadron exchanging with it a 
virtual photon 7* (see Figure fOl) : 

H{P) + e{k)-^e{k')+X{K), (60) 

where H is typically a proton with momentum P, e is the scattered electron and X 
is any collection of hadrons. The standard parametrization of DIS is done in terms 
of three relevant parameters: 

q2 ^ _q2^_{^k-k'f (61) 

y = o<y<i (62) 

- - ^-^■ = ^=(p + ,Vq^' ^''^ 

where in the last line we have neglected the proton mass. is the virtuality of the 
photon exchanged between the electron and the proton and y is the fraction of energy 
that the incoming electron transfer to the proton. The Bjorken scaling variable x has 
a simple physical interpretation: it is the fraction of longitudinal momentum of the 
LO incoming quark of the partonic subprocess. 

Indeed, the most general parametrization of the differential cross section is 
given by: 

da , „o X Ana^ ,^ so.„, (I — y) 



{x,Q ,y) 



dQ^' '-^ Q4 



[1 + (1 - yy]F,ix, + }__J^l^F,ix, Q') - 2xF,ix, Q')) 



(64) 



The functions -^1(2) are called structure functions and contains the information about 
the structure of the proton. In fact they are determined by the photon-proton sub- 
process in this way: 



Fii^^Q') = -r^Wj:{l*P) + aLil*P)] (65) 
F2{x,Q^) = 2xFiix,Q^)+FL{x,Q^) (66) 

Fl{x,Q') = -^az.(7*P), (67) 

where crs(7*-P) and aj]{'y*P) are the cross sections of the photon-proton process deter- 
mined summing over all the virtual photon polarization and over only the longitudinal 
one respectively. At LO 

Fiix,Q') = ^5^g;jg.(x,Q2) + g,(x,g2)] (68) 

i 

F2{x,Q^) = 2a;Fi(x,/i'), (69) 

where qi and qi are the parton densities. In general the structure functions should 
depend on both x and Q^, because these are the relevant kinematic variable of the 
photon-proton sub-process. Now we want, as we have done for the Drell-Yan case, 
rewrite the structure functions in terms of parton densities and of a coefficient function 
that can be computed in perturbative QCD. If we use the identity 

S{y -x)= ! dz (5(1 - z)5{yz - x), (70) 
Jo 

we have for the LO structure functions F2 and Fl. 

F2{x,Q^) = / dydz [qi{y,fi^) + qi{y,fi^)]QlCg{z)6{yz-x) (71) 

i Jo 

= ^ E - [^^(^' + i^^y^ /^')] QIC. (-) (72) 

j Jx y \y / 

C,{z) = 5{l-z) (73) 
Fi(x,Q2) = 0, (74) 

where Cq{z) is the LO DIS coefficient function for F2. From Eg. (1711) . we see that the 
variable z is in general given by 

z = -. (75) 

y 

This means that at the partonic level, z can be viewed as the longitudinal momentum 
of the incoming parton before it scatters with the virtual photon. At LO it is clear 
that z = 1 as can be explicitely seen from Eq. fl73p . because there is no extra emissions. 

Beyond the LO the extra radiated partons in the final state can carry some energy 
(so z < 1) and also the gluon channel contributes. At NLO we have ultraviolet, 
infrared and coUinear singularities. They must be regularized and treated as in the 



Drell-Yan case. We report the NLO structure functions (see e.g. [271 EH]) with the 
renormahzation scale equal to the factorization scale: 



(76) 



where, in the MS scheme, 
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(78) 



(79) 



is factorization scheme independent at the lowest non-vanishing order. 



1.4 The DGLAP equations 

The coefficient function and the parton densities depend on the factorization scale 
in such a way that the resulting hadronic cross section is //^-independent. The 
equations that fix the ^^-de'pendence of parton densities (the so called Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations) can be found imposing the ix^- 
independence of the DY cross section or of the DIS structure functions. For example. 



imposing this condition to the exphcit expression for the NLO F2 (see Eqs. fl76|77f78l) 
we find the LO DGLAP evolution equations for the quark parton densities: 



/i 



,dqi{x^ 
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y '''\y 
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where 



and 







3 L[l — z\+ 2 



with Nf the number of active fiavors. The functions Ppp}{z) are called LO splitting 
functions. They can be viewed as the probability per unit of ln(/i^/Q^) to find a 
parton p in a parton p'. The LO evolution equation for the gluon can be calculated 
from the LO splitting diagrams for a quark into another quark and a gluon and 
for a gluon into two gluons. Furthermore, we simplify the notation introducing the 
convolution product 0, defined in this way: 



(80) 

(81) 
(82) 

(83) 
(84) 
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We report here the full result for the DGLAP evolution equations: 



/i 



E 



Pg^g^iz, fj.^) Pg,g{z,fj.^) 

Pgg,{z,fl^) Pgg{z,fl^) 



i5) 



where can be also a quark or anti-quark and where the splitting functions Ppp' have 
the following perturbative expansion: 



q^qj\■^^M' J "ij^ qq 
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Pggiz,fJ^^ 



E 

fc=i 



fc=i 

477 



99 



(89) 
(90) 
(91) 



where Nf is the number of active flavors. Eq.( l86i) represents a system of 2Nf + 1 
integro-differential equations equations. The solution to this system however can be 
calculated analytically for a certain fixed-order. In fact it can be translated into a 
system of ordinary differential equation performing a Mellin transform: 



Fp(iV,/i2 



where p, p' = qi, qj, g and 



7//(iV,/i') 
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dzPppi{z,^^). 



(92) 
(93) 

(94) 



After that, these equations can be decoupled searching linear combinations of par- 
ton densities that depends on the independent splitting functions of Eqs. fl68ti77p and 
that diagonalize the system. For example, at LO, there are 4 independent splitting 



functions which are Pc 



V(0) 
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j(0) p(0) 
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Pqq' 1 Pqg' (glvcu lu Eqs.( l83|8^ respcctivcly) and 
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Nf 6{l-z). 



(96) 



At NLO there are 6 independent splitting functions which are for example P^q{z), 
Pgq{z), P^g{z), Pqg{z), Pgq{z) and Pgg{z). They are given in Ref. [27] pages 111 and 
112. The LO and the NLO solution to the DGLAP equations ( l93i) in Mellin space is 
computed in the next Section. 



1.5 NLO solution of the DGLAP evolution equa- 
tions 

In this Section, we want to solve the NLO DGLAP equations (Eq.f l93l) of section [L4l) : 



(97) 



P' 



where all the splitting functions defined here (and in the following) have the same 
perturbative expansion: 

T„.(iV.A.^) = ^-Cm + (^)'-rS(iV) + 0(al). (98) 

This is a system of 2Nf + 1 coupled equations with Nj the number of active flavors. 
At NLO, there are 6 independent splitting functions defined through the following 
equations: 



Igq. = Igq. = Igi (99) 
%.9 = 7,-,p = l/iV/7,s (100) 

Igilk ~ 7(?i(?fe = ^iklqq + 1 qq (101) 
Imk = Imk = ^iklqq + Iqq (102) 

Igg = Igg^ (103) 

where i, k are a fiavor index. We omit the dependence on and /x^ for brevity 
of notation. Note that beyond the NLO, there is one more independent splitting 
function. In fact in Eq. (11021) we should substitute 7^^ with 7;^^ which are different 
beyond the NLO [221 EDj- We note also that at LO 7^^ = 7^. = 7^^ = and hence at 
LO there are only 4 independent splitting functions. 

Now, we define the 2A^/ — 1 so called non-singlet (NS) combinations 

k 

^tNS)k = J2^qi ± q^) - k{qk ± qk); k = 2,...,Nf (104) 

i=l 

^Ins) = "^(qi-qi) (105) 

i=l 

and the 2 so called singlet (S) combinations: g and 

q(s) = J2^q^ + q^)■ (106) 



i=l 



With this definitions, from Eq. (l97|) and Eqs. fl99til03p . we find that for the non-singlet 
combinations 



2^^)k _ 

^ ~7h:^ = ^ ^(^s). (108) 



/^^^ = 7^^f..). (107) 
9/i 



where 



(109) 



For the 2 remaining singlet combinations, we find in the same way that 



9{s) 
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Iqq Iqg 



9 



;iio) 



where 

Iqq = 7^ + 7P5, IPS = '^Nf-f^g. (Ill) 

The NLO Melhn sphtting functions can be found in Ref.[3I] written in terms 
of harmonic sums. In many cases, however, their analytic continuation to all the 
complex plane is useful (see e.g. [32l l33] ) For the NNLO solution of the DGLAP 
equations and the NNLO splitting functions we refer to [29],[30j. The techniques for 
the analytic continuations of the NNLO splitting functions can be found in Ref . [34] . 

For the NS combinations, the solution is easy to obtain. Indeed, making the 
change of variable 



where P{as) is the /5 function defined in section [T72| we get: 
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where we have omitted the A^-dependence of the splitting functions for brevity. For 
the S combinations, some linear algebra is needed. We, first, define the singlet vector 
and the splitting matrix: 



Qs 
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Iqq Iqg 
739 Igg 



(115) 



Using the NLO splitting matrix and the change of variable Eq. flll2l) . we find im- 
meditely the formal solution, which is 



where 
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(117) 



The two matrices Rq and -Ri in Eq.f lll6p cannot be diagonalized simultaneously, as 
they do not commute. Hence, in order to extract the NLO solution from Eq. flll6p . 
we use the following Ansatz: 



as(/io) 



118) 



where the matrix U has the perturbative expansion: 



(119) 



The condition that the matrix Ui should satisfy can be easily obtained imposing that 
the derivative with respect to as(/i^) of Eg. (11 161) and of Eg. (II 181) are egual at NLO. 
Thus, we get 

[U,,Ro] = U, + R,. (120) 



We write Rq in terms of its 2 eigenvalues 
R^ 



1 

Wo 



(121) 



and of the 2 corresponding eigenspaces projectors P±: 

Ro = R+P+ + R-P^. 



(122) 



The explicit expression for the projectors can be obtained using the completeness 
relation P. + P_ = 1. We find 



P4 



1 



■[Ro-R^]. 



Now, writing the matrices Ui and Ri in terms of these projectors 

f/i = P^UiP^ + P^UiP+ + P+UiP^ + P+UiP+ 
Ri = P.RiP- + P-RiP+ + P+RiP- + P+RiP^ 



(123) 



(124) 
(125) 



substituting them and Eg. (11221) into Eg. (11201) and comparing each matrix element, 
we find 



[P-RiP- + P+RiP^ 



P+RiP- ^ P-RiP, 



R- -R+ -1 R+-R- 



(126) 



Thanks to this result, we can now write the NLO solution of the singlet doublet in 
a form which is useful for practical calculations. Indeed, if we substitute Eg. (11261) in 
Eg.([Il8]), we get (at NLO) 



9s (/i^ 
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(127) 



After the evolution of the NS and S combinations has been performed from a certain 
scale /ig to the scale /x^, we need to return to the parton distributions for all the quarks 
but the gluon. These are obtained straightforwardly with the following relations 



qk + qk = -j^qs - plNS)k + -(—J)'ltNS)v k = l,...,Nf (128) 
■f i=k+i ^ ' 

1 1 ^ \ 

qk-qk = -j^qlNS) - piNS)k + Yl ^u_if iNS)v k = l,...,Nf. (129) 

f i=k+i ^ ' 

However, Eqs . f l 1 1 3ll 1 1 4| 1 2 71) represent the NLO solution of the DGLAP equations 
Eq.( l97I) in the case when the number of active flavors Nf has been kept fixed. This 
is the so called fixed flavor scheme solution. If we want to take into account the 
thresholds of the heavy quark flavors, we can evolve up the NS a S combinations from 
the scale /Iq (with a certain number Nf of active flavors) to the scale of production of 
a new flavor. Then, we can take the result of this evolution as the starting point of a 
second evolution (with A'^^ + 1 active flavors this time) above the production scale of 
the new flavor, assuming that the new flavor vanishes at threshold. This is the most 
simple way to generate dynamically a new flavor. 

Finally, we note that the procedure outlined in this appendix can be easily gen- 
eralized beyond the NLO order [33]. Furthermore, in many cases, it is interesting 
to study the dependence on the renormalization scale, in order to estimate the the- 
oretical error of the evolution. Here the renormalization scale //^ has been chosen 
equal to the factorization one //^ for simplicity. To restore the implicit /i^-dependence 
in parton densities, we need only to rewrite the running coupling constant ^^(/i^) 
in terms of /i^ (see Eq. (l41l) in section \UI\j in the splitting functions. Making this 
substitution, we have that the perturbative expansion of a generic splitting function 
Eq.f ll30p becomes 

7„,(/V. ,\ k') = ^^^^oiJkA') + (^^) ' (I'^m + /3o7S! 'n + 0(4), 

(130) 

where k' = Hence, in Eqs. (11 131)1 14fl27p . we should perform the following 

substitutions: 

^(1) ^ +/3o7(o)lnfc', asifi') ^ «.(fcV), a.(/i^) ^ «.(A;Vo) (131) 
and use k'fi^ as reference scale for new flavors production. 



Chapter 2 

High order QCD and resummation 



2.1 When is NLO not enough? 

In section 11.31 we have discussed briefly the analytic NLO calculation of the full 
inclusive DIS and DY cross sections. However, in many cases, the NLO pQCD com- 
putation turns out not to be enough. This is, for example, often the case at LHC 
where the Higgs boson production has to be distinguished from the background. A 
computation beyond the NLO is needed also when the NLO corrections are large 
and higher-order calculation permit us to test the convergence of the perturbative 
expansion. In figure 12.11 the total cross section of the production of the Higgs boson 
at LHC [36] is plotted and we note convergence in going from LO to NLO and to 
NNLO. 

This can also happen when a new parton level subprocess first appear at NLO. 
This is the case for example for the rapidity DY distributions at Tevatron (shown in 
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Figure 2.1: Total cross section for the Higgs boson production at LHC at (from bottom to top) 
at LO, NLO, NNLO in the gluon fusion channel Hg. 
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Figure 2.2: DY rapidity distribution for proton anti-proton collisions at Tevatron at (from bottom 
to top) LO, NLO, NNLO, together with the CDF data [37]. 
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Figure 2.3: DY rapidity distribution for proton proton collisions at fixed-target experiment 
E866/NuSea at (from bottom to top) LO, NLO, NNLO, together with the data [HIIH]. 



figure EH]) and at the fixed-target experiment E866/NuSea (shown in figure . The 
agreement with the data of figure EH] has represented an important test of the NNLO 
splitting functions [291 EH]- We note also that going from the LO to the NNLO the 
factorization scale dependence is significantly reduced. 

Calculations beyond the NLO can be important also in processes which involve 
large logaritms when different significant scales appear. In these cases, these large 
logarithms should be resummed and this is the topic of this thesis. A first example of 
these large logarithms has appeared in section [L3l In fact, from Eqs. fl57|77p of section 
IL31 we see that there are contributions that become large when z ^ 1 from the 
quark- ant iquark channel in the DY case and from the quark channel for the structure 



function F2 in the DIS case. These are the terms proportional to 



'log(l - z) 
1 - z 



a. 



The terms of the type of Eq. ([T]) arise from the infrared cancellation between virtual 
and real emissions. It can be shown that enhanced contributions of the same type 
arise at all orders. In fact, at order there are contributions proportional to 

[DEIE]: 

'log'"(l - z 



a. 



m < 2n 



(2) 



These terms become important in the limit z —>■ 1 spoiling the validity of the per- 
turbative fixed-order QCD expansion and, thus, should be resummed to all-orders of 
pQCD. 

The limit z —>■ 1 corresponds in general to the kinematic boundary where emitted 
partons are all soft (as it happens in the DY case) or coUinear (as in the DIS and the 
prompt photon case as we shall see in section 14. 3p . 

In fact, in the DY case, if we consider a contribution to the coefficient function with 
n radiated extra partons with momenta ki, . . . ,kn, the squaring of four-momentum 
conservation {pi + P2 = Q + ki + ■ ■ ■ + kn) implies 



XiX2S{l — z) = ki ■ kj + 2 Q ■ ki 



(3) 



i=l 



J2 k^A^ - + 2 J2 k'WQ' + \Q\'- \Q\ cosed, (4) 



i=l 



where 6ij is the angle between ki and kj and 6i is the angle between ki and Q. 
Since all the terms in the first sum of Eq.(jl]) are positive semi-definite and the terms 



(y + IQp — \Q\ cosOi) in the second sum are positive for all possible values of 6i, 
we have that the limit 2; = 1 is achieve only for kf = for all i. This means that 
when z approaches 1 all the emitted partons in the Drell-Yan processes are soft and 
that we have reached the threshold for the production of a virtual photon or a real 
vector boson. 

In the DIS case, at the partonic level, we have 



p + q = ki^ \-kn + kn+i, 

where kn+i is the LO outgoing parton. If we square this last equation, we get 



(5) 



(6) 



where 9ij is the angle between ki and kj. Eq.(l6]) tells us that in the z — > 1 limit, there 
can be not only soft radiated partons in the final state, but there can be also a set 



LHC 



NNLL 
NWT.O 



30 




Fi gure 2.4: Total cross section for the Higgs boson production at LHC at (from bottom to top) 
at NNLO and NNLO improved with NNLL resummation in the ghion fusion channel [39 . 
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Figure 2.5: Total cross section for the transverse momentum Higgs boson production at LHC at 
LO and LO improved with NLL resummation in the gluon fusion channel [ID] . 



of partons collinear to each other. However, in Section 13.21 we will show with a more 
detailed analysis of the DIS kinematics and phase space that the collinear partons 
are also soft in the z —>■ limit for the deep-inelastic process. 

An example of the impact of resummation can be seen in figure 12.41 There, the 
total cross section for the Higgs production at LHC is plotted at NNLO with its 
NNLL resummation improvement improvement [39]. The scale uncertainty reduced 
to about 15% at NNLO is further reduced to 10% by the NNLL resummation. The 
resummed large logaritms in this case are of the class of the DY-like soft emissions. 

Resummation of another class of large logarithms plays a crucial role in transverse 
momentum distributions. Indeed, in figure (231 we observe that resummation changes 
substantially the behavior of the cross section for the production of the Higgs boson 
at small transverse momentum. In these case the large logarithms of q^/Mfj with 
Mh the Higgs mass are resummed. 

The state of art of QCD predictions for Higgs boson production at LHC is reported 
in figure 12.61 as it was summarized by Laura Reina at the CTEQ summer shool 2006 
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Figure 2.6: State of art of QCD predictions for Higgs boson production at hadron colliders. 



on QCD analysis and phenomenology, where also the Monte Carlo event generators 
are indicated. 

Furthermore, at LHC, multi-particles/jet production will be the inescapable back- 
ground to Higgs searches and searches for new physics. This means that we should 
have a precise knowledge of the QCD background. As seen previously, we know many 
QCD processes up to the NNLO. However, we have at the moment limited NLO 
knowledge of some important final states that will constitute background. They are 

W/Z +iets (2j) 
WW/ZZ/WZ + jets (Oj) 
WWW/ZZZ/WZZ + jets (Oj) 

gg+ jets (oj) 

7+ jets (Ij) 
77 + jets 
Z77 + jets, 

where in parenthesis is indicated the NLO knowledge. 

Finally, we also note that in higher order contributions to the splitting functions 
{Pgg,Pgg foi examp^), it can be shown that there can appear also terms proportional 

to 

<ln™-; m<n. (7) 

z 

These contributions spoil the convergence when 2; — and, in order to have reliable 
predictions, must be resummed. The inclusion of the terms with m = n defines a LL^ 
resummation, the inclusion of also the terms with m = n — 1 defines a NLLz resum- 
mation. This resummation is realized by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 
equation. Anyway, in this thesis, we will not concentrate on this resummation. We 



will give a briefly description of the various techniques to resum the large soft logs 
giving attention to the renormalization group approach and studying in detail its 
applications. 

2.2 The renormalization group approach to resum- 
mation 

The aim of resummation is to include all the logarithmic enhanced terms of the form 
of Eq. ([2]) of Section 12.11 with a certain hierarchy of logarithms that we shall deflne in 
the current section. 

From Eqs.( l56|76l) of Section [T73] we see that the QCD cross section (up to dimen- 
sional overall factors) can in general be written as a convolution of the parto densities 
F^*(xj, /i^) and of the dimensionless partonic cross section, i.e. the coefficient function 
C(z,QV/.^a.(/i2)): 

aDY(x,Q^) = J2 [i^i'H/i') ® ® C'.^WV/^',".!/^'))] (x), (8) 

a,b 

for the DY case; and 

anisix, Q') = J2 [Fail^') ® CaiQV^', «.(/^'))] (^), (9) 

a 

for the DIS case. The convolution product ® has been defined in Eq. (!85|) of Section 
11.41 Performing the Mellin transformation 

(t{N,Q^)= [ dxx^-^a{x,Q^) (10) 
Jo 

we turn the convolution products of Eqs. llHllQl) into ordinary products: 
aDY(iV,Q') = 5^a,fe(iV,g2) 

a,b 

= ^Fi^niv,/i')i^f^(A^,/i')a.(iv,QV/i',«.(/i')), (11) 

a,b 

^DIs(iV,g') = $^a,(iV,Q2)=^F^^(iV,^2)^^^^^g2/^2^^^(^2))^ ^^2) 



where 




Fa^^iN,^^') = /'dxx^-iFi^^(x,/i2), (14) 

and where the second index in brackets (6) is involved only when there are two hadrons 
in the initial state as the the DY case. 



The large logs of 1 — 2; of Eq. ([2]) in Section 12.11 are mapped to the large logs of 
by the Mellin transform. This fact and the relations between the large logs of 1 — z 
and the large logs of N are shown in detail in Appendix 1X1 

Whereas the cross section a{N, Q"^) is clearly /i^-independent, this is not the case 
for each contribution aa{b){N,Q'^). However, the iJ? dependence of each contribution 
to the sum over a, (6) in Eqs.( |TT|T2l) is proportional to the off-diagonal anomalous 
dimensions 7^^ and 7^5. In the large N limit, these are suppressed by a power of 
in comparison to 7gc, and 7^^, or, equivalently, the corresponding splitting functions 
are suppressed by a factor of 1 — 2; in the large z limit (see for example Eqs.( 



in Section ri.41) . Hence, in the large limit each parton subprocess can be treated 
independently, specifically, each Ca(b) is separately renormalization-group invariant. 
Because we are interested in the behaviour of Ca(b)(A^, Q^Z/x^, as(yU^)) in the limit 
A^ — > cxo we can treat each subprocess independently. 

Because resummation takes the form of an exponentiation, we define a so-called 
physical anomalous dimension defined implicitly through the equation 



The physical anomalous dimensions 7a(b) Eq. ( JTSl) is independent of factorization scale, 
and it is related to the diagonal standard anomalous dimension 7^^, defined by 

^.dF^iE^ = 7er(A^,«.(/^^))^c(iV,/.2), (16) 

according to 

la{b){N,as{Q )) = dh^^ = laaiN,as{Q )) (17) 

+7ibb)iN,as{Q )) + . (18) 

We recall that both the standard anomalous dimensions (Altarelli-Parisi splitting 
functions) and the coefficient function are computable in perturbation theory. Hence, 
the physical anomalous dimensions differs from the standard anomalous dimensions 
only beyond the LO in as as can be seen directly from Eq.( fT8l) . In terms of the 
physical anomalous dimensions, the cross section can be written as 

a{N,Q') = Y,Ka(b){N-QlQ^)aa(b){N,Ql) (19) 

a,(fe) 

= ^exp [Eaib){N;QlQ^)\ aaib){N,Ql), (20) 

a,{b) 



Q' " ' = 7aw(iV, asiQ')) a,(,)(Ar, Q'). (15) 



where 



Ea{b) 



{N-QlQ') = I _7,(,)(Ar,a,(fc2)) (21) 







k 



+ lnQ,)(Ar,l,a,(Q^))-lnQfe)(Ar,l,a,(Q^)). (22) 



We now concentrate on the single subprocess with incoming partons a, [b). Re- 
summation of the large logs of N in the cross section is obtained performing their 
resummation in the physical anomalous dimension: 



This shows how in general the large logs of N can be exponentiated. For the DY case 
only the quark-anti-quark channel should be resummed and in the DIS case only the 
quark one. This is a consequence of the fact that the off-diagonal splitting functions 
are suppressed in the large N limit as discussed before. 

The accuracy of resummation here depends on the accuracy at which the physical 
anomalous dimension 7 is computed. We say that the physical anomalous dimension 
is resummed at the next*^" ^-to-leading- logarithmic accuracy {N^~^LL) when all the 
contributions of the form 



are included in its determination. The goal of resummation is to determine the re- 
summed physical anomalous dimension from at most a finite fixed-order computation 
of it. Clearly, once the resummed physical anomalous dimension is determined, it 
can then predict the leading, next-to- leading... logarithmic contributions to the cross 
section at all orders. 

Here, in this Section, we shall only outline the key ideas of the renormalization 
group approach to resummation. However, throughout this thesis we shall show in 
detail how this method works and how the resummed physical anomalous dimensions 
obtained with such an approach can be fully determined expanding it to a certain 
finite fixed-order and comparing this expansion with Eq.l fTSl) obtained from explicit 
computations. 

The renormalization group approach to resummation is essentially divided in two 
steps. The first is to analyze the generic phase space measure in = 4 — 2e dimensions 
thus finding where the large logs are originated in the coefficient function and in the 
physical anomalous dimension. The second consists in resumming them imposing the 
renormalization group invariance of the physical anomalous dimension. 

So, let's consider a generic phase space measure d(j)n for the emission of n massless 
partons with momenta pi, . . . ,Pn- In Appendix |Bl we show that this phase space 
can be decomposed in terms of two-body phase space. Roughly speaking, the phase 
spaces measure for the emission of n partons can be viewed as the emission of two 
partons (one with momentum p„ and the other with momentum Pn = pi + ■ ■ ■ +n-i 
and invariant mass times the phase space measure where the momentum P„ is 
incoming and the momenta pi, . . . ,pn are outgoing. The price to pay for this is the 
introduction of an integration over the invariant mass P^. Then using recursively 
this procedure, we obtain that the n-body phase space measure is decomposed in 
n — 1 two-body phase space measures. This means that we have reduced the study 
of the soft emission of n-body phase space measure to the study of the soft emission 
of two-body phase spaces. 




(23) 




n = 0, . . . , A; — 1 



(24) 



The two-body phase space with an incoming momentum P and two outgoing 
momenta Q and p in d = 4 — 2e dimensions (derived exphcitly in Eq.( fT6ll of Section 
|B|) is given by 

l-2e 



.XP;Q,p) = Nie)iP 



2\-e 



1-91 

p2 



iV(e) 



1 



(25) 



2(47r)2-2e' 

where dVld-i is the sohd angle in d — 1 dimensions. For definiteness, we can think 
that this is the phase space measure for a single DIS-like or DY-like soft emission 
with with momentum p. Thus, we have 



DIS-like emission: 
DY-like emission: 



oc (1 - 



Q' = 



1 



p2 



oc d 



p2 



(26) 
(27) 



where z^^^g, z^^ are close to one for a soft emission. Hence, we have that the two-body 
phase space measure for a single soft emission contributes with a factor (1 — z)^"^ 
with a = 1 for a DIS-like emission and a = 2 for DY-like emission. In the case of 
the prompt-photon process, we will see in Chapter H] that there are both types of 
emission. The large logs of 1 — z are originated by the interference with the infrared 
poles in e = in the square modulus amplitude in the e — > limit. For example 



e 



- - ln(l - zY + 0(e) 



(2^ 



for the case of interference with a pole of order 1. 

Now, since each factor of (1 — z)~°-^ that comes from the phase space measure 
is associated to a single real emission then it will appear in the coefficient function 
together with a power of the bare strong coupling constant a^. In d-dimensions, the 
coupling constant is dimensionful, and thus on dimensional grounds each emission is 
accompanied by a factor 

ao [Q\l-zY]-\ (29) 

where is now the typical perturbative scale of a certain process. Upon Mellin 
transformation, this becomes 



(30) 



Furthermore an analysis of the structure of diagrams shows that in the soft (large 
A^) limit, all dependence on appears through the variable Q'^/N"' also in the am- 
plitude. Finally, a renormalization group argument shows that all this dependence 
can be reabsorbed in the running of the strong coupling. Indeed, the first order 
renormalization of the bare coupling constant at the renormalization scale /i 



(31) 



and the renormalization group invariance of the physical anomalous dimension imply 
that 

'QT . 2.1 Q' ' ' 



ao 



AT" 



+ 0{ai) 



(32) 



where ao is the bare couphng, ^^(/i^) the renormahzed couphng and the higher order 
terms contain divergences which cancel those in the cross section. Following this 
line of argument one may show that the finite expression of the renormahzed cross 
section in terms of the renormahzed coupling is a function of as(Q^) and a<j(Q^/iV") 
with numerical coefficients, up to 0(1/ N) corrections. We shall see this in detail in 
Chapter [3l 



2.3 Alternative approaches 

The exponentiation of the large soft logs and their resummation has been demon- 
strated in QCD with the eikonal approximation [2J or thanks to strong non-standard 
factorization properties of the cross section in the soft limit [1]. Recently, also the 
effective field theoretic (EFT) approach has been applied to QCD resummation in 
Refs. [9l [To] for DIS and DY and in Ref. [H] for the B meson decay B Xg'y. In this 
Section, we shall only give a brief description of these alternative approaches to the 
resummation of the large perturbative logarithms. 



2.3.1 Eikonal approach 

We first consider the simpler case of QED. In QED the exponentiation of the large soft 
logs has been proved thanks to the eikonal approximation in Ref. [7]. We report the 
basic steps of the proof for the QED case and a brief description of the generalization 
to the QCD case. 

Consider a final fermion line with momentum p' of a generic QED Feynman dia- 
gram. We attach n soft photons to this fermion line with momenta ki, . . . ,kn- For 
the moment we do not care whether these are external photons, virtual photons con- 
nected to each other, or virtual photons connected to vertices on other fermion lines. 
The amplitude for such a diagram has the following structure in the soft limit: 

u{p')HeY^)7^Hein ^ , .f',, ■ ■ ■ (-^e7/in) ^ , f ^ , . 2.M,,(33) 
2p' ■ fci 2p' ■ (fci + k2) 2p' ■ (fci H h kn) 

where e = — |e| is the electron charge and iAih is the amplitude of the hard part of 
the process without the final fermion line we are considering. We note that here we 
have neglected the electron mass. Then we can push the factors of to the left and 
use the Dirac equation u{p')f = 0: 

uip'h'^'fY^f ■ ■ ■ = u{p'W"'l^^f ■ ■ ■ = u{p'W^"^P'^^ ■ ■ ■ V^"- (34) 
Thus Eq. (!33l) becomes 



Still working with only a final fermion line, we must now sum over all possible order- 
ings of momenta ki, . . . ,kn- There are n\ different diagrams to sum, corresponding to 



the n\ permutations of the n photon momenta. Let P denote one such permutation, 
so that P{i) is the number between 1 and n that i is taken to. Now, using the identity 



y _^ i i _ _^ ^_ .ggX 

p P ■ kp(i) p ■ (A;p(i) + A;p(2)) p ■ (A;p(i)) H h kp^n)) P ■ ki p ■ k^ 

the sum over all the permutations of the photons of Eq. flH^ is: 

At this point, we consider an initial fermion line with momentum p. In this case 
the photon momenta in the denominators of the fermion propagators have an opposite 
sign. Therefore, if we sum over all the diagrams containing a total of n soft photons, 
connected in any possible order to an arbitrary number of initial and final fermion 
lines, Eq.(1H71) becomes: 

«"'^»nEr^. (38) 

-I . Pi rvr 

where iM.Q is the full amplitude of the hard part of the process and where the index r 
runs over the radiated photons and the index j runs over the initial and final fermion 
lines with 

{1 for a final fermion line , 
— 1 for an initial fermion line 

If only a real soft photon is radiated, we must multiply by its polarization vector, 
sum over polarizations, and integrate the squared matrix element over the photon 
phase space. In the Feynman gauge this gives a factor 




TTT^eM) (40) 



in the final cross section. If n real photons are emitted, we get n such Y factors 
Eq. lHOl) . and also a symmetry factor 1/n! since there are n identical bosons in the 
final state. The cross section resummed for the emission of any number of soft photons 
is therefore 

^-(, _ /) = ^f)J2 = ao{^ - /)e^, (41) 

n=0 

where ao{i — > /) is the cross section for the hard process without extra soft emissions. 
This result shows that all the possible soft real emissions exponentiate and that only 
the single emission contributes to the exponent. However, this is not the end of 
the story, because the exponent Y Eq.l HOl) is infrared divergent. Indeed, to obtain a 
reliable finite result, we must include also loop corrections to all orders. For a detailed 
analysis about the inclusion of loops see for example Ref.^2]. Here, we just give the 
final result which reads: 

^-(, _ /) = ^,(, ^ /)e-'^\ (42) 



where a^^^ is the cross section relative to the single soft emission from the hard process. 
Clearly, the accuracy of this resummation formula for soft photon emission Eq. ( H2l) 
depends on the accuracy at which the exponent for the single emission is computed. 

In Ref.p] the exponentiation of the soft emissions, here outlined for QED, is gen- 
eralized to the QCD case. Differently from QED, QCD is a non-abelian gauge theory 
and this implies that this generalization is highly non-trivial. Indeed, the gluons 
can interact with each other. This fact makes the exponentiation mechanism much 
more difficult since the three gluon vertex color factor is different from that of the 
quark-gluon vertex. In order to exponentiate the single emission cross section (as it 
happens in QED), one should prove that these gluon correlations cancels out order by 
order in perturbation theory. This is shown for example in Ref . |13] . According to this 
result, it has been shown in Ref. [2] how the the exponentiation of soft emission works 
in QCD resummation. We report here the result for the NLL resummed coefficient 
function in Mellin space for inclusive DIS and DY processes in the MS scheme in a 
compact form: 



CNLL(iV,Q7/^ ,«s(/i')) 



exp < a 



dx 



X 



N-l 



X 



+B^-){a,{Q\l-xY)) 



(43) 



where 



A{as 



AiUs + A2al + 



(44) 
(45) 



with 

IT 



A, 



27r2 



[67 



-Nf 
9 ^ 



Bl 



(a) 



(2 - a)3C^ 
A-n 



(46) 



Here a = 1 for the DIS structure functio F2 and a = 2 for the DY case. How this 
result is strictly connected to the resummed results that can be obtained with the 
renormalization group approach will be discussed in Section 13.41 



2.3.2 Resummation from strong factorization properties 

This is the approach of Ref.pj. Also in this approach the results given in Eq. fHHl) 
of Section 12.3.11 are recovered. Here we give only a rough description of this method 
based on strong factorization properties of the QCD cross section. 

It is essentially assumed that at the boundary of the phase space, the cross section 
is factorized in a hard and in a soft part and eventually in an other factor associated 
to final collinear jets as in the DIS case where there is an outgoing emitting quark. 
The final result is then obtained exponentiating the soft and collinear factors. This 
is done solving their evolution equations. 

In Ref. [1] it is shown that the semi-inclusive cross section can be factorized in three 
factors relative to the three different regions in the momentum space of the process: 



the off-shell partons that participate to the partonic hard process, the collinear and 
soft on-shell radiated partons. The cross section is given by 

. N rj f Pi P2 y.\ f dwi dw2 dws f pi ■ Ci f QY 
a{w) = H[—,—Xi] ^1 ,wi[ — ] 

J2 {^~^^ ^2 ^ 5* {^s^^ Ci^ 5{w - wi - W2 - W3), (47) 

where a is the number of hadrons in the initial state, /x is the factorization scale and 
(i are gauge- fixing parameters; the integration variables Wi, W2 and W3 are referred 
to the two collinear jets and to soft radiation respectively. Each factor of Eq. fHTl) is 
evaluated at the typical scale of the momentum space region which is associated to. 
The delta function imposes that 

f 1 — XBi, for the DIS case 
w = w, + W2 + W3 = ^^_ Q2)s for the DY case ^^^^ 



The convolution of Eq. fl47j) is turned into an ordinary product performing the Mellin 
transform: 



00 



Each factor H, Ji, S satisfy the following evolution equations 



/i^— Ini/ = -7^(a,(/x2)), (50) 
d 



/^^^In^ = -7s(as(/i^)), (51) 
/x^AinJ, = -^j^^a.ifi')), (52) 

where the physical anomalous dimensions 7H(as), 7j,;(«s) and 75(0^) are calcula- 
ble in perturbation theory and must satisfy, according to the renormalization group 
invariance of the cross section the relation 

2 

Inias) +7sias) + ^l.hi<^s) = 0. (53) 

i=l 

Solving Eqs. ( I50ll52l) and imposing renormalization group invariance Eq.( l53il . the re- 
summed section can be written in the form 

.(iV) ^ exp{z.,«,g^)).Z..(«.(|;))-^^^fln(f).(«^ 

(54) 



-2 



' 'Q\ ..2^^ ..2 



where the functions A{as), B{as), Di{as) are determined in terms of the anomalous 
dimensions and the beta function. Finally, it can be shown that this result can be 
casted in the form of Eq. (H3i) in Section 12.3.11 for the resummed coefficient function. 



2.3.3 Resummation from Effective Field Theory 



This is the approach of Refs.[9l ^U\. This EFT methodology to resum threshold 
logarithms is made concrete due to the recently developed "soft collinear effective 
theory" (SCET) [3H |3Sl HHlllZI • The SCET describes interactions between soft and 
collinear partons. 

The starting point (considering the DY case as an example) is the coUinearly 
factorized inclusive cross section in Mellin space: 



(55) 



where cxo is the born level cross section and Fi{N,^'^) are the parton densities at 
the factorization scale /z. Here, the basic idea is to write the coefficient function 
C{N, Q"^ / ji^, asi^"^)) in terms of an intermediate scale /x/: 



C{N, QV/^', «.(/^')) = C{N, gV/i?, a^ifij)) exp 



-2 



^2 7,f (iV,«.(fc2)) 



(56) 



and then to compute C(A^, ^^(/ij)) with the "soft collinear effective theory" 

with the intermediate scale equal to the typical scale of the soft-collinear emission, 
i.e. fii = Q^N^ in the DY case and /if = Q^N in the DIS case. In SCET, 
C{N,Q'^/ fj,j,as{^'j)) has the following general structure: 



C ( A^, %,a^(/ij; 



c 



2 1 '^sif^Ij 



M{N,as{^A)), 



(57) 



Here C^Q"^/ fi'j,as{f^'j)) is the effective coupling that matches the full QCD theory 
currents 

Jqcd = CiQ^i^j, a,{fij))J,ff{i2j). (58) 

We note that the effective coupling C contains the perturbative contribution between 
the scale and fij and Jeff{iJ>j) contains the soft and collinear contributions below 
the scale fij. Then Ai{N,as{fij) is the matching coefficient that guarantees that the 
EFT used generates the full QCD results in the appropriate kinematical limit. In 
SCET the matching coefficient A4 can be computed perturbativly and is free of any 
logarithms. 

The effective coupling C satisfy to a certain evolution equation 



/i 



_d_ 



-7i(a,(/i^)) 



(59) 



where the physical anomalous dimension 71 is computable perturbativly in SCET. 
Finally, solving the evolution equation Eq. lj^Hjl one finds that Eq. (IFr|) becomes 



C(iV,^,«,(/if) 



C'(l,a.(Q')) 



^Ei{Q^/li'j,alphas{Q^)) ^ 



xA^(iV,a,(Q2))e^2(QVM?,-4Q^))^ 



(60) 



where 



^1 

with P{as) the beta function of Eq. (H3l) of Section [L2l C(l, as(Q^)) contains the non- 
logarithmic contribution of the purely virtual diagrams and the exponent Ei contains 
all the logarithms originating from the same type of diagrams. E2 encodes all the 
contributions due to the running of the coupling constant between the scale fij and 
Q^. All the logarithms appear only in the exponents (see Eqs. ( 156)1601) ) and the term 
1(7(1, as{Q'^))\'^M{N, asiQ"^)) is free of any large logarithms. 

The various approaches can be related one to the other according to factorization 
properties of the QCD cross section in the soft limit. In this way, it is possible to show 
that all the approaches are equivalent except for the renormalization group approach 
that produces correct but less restrictive results. 



1 dif 
4 7^2 ^ 



7i(as(^^)), 



dk 



dlnMjN, a,(P)) 
dlnas{k'^) 



(61) 
(62) 



Chapter 3 



Renormalization group 
resummation of DIS and DY 



In this chapter, we analyze in detail how the renormalization group approach to 
resummation works in the case of the all-inclusive Drell-Yan (DY) and deeep inelastic 
scattering (DIS). We recall that this is done only for the quark-anti-quark channel in 
the DY case and only for the quark channel in the DIS C3iSG clS discussed in Section 
12.21 First, we determine the N dependence of the regularized coefficient function 
in the large- limit. Then we show that, given this form of the A^-dependence 
of the regularized cross section, renormalization group invariance fixes the all-order 
dependence of the physical anomalous dimension in such a way that the infinite 
class of leading, next-to-leading etc. resummations can be found in terms of fixed 
order computations. This approach will lead us to resummation formulae valid to all 
logrithmic orders. 

3.1 Kinematics of inclusive DIS in the soft limit 

In the case of deep-inelastic scattering, the relevant parton subprocesses are: 



where Q is a quark or an anti-quark, ^ is a gluon and X is any collection of quarks 
and gluons. We are interested in the most singular parts in the limit z — > 1 of the 
phase space and of the amphtude for the generic processes Eqs. (ITpll . We treat first 
the tree level processes and then we will introduce the loops. 

Using Eq. (1121) in Appendix [B] with m = 1 recursively, we can express the phase 
space for a generic process in terms of two-body phase space integrals. For the DIS 
processes Eqs.( IT|2l) with n extra emissions {K = ki + ■ ■ ■ + kn) we have 



Qip') + X{K) 



(1) 

(2) 




d(j)n+i{p + q;ki,..., kn,p') 

/ -7r^d(j)2{p + q; kn, Pn)d(f)n{Pn] ki, . . . , kn-i,p) 
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MS 



27r 



^2iPn; kn^l, Pn~l)d(f)n-l{Pn-i; h, ■ ■ ■ , kn-2,P') 



ml 

27r 



^2(-fn+l j Pn 



M. 



27r 



Pn-l) 



X ■ ■ ■ X 



—^6/02(^3; k2, P2)dMP2; kl, Pi), 



(3) 



where we have defined = Pn, Pi = p', Pn+i = p + q and 

(1 - z) 



z 2p ■ q 

Now, according to Eq. (fTB|) in Appendix [B] we have for each two-body phase space 



(4) 



dh{P^+l■,ki,P^) = N{e){M, 



-2e 



where 



N{e) 



1 - 



Ml, 



2 \ l-2e 



2(4vr; 



2-2e 



dVLi\ i = 1, . . . n, (5) 



(6) 



and Qi is the sohd angle in the center-of-mass frame of Pj+i- We perform the change 
of variables 

= Tr~'^ Mi = sZn...Zi; i = 2,...,n. (7) 

From the fact that Mf_^'^^ > Mf (we have one more real particle in Pj+i than in Pj), 
it follows that 

< z,- < 1. (8) 



From Eq.© we get 



dMldMl_^ ■ ■ ■ dMl = det 



ml 

dzj 



dZjidZn—l ' ' ' dz2 



where 



det 



dzA 



dz„. dzo 



ra-l n-2 n-3 



Furthermore, in these new variables the two-body phase space becomes 
dMPi+u ki, Pi) = N{e)s~'{znZn.i ■ ■ ■ Zi+,)~'{1 - Zif-^'d^i. 
Substituting Eqs. lpfTOllTTl) into the generic phase space Eq.(j3]), we finally get 



(9) 



(10) 



(11) 



d(t)n+i{p + q;ki,..., kn,p') = 27r 



2tt 



n—l—ne 



dQn ■ ■ ■ dQi 



X / dZnZ^r'^-^-'^'i^-ZnY 



-2e 



dZ2Z2 ^(1 - Z2) 



l~2e 



(12) 



The dependence of the phase space on 1 — 2: comes entirely from the prefactor of 
^n-i-ne according to Eq.(jll). Indeed the dependence on z and has been entirely 
removed from the integration range thanks to the change of variables of Eq.([7]). 

Now, the amplitude whose square modulus is integrated with the phase space 
Eq.( fT2l) is in general a function: 

Ql,...,Qn). (13) 

The number of independent variables for a process with 2 incoming particle (one on- 
shell and the other virtual) and n + 1 outgoing real particles is given by the number 
of parameters minus the on-shell conditions and the ten parameters of the Poincare' 
group: 

4(n + 3) - (n + 2) - 10 = 3n. (14) 
These 3n variable correspond in this case to 

Q^, s, Z2, . . . , Zn, ^i, . . . , ^n, (15) 

where an azimutal angle is arbitrary. In the z ^ 1 limit, s — and the dominant 
contribution is given by terms which are most singular as s vanishes. Because of 
cancellation of infrared singularities 1^ . ~ g-n+o{e) ^j^gj^ s — 0. Indeed, 

a stronger singularity would lead to powerlike infrared divergences and a weaker 
singularity would lead to suppressed terms in the z —>■ 1 limit. 

Hence only terms in the square amplitude which behave as contribute in 

the z — > 1 limit. In d dimensions, these pick up an g-'^-^^+oie) prefactor from the 
phase space Eq. ffT^ . Let's consider the simplest case, that is the tree level case where 
we have only purely real soft emission. In this case 0(e) = and thus we get that 
the contribution to the coefficient function from the tree level diagrams with n extra 
radiated partons behaves as: 

\An+l\'d(Pn+l ~ S~'~^' f dZn ■ ■■dZ2Z^:-''^-^^-^^\l - Z^f-''' ■ ■ ■ Z^ {I ' ^2)'"'^ 

X / dVli---dVtn. (16) 



We note that each z integration can produce at most a 1/e pole from the soft region 
and that each angular integration can produce at most an additional 1/e pole from 
the coUinear region (see Eq.f l23l) in Appendix [B]) . Therefore, from the contribution of 
tree level diagrams with n extra radiated partons there come at most 

1 1 1 

(17) 



gn— 1 gTi g2n— 1 

poles in e = 0. All this means that we can write the 0{a^) to the bare coefficient 
function in d dimensions in the following form: 

CT\z, Q\ e) = m-'^^^ (1 - zr'-""', (18) 



where the factor (Q^)""*^ is due to elementary dimensional analysis, Cnn(e) are coef- 
ficients with poles in e = of order at most of 2n and the F function factor r~^{—ne) 
has been introduced for future convenience. For the LO tree level case (see Eg. (1731) 
of Chapter [1]) we have that 

C',^'\z,Q\e) = 6{l-z). (19) 
Hence, the tree level coefficient function Ctrcc ^^e z 1 limit has the form: 

oo 

cSl{z,Q\ao,e) = Y,a",C'^'\z,Q',e) (20) 



n=0 



V(0), 



, Cnn (c) , , 

Z] 



n=l 

+ 0{{l-zf), (21) 

where is the bare i.e. the unrenormalized strong coupling constant. 

We will now study how the result of Eq. (l?T|) is modified by the inclusion of loops. 
To this purpose, we notice that a generic amplitude with loops can be viewed as a 
tree-level amplitude formed with proper vertices. Contributions to the dimensionless 
coefficient function with powers of can only arise from loop integrations in the 
proper vertices. We thus consider only purely scalar loop integrals, since numerators 
of fermion or vector propagators and vertex factors cannot induce any dependence 
on s'^. Let us therefore consider an arbitrary proper diagram G in a massless scalar 
theory with E external lines, / internal lines and V vertices. It can be shown (see 
e.g. section 6.2.3 of [50] and references therein) that, denoting with P the set of 
E external momenta and Pe the set of independent invariants, the corresponding 
amplitude Ag{Pe) has the form 

Ag{Pe) = K{27r)%{P)AG{PE), 

^I-L(d-l) 

(47r)'^V2 



I 



1=1 



d(3i 



[P^(/?)]4I^+l)/2-/[Z^^(/5,P^)]/-'^V2- 



(22) 



Here, f3i are the usual Feynman parameters, PciP) is a homogeneous polynomial of 
degree L in the Pi, Dg{P, Pe) is a homogeneous polynomial of degree L -|- 1 in the Pi 
with coefficients which are linear functions of the scalar products of the set Pe, i.e. 
with dimensions of (mass)^, and K collects all overall factors, such as couplings and 
symmetry factors. 

The amplitude Ag{Pe) Eq.(l22ll depends on s only through Dg{P, Pe), which, 
in turn is linear in s. We can determine in general the dependence of Ag{Pe) by 
considering two possible cases. The first possibility is that Dg{/3, Pe) is independent 
of all invariants except s, i.e. Dg{P, Pe) = sdG^P)- In such case, Ag{Pe) depends on 
s as 



AGiPE) = - «G = - ac, (23) 



where ac is a numerical constant, obtained performing the Feynman parameters in- 
tegrals. The second possibility is that Dg{P, Pe) depends on some of the other 
invariants. In such case, Ag{Pe) is manifestly an analytic function of s at s = 0, and 
thus it can be expanded in Taylor series around s = 0, with coefficients which depend 
on the other invariants. In the former case, Eq.( !23l) implies that the s dependence 
induced by loops integration in the square amplitude is given by integer powers of 
s~^ In the latter case, the s dependence induced by loops integrations in the square 
amplitude is given by integer positive powers of s. 

Therefore, we conclude that each loop integration can carry at most a factor of 
s"" and that Eq. (l^TI) . after the inclusion of loops, becomes: 

oo 

CW(^,Q^«o,e) = Y.<^n\z,Q\e), (24) 

n=0 



C^\z,Q\e) = {Q 



2\-ne 



+ 0((l-.)°), (25) 

where again for future convenience we have defined the coefficients of (1 — z)~^~'^^ in 
terms of T~^{—ke) and where 0((1 — 2;)°) denotes terms which are not divergent as 
z ^ 1 m the e — s> limit. 
Using the identity 

and the Stirling expansion Eq.Q of Appendix |X] we get that the Mellin transform of 
Eqs.f l55|l231) in the large- limit is given by 



&'\N,Q\ao,e)=Y,<C^n\N.Q\e), (27) 

n=0 

The content of this result is that, in the large- limit, the dependence of the regular- 
ized cross section on A^ only goes through integer powers of the dimensionful variable 
(QVAT)-. 



3.2 Kinematics of inclusive DY in the soft limit 

In the Drell-Yan case the argument follows in an analogous way with minor modi- 
fication which account for the different kinematics. In this case the relevant parton 
subprocesses are: 

Q{p) + Q{p) - i\Q) + x (29) 
Q{p) + Q{p') 7*(Q) + ^ (30) 
G{p) + G{p) - 1*{Q) + X- (31) 



The recursive application of Eq.( JT2l) in Appendix [B] with m = 1, in this case gives: 

{p + p; Q,ki,..., kn) 
—-^d4>2{p + p'; kn, P„)rf0„(P„; Q,ki,...,k 

n—l) 



dMl 
27r 



d(t>2{p + p; kn,Pn) 



X / T:^d(j)2{Pn;kn-uPn-l)d(j)n~l{Pn-i;Q,ki,...,kn-2) 
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2(-fn+li kn, Pn) / „ d(f)2{Pm kn—1, Pn-l) 



27r 



X ■ ■ ■ X 







27r 



—^dMP^; h, P2)dHP2] ki. Pi), 



(32) 



where now we have defined P„+i = p + p', so M^^-^ = s and Pi = Q. The change of 
variables which separates off the dependence on {1 — z), where now z = Q'^/s, is 



M2-g2 = ^s-Q^)zn---z,. 



i = 2, . . . ,n 



(33) 
(34) 



Also here all Zi range between and 1, because Mf < M'^j^^ < Q"^. From Eq. flHHl) . we 
get: 



dMldMl_i ■ ■ ■ dMl = det ( — — ) dzndzn-i ■ ■ ■ dz2. 



where 



dzi 



dzi 



(35) 



dz„ 



dzo 



(36) 



In this case in the new variables Eg. (1331) the two-body phase space becomes 

d(f)2{Pi+i; ki,Pi) 

= N{eml,)-'+^ [(Mil - Q') - (Mf - Q')Y~''dQ, 

= N{e){Q'r'^^{s - Q'y-'^{zn ■ ■ ■ z^^^y-'^i - z,y~'^dn„ (37) 

where in the last step we have replaced (Mj^_,_^)~^+'^ by (Q^)^^"*"^ in the z ^ 1 limit as 



can be seen comparing with Eq. (!M|) . Now, substituting Eqs.( 
we finally get 



into Eq. 



d(f)n+i{p + p';Q,ki, . . . , kn) = 2n 



27T 



2\-n{l-e). 



s-Q 



2\2n-l-2ne 



X / dZnzl:^-'^-^--'>{l-ZnY-'' 

Jo 



dQn ■ ■ ■ dQi 



dZ2Z2'il- Z2y-''.i38) 



The dependence on 1 — z is now entirely contained in the phase space prefactor 



^l-2n+2ne 

m-^'^'\s-Qr-'-'"'= [g^(l-.)T'~"^ (39) 

As before, this proves that the coefficient function for real emission at tree level is 
given by 

oo 

cSL{z,Q',ao,e) = Y.alC'^'\z,Q\e) (40) 

n=0 

= 5(l-^) + 2^«o(Q ) Y(-2ne) ^ ' 

n=l ^ ' 

+ 0((l-^)°), (41) 

In this case the introduction of loops requires more care than for the deep-inelastic- 
scattering case. We shall now show the main difference between the DIS case and the 
DY far as the introduction of loops is concerned. The two-body kinematics 

(see Eq.( !T^ in Appendix [B] together with Eqs. (!H|7|) states that the radiated partons 
in the DIS case are all soft: 

*° = l<.<n-l (42) 

,0 _ ^^(^l_Mj=^(l_,„); ,, = Q^i-i. (43) 

This confirms the validity of the argument of Section [3TT] for the introduction of loops, 
because all the invariants that can appear in the function Dg{(3,Pe) in Eq. (!22|) can 
be expresses in terms of the following ones 

= -Q' (44) 
p2 = p'2 = A;2 = (45) 
p-p' ^p-ki^Q^ (46) 
ki-kjr^Q\l-z), (47) 

which are either constant or proportional to 1 — 2;, i.e. to s. Here we have used 
Eqs.( H2]|43i) of this Section, Eq.(l6]) in Chapter [2] and the fact that the definition of z 
in the DIS case Eq.(jl]) implies that 

ipr = T^-y (48) 
4z(l — z) 

In the DY case things are quite different, because in this case two-body kinematics 
(see Eq.f lT^ in Appendix [Bl together with Eqs. ( j32]l33l) gives 



ki = ^{l-z)zn---z,+^(l-Zi) + 0{{l-zY); l<i<n-l (49) 

kl = ^(1-^)(1 -;.„); s=^ (50) 



and this implies that all the invariants that can appear in the function Dg{(3, Pe) in 
Eq. (l22l) can be expresses in terms of the following ones 



V ' V = /2 

P ■ ki p' ■ ki S{1 — Z) 

ki ■ kj ~ s(l — zY- 



(51) 
(52) 
(53) 
(54) 



Hence, we see that, in general, both odd and even powers of (1 — z)""^ may arise 
adding the loops contribution to the tree level coefficient function Eq. fHT]) . Here, in 
this thesis, we will assume that odd powers of (1 — z)"*" do not arise, because it can 
be shown by explicit computations that it is the case up to order O(a^). However, 
there are possible indications that at higher orders this assumption could not be true. 
Anyway the investigations of these aspects is beyond the aim of this thesis. 

Thus, after the inclusions of loops and with our assumptions, Eq.l HTl) becomes 



oo 



n=0 

^2\-ne 



+ 0((l-.)°). 



n 



(55) 



k=i ^ ' 



-l-2ke 



(56) 



Its Mellin transform can be written in a compact way together with that of the DIS 
case Eq.([2Hl): 

oo n 
n=0 A:=0 

where a = 1 for the DIS case and a = 2 for the DY case and where the coefficients 
are those that could be obtained from the parton-level cross sections for the partonic 
subprocesses that contribute to the given process. 



ao 



+ 



(57) 



3.3 Resummation from renormalization group im- 
provement 

In this section, we want to impose the restrictions that renormalization group in- 
variance imposes on the cross section. Our only assumption is that the coefficient 
function can be multiplicatively renormalized. This means that all divergences can 
be removed from the bare coefficient function C^'^^N, Q^, ao, e) Eq. (l571) by defining a 
renormalized running coupling ^^(/i^) according to the implicit equation 



ao(/i', as(/i'), e) = /i'^a.(/i')Z("^Ha.(/i'), e) 



(58) 



and a renormalized coefficient function 

C (^N, ^, a,{^,\ e)^ = Z(^)(iV, e)C^'\N, Q\ a,, e), (59) 

where /x is the renormahzation scale (here chosen equal to the factorization one) and 
Z(°=^(as(/i^), e) and Z'^^N, as{jJ^)-, e) are computable in perturbation theory and have 
multiple poles in e = 0. The renormalized coefficient function Ci^N^Q"^ / jj? ^asi^jJi^ ,e)) 
is finite in e = and it can only depend on through Q"^ jjj?, because as(/U^) is 
dimensionless. 

The physical anomalous dimension is given by 



7(iV,a,(Q2),e) 



-6(«oQ-^^) g^^^^_2,) lnCW (iV,g2,ao,e) , (60) 



where we have exploited the fact that Eq. (!Fn) depends on Q"^ through the com- 
bination aoQ"^*^. This implies that the physical anomalous dimension 7 has the 
following perturbative expression: 



7(iV,a,(Q2),e) = 5^7,,.(e)[ao(Q')-f ^"^'^ 

i=0 j=0 



(61) 



The renormalized expression of the physical anomalous dimension is found ex- 
pressing in this equation the bare coupling constant in terms of the renormalized one 
by means of Eq.( l58i) . Now, the functions 

(Q')-^«o = as{^?)Z'^^^\as{^^^),e) (62) 



/i2 



. e) (63) 



are manifestly renormahzation group invariant, i.e. //^-independent. Thus, it follows 
that 

{Q^Y^ao = a.(Q2)Z("=)(a,(Q2),e) (64) 

{Q^N'^r'ao = as{QyN'')Z^"^\as{QVN''),e). (65) 

The renormalized physical anomalous dimension is then found by substituting Eqs. 
into Eq. (l6ip and re-expanding Z^"'"'^ in powers of the renormalized coupling. We ob- 
tain: 



00 rn / 1 \ 

7(iV,a,(Q2),e) = 7l(e)«r"(Q')<(QViV") + ( - . 

m=l n=0 ^ ^ 



(66) 



At this point, we cannot yet conclude that the four- dimensional physical anoma- 
lous dimension admits an expression of the form of Eq.( |66l) . because the coefficients 



7mn(^) necessarily finite as e ^ 0. In order to understand this, it is convenient 

to separate off the A^-independent terms in the renormahzed physical anomalous di- 
mension, i.e. the terms with n = in the internal sum in Eq. flUUl) . Namely, we 
write 

7(iV, a,(Q2), e) = 7^')(«.(Q'), o^siQ'/N''). ^) + l^''\o^s{Q'), e) + O (^^^ , (67) 
where we have defined 

oo oo 

^^\as{Q\as{Q'/N^),e) = Y^H^^^nS^^^^iQ^XiQ' /N^) (68) 

m=0 n=l 

oo 

l^^\as{Q%e) = 5^7l(6)ar(Q')- (69) 



m=l 



Whereas 7(A^, as(Q^), e) is finite in the limit e — > 0, where it coincides with the four- 
dimensional physical anomalous dimension, 7*^'-' and 7'-'^'' are not necessarily finite as 
e ^ 0. However, Eq. (l67|) implies that 7*^'^ and 7'^'^-' can be made finite by adding and 
subtracting the counterterm 

zW(a.(Q2),6) =7(')(a,(Q^),a,(Q2),6). (70) 

In this way the physical anomalous dimension Eg. ( 1671) becomes 

7(iV, a,(Q2), e) = ^'''\as{Q^), a^iQ^N^), e) + 7^^H«s(Q'), e) + , (71) 



where 



l^'\a,{Q^),a,{QyN'^),e) = 7^')(a,(g2), (..(QViV'^), e) + 

-7(')(a,(Q^),a.(Q2),e), (72) 
j^^\as{Q'),e) = ^^^\as{Q'),e)+^(^\a,{Q'),a,{Q'),e). (73) 

Now, 7*^^^ is clearly finite in e = 0, because at = 1 7*-'^ vanishes and it is A^- 
independent. This also implies that 7'-'^ is finite for all A^, because 7 shoul be finite 
for all A^. Therefore, 7^^^ provides an expression of the resummed physical anomalous 
dimension in the large A^ limit, up to non-logarithmic terms: 

7(Ar, a.m, e) = 7^'H«.(Q'), c^s{QVN^), e) + 0{N'). (74) 

It is apparent from Eq. (!73|) that 7'-'^'' is a power series in as(Q^) with finite coefficients 
in the e — limit. In order to understand the perturbative structure of 7^'-* as well, 
define implicitly the function g{as{Q^) , asiQ"^ /n) , e) as 

i'\as{Q'),a,{Q^/N'^),e)= / — ^^(^.(Q^), a.(QVn), e), (75) 



where 

^7(a,(Q2),a,(/i2),e) = -^^l-A^W^asm,as{fi'),e) (76) 

= -/3W(a,(/x2),e)^-^7^^^ (77) 

with p^'^^as) is the (i-dimensional beta function 

/?W(«,(/i2), e) - ea,(/i2) + /^(^.(/i^)) (78) 

and where we have performed the change of variable n = Q'^/fi'^. It immediately 
follows from Eqs. fl68ti77|) that g is a power series in as(Q^) and asi^"^) with finite 
coefficients in the limit e ^ 0: 



limfi'(as((5^),«s(/W^),e) = fi'(as(Q^), «s(/i^)) 

e—>IJ 

oo oo 



m=0 n=l 



Hence, our final result for the four-dimensional all-rder resummed physical anomalous 
dimension is given by 



dn 

Y'^^N, asm) = I — asiQ'/n)) + 0{N'). (80) 



n 



This result can be compared to the all-order resummation formula derived in 
Ref ■ [5T\ . This resummation has the form of Eq. (l80|) . but with g a function of as(/i^) 
only, i.e. with all gmn = when m > and our result is thus less predictive in 
the sense that it requires a higher fixed-order computation of the physical anomalous 
dimension in order to extract the resummation coefficients gmn- The predictive power 
of the resummation formulae is analyzed in detail in Chapter [71 According to the 
more restrictive result of Ref. [51] , Eq. (IHOil becomes 



dn 

r'{N,as{Q')) = / -g{as{QVn)) + 0(iV°), (81) 

where 

oo 

g{as{fi')) = Y,9ona:{fi'). (82) 

n=l 

The conditions under which the more restrictive result of Ref. [SI] holds can be 
understood by comparing to our approach the derivation of that result. The approach 
of Ref. [51] is based on assuming the validity of the factorization formula Eq. (H71) of 
Sec J2.3.2l which is more restrictive than the standard coUinear factorization. This 
factorization was proven for a wide class of processes in Ref. [11] , and implies that the 
perturbative coefficient function Eq. flFT]) in the large N limit can be factored as: 



C(°)(iV, Q\ ao, e) = C(°'')(QViV^ ao, e)C^'''^\Q', ao, e). (83) 



We notice that this can happen if and only if the coefficients (e) in Eq. fl^ can 
be written in the form 

C^i°2(e) = ^.(e)G'n_.(6). (84) 

The vahdity of factorization Eq. lHTl) of Sec J2.3.2l to all orders and for various processes 
is based on assumptions whose reliability will not be discussed here. Anyway, Eq. flHHl) 
implies that the physical anomalous dimension Eq. flUUj) becomes 

7(iV, «.(Q2), e) = 7«(a.(gViV'^), e) + 7^^H«.(g2), e). (85) 

Then, proceeding as before, one then ends up with the resummation formula Eq. fl8ip . 

3.4 NLL resummation 

In this Section, we shall give explicit expressions of the reummation formulae at NLL 
for the deep-inelastic structure function F2 and for the Drell-Yan cross section. These 
explicit expressions are useful for practical computations and we shall use them in 
Chapter O 

The expression of the resummed physical anomalous dimension in Eq. fl80p can be 
used to compute the resummed evolution factor K^^^i^N^Ql.Q'^) Eq.(fT9l) in Section 
12:21 At NLL, we get 

i^NLL(A^; Ql Q') = exp [Enll(A^; Ql Q')] = exp / ^ -^7n£l(A^, «.(^')), (86) 
where 

Tjr^jAT, as{k^)) = / — [goia,{e/n) + go2al{e/n)] ; g^, = 0. (87) 
Ji ^ 

The fact that the resummation coefficient gn = for both the deep-inelastic and the 
Drell-Yan case can be shown by explicit computations of the fixed-order anomalous 
dimension (see Ref. [8] in Section 4.3). This shows that at NLL level Eq.i lHTj) holds. 
However this does not mean that this is the case at all logarithmic orders. 

Many times in literature, the resummed results are given in terms of the Mellin 
transform of a resummed physical anomalous dimension in z space. To rewrite Eq. flHTl) 
as the Mellin transform of a function of x, we can use the all-orders relations between 
the logs of and the logs of 1 — z that are given in Appendix [XI In particular, we 
can use Eqs. ( 123|2^ in Appendix [A] to rewrite Eq. ( IHTI) at NLL in the following form: 

rl N-l _ 1 

YnldN, = a dx — [goias{k\l - xT) + go2a',{k\l - xT)] , (88) 

Jo i — X 



where 



goi — —goi, 



§02 = -{go2 + ajEbogoi). 



(89) 



and where we have used the definition of the beta function Eg. (1431) in Section [1.21 
As a consequence, the NLL resummed exponent in Eq. (!86l) can be rewritten in the 
following form 



N-l 



1 /'^'(l-^)'' rfP 



where 



(90) 



(91) 



Beyond leading order the standard anomalous dimension differs from the physical 
one, so gQ2 receives a contribution both from the standard anomalous dimension and 
from the coefficient function. It is thus natural to rewrite the resummation formula 
Eq. flUUl) separating off the contribution the contribution which originates from the 
anomalous dimension Eq. (1181) of Section 12.21 This is done defining two functions 
of as, A{as) and B"-{as) in such a way that 

g{as)=A{as) + ^^^^^^^^, A{as) = A,as + A2al B^^\a s) = Bf"^ a s. [92) 

It is clear that the constant Ai are obtained directly form the coefficients of the 
1/[1 — terms of the i-loop quark-quark splitting functions and that the coefficients 
Bf depends on the particular process (a = 1 for DIS and a = 2 for DY). To the NLL 
order, we find that 

901 = A,, go2 = A2 - boBi"^ (93) 

and that 



a dx 



X 



N-l 



1 — X 



Aiasik')) 



+5(^)(a,(g2(l - xT)) - B^^\asiQl{l - x)'^)) 

We can then rewrite the resummed cross section 

ctnll{N, Q^) = exp [ENLL(iV; Ql g')] aNLL(iV, Qo) 



(94) 



(95) 



in a factorized form according to Eqs. ( |TT][T2i) in Section 12.21 by collecting all Q^- 
dependent contributions to the resummation Eq. (l94|) into a resummed perturbative 
coefficient function Cnll: 



where 

CNLL(iV,QV/i',«.(/^')) 



aNLL(iV,Q') = CNLL(iV,QV/^',as(/x'))F(iV,/i') 



(96) 



exp {a I dx 



1 — X 



^A{as{k^)) 



+B'^^\as{Q\l~xr)) 



(97) 



which ha the same form of the resummed results discussed in Section I2.3.1[ The 
precise definition of the parton distribution F and the factorization scale /i^ will 
depend on the choice of factorization scheme: according to the choice of scheme, the 
resummed terms will be either part of the hard coefficient function Cnll? oi of the 
evolution of the parton distribution F. In the MS scheme the NLL coefficients Ai, 
A2 and B^f^ are given in Eq. ( H6i) of Section 12.3.11 and the NNLL ones are given for 
example in Ref.[52]. These coefficients can also be obtained comparing a fixed order 
computation of the physical anomalous dimension with a fixed order expansion of 
Eq-dHOD as is shown explicitly in Section 4.3 of Ref.[H]. 

To compute explicitly Eq. (!U7|) . we first exploit Eq. (!2Ul) in Appendix 1X1 at NLL 
level, thus finding 



Cnll(A^, QV/^', as(yu')) = exp 



AT" 



dn 
n 



dk 



— A{as{k^/n)) + B^'^\as{Q'/n)) 



(98) 

where B^"-\as) = B^"-\as) — ajEAiag- The exphcit expression of Eq. (l9711 is then 

dasiQ'^/n) 



obtained performing the changes of variables, 
dk"^ das{k'^/n) dn 



(3{as{k^/n)y n 



/5(a,(QVn))' 



(99) 



to evaluate the integrals in Eq. ( l98l) and using the two loop solution of the renormalitazion- 
group equation for the running of given in Eq. fl4ip of Section II. 2[ Now, after some 
algebra we find for the integral in Eq. ( 



AT" 



dn 
n 



dk^ 



^2 

^ + A2aj^ 



n 



n 



Aia 

log Ngi{X, a) +5(2, (A, a) 



B\ 'as 



^ n 



where A 



bQas{fil)\ogN and 
A, 



[aX + (1 - aA)log(l - aX)] 



(a) 



Aia-jE - B\ 

bo 
A2 

■--2-[aA + log(l - aX)] + log 
"0 



Aih 



■log(l-aA) + -^[aA + log(l 



aX) + -\og\l 



(100) 

(101) 
■aX)] 



1-log 1 

bo 



aX) 



+ log 



/i^A Ai 



111 



aA, 



(102) 



where a = 1 for the DIS case and a = 2 for the DY case. Evidently, in Eqs. 
there is a dependence on the renormalization scale. To obtain the desired result, we 
simply have to keep the renormalization scale equal to the factorization scale /i 
Thus, for the explicit analytic expression of Eq.( l98i) . we get 



/i 



Cnll ( N, ^,asin^] 



[exp{logiV5(i(A,a) +5(2 (A, a)}] 



(103) 



with the resummation coefficients given in Eq. (H6l) in Section [2.3.1l for the MS factor- 
ization scheme choice. However, a general analysis of the factorization scheme choices 
and changes for the resummation formulae is given for example in Section 6 of Ref. [8j. 



Chapter 4 



Renormalization group 
resummation of prompt photon 
production 

In this chapter, we prove the all-order exponentiation of soft logarithmic corrections to 
prompt photon production in hadronic collisions, by generalizing the renormalization 
group approach of chapter [31 Here, we will show that all large logs in the soft 
limit can be expressed in terms of two dimensionful variables. Then, we use the 
renormalization group to resum them. The resummation formulae that we obtain 
are more general though less predictive than those that can be obtained with other 
approaches discussed in chapter [21 



4.1 Kinematics and notation 

We consider the process 

H,{P,) + H2{P2)^l{p^)+X, (1) 

of two colliding hadrons Hi and H2 with momentum Pi and P2 respectively into a real 
photon with momentum and any collection of hadrons X. More specifically, we are 
interested in the differential cross section p\-^{x±,p\), where p± is the transverse 
momentum of the photon with respect to the direction of the colliding hadrons Hi 
and H2, and 

XX = S = (Pi + P2f. (2) 

The scaling variable x can be viewed as the squared fraction of transverse energy 
that the hadrons transfer to the outgoing particles (hence < xj, < 1) and 5* is 
the hadronic center-of-mass energy. We parametrize the momentum of the photon in 
terms of its partonic center-of-mass pseudorapidity and its transverse momentum p^. 
The pseudorapidiry of a massless particle is defined in terms its scattering angle 9 in 
the center-of-mass frame as follows 

7] = -ln(tan(^/2)). (3) 
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So, in the partonic center-of-mass frame, we can write: 

Pf = ip± cosh fly, p± , p± sinh fj^). (4) 
In the same frame, the incoming partons' momenta can be written as 

Pi = xiP, = ^(l,Ox, 1), P2 = X2P2 = ^(1,0±, -1), (5) 

where Xi(2) are the longitudinal fraction of momentum of the parton 1(2) in the hadron 
ifi(2) and s = {pi +P2Y = X1X2S is the center-of-mass energy of the partonic process. 
The relation between the hadronic center-of-mass pseudorapidity and the partonic 
one is obtained performing a boost along the collision axis: 

f)y = 7]y-^ln^. (6) 
^ X2 

The factorized expression for this cross section in perturbative QCD is 

d(j 

p\- — {x±,p]_) = / dxidx2dzxiF^\xi,iJ?)x2F^^{x2,n'^) 

0,0 



Q 



2 



xCab — , asi/jL )j S {X± - ZX1X2) , (7) 

where F^^[xi, ^'^), F^^[x2.,n'^) the distribution functions of partons a,b in the 
colliding hadrons. Here we have defined the perturbative scale and the partonic 
scaling variable z as the squared fraction of transverse energy that the parton a, b 
transfer to the outgoing partons {0 < z < 1): 

= 4pi, (8) 

, = = (9) 

S X1X2S 

The coefficient function Cab{z, ^,q;s(//^)) is defined in terms of the partonic cross 
section for the process where partons a, b are incoming as 

C„,(^.,^,«>^))=pi^. (10) 

4.2 Leading order calculation 

The hard-scattering subprocesses that contribute to the amplitude of the prompt- 
photon prodution at the leading order are: 

q{pi) + q{p2) 9{p')+l{P^) (11) 

q{pi) + 9{p2) q{p') + i{p^) (12) 
q{pi) + 9{p2) q{p') + i{p^) (13) 



iM 






Figure 4.1: Feynman graphs that contribute to the amphtude of the partonic process 
at the leading order 



The corresponding Feynman graphs are shown in figure 14.11 

We now want to obtain the coefficient function for these two elementary sub- 
processes at the leading order. Using the QCD Feynman rules to evaluate the first 
amplitude in figure 14.11 we have: 



iM, 



-iv{p2)Qqeg 



2pi ■ 



2pi ■ p' 



<(P7)e;(p')iXpi), 



(14) 

where e = |e| is the electrical charge, g is the strong charge and Qq is the electrical 
charge of the quark q in units of e. Taking the square modulus of equation ffT^ and 
averaging over the two quarks polarizations and colours we obtain: 



--y 

A All. ^ 



\M - ? 



C 



poljCol 



^ Pi 



p' 



Pl-P' Pl-P-y 



(15) 



where Nc is the number of quark colours and Cp = {N^ — l)/2Nc is the Casimir 
operator with respect to the colour matrix 

Proceeding in the same way for the Feynman graphs contributing to the second 
amplitude in figure 14.11 we could obtain the averaged square modulus for this other 
subprocess, but, we note that it can be immediately obtained using crossing simmetry. 
More precisely, we only have to substitute p' with p2 and take into account the fact 
that we must now average not over two quarks' colours but over the colours of a quark 



and a gluon. Therefore we arrive at the following expression for the averaged square 
modulus of the second amplitude of figure 14.11 



4 2NcCf 



1 



poljCol 



4 2NcCf 



poljCol 



P'-*P2 



C 



Pl-Pj ^ Pi ■P2 

lPi-P2 Pi-P-ri 



(16) 



We will now rewrite equations fllSI) and fll6l) in terms of the kinematic parameters 
defined in section 14. 1[ For this purpose we first note that momentum conservation 
and the parametrizations defined in Eqs. (H]l5ll imply (in the centre of mass of the 
incident partons): 



p' = {p± cosh r]j, —p± , —p± sinh 17^) 

s = iPi + Pif = (Py + p'Y = ^p]_ cosh^ f]^ 
Apl 1 



(17) 
(18) 

(19) 



s cosh 17^ 

From the last equation we see that in the limit 2; — > 1, 17^ ^ 0. Physically this is 
because in this limit all the centre of mass energy is transverse and so the photon 
cannot have a non zero pseudorapidity. Now, using again Eqs .( 141151) and the equation 
for p' (1T7|) we have: 



Pi ■P2 



Pi -P 



Pi -Pi 



s 
2 

2 



P±e 



(20) 



-p±e 



A combination of equations (IT^ and 
square modulus of the amplitudes: 



^ poljCol 



yields our final result for the averaged 

(21) 



4^nVq2 (2-^) 



1 1 



4 2NcCf Nc 



1(1)9^11(9) I 



poljCol 



2iVc^' ^ 



l± VT 



1 + ' ^22) 

1 ± Vl - 

where the plus sign has to be chosen for positive values of the pseudorapidity f/^ and 
the minus sign for negative values. The two-body phase space is: 



d(f){pi+p2;pj,p) 



d^p^ 



d'^p' 



(27r)32Ep^ {2nf2Ep, 
d^p^ 1 



(27r)^5(4)(p' -pi -P2) 



4EpEp> (27r)2 



6(i){Ep' + Ep^ - Ep^ - Ep^). 



(23) 



Imposing the conservation of spatial momentum {Ep^^ = Epi = the two-body 
phase space becomes: 



d(f){pi+p2;p^,p) 



d?p.^ 



4|j9;|2(27r)2 

1 



^{1)^1 



167r 



d cos9d\p^\8{i){\p^ 



sl2) 



Y^'^(i)(P-L cosh 17^ — \fl/2)d COS 9 d\p. 



71' 



(24) 



where 9 is the scattering angle of the photon with respect to the collision axis. Because 
of the fact that we want to integrate over the pseudorapidity of the photon 77^ at fixed 
p_L, we must perform a change of variables. In particular we must rewrite the two- 
body phace space, expressed in terms of cos^ and in terms of the new variables 
cosh 57^ and j9_l. This change of variables is given by the equations: 

cos 9 = tanh fj^ 



p± cosh 7}^. 



The determinant of the Jacobian matrix is easily obtained and is given by: 



\J\ 



d{cos9, \pl^\) 



9(cosh fiy,p±] 



1 



(25) 



cosh 17^ A / cosh 17^ — 1 



Thanks to equation (fT9l) . we obtain this determinant in term of z: 

z 



Using this last result, equation (12^ becomes: 

1 



(26) 



d(f){pi+p2;p^,p') 



— 7=^=5(1) (cosh J7^ - ^/s/2p±)^^dcoshTi^ 

lOTT y/1 - Z P± 



- v+) + ^wiv^ - V-)] 



dp±,. 
■df]^, 



X- 



P± 



(27) 



where 17+ and r}- are the two solutions (one positive and one negative) of the equation 
imposed by the delta function p± cosh 57^ = which are: 



ri± = In 



± 



1 



2p± V ^pI 

The flux factor $ is immediately obtained from equations ([5]) and (fT9|) : 



(2^ 



$ = A{pi ■ P2) 



2s 



ipl 



(29) 



Remembering the definition of tfie QED and QCD coupling constants, 



and putting togetlier expressions (^2), (122]), (I2ZI), and performing the integration 
over fj^, we obtain our final result for the coefficient function at the leading order for 
the two subprocesses in figure 14. H 



r^^^^ (z cy ') 



[l-z] 



" ' 2Nc 



1 + 



(30) 
(31) 



4.3 The soft limit 

We will study the cross section Eq. ([7]) in the threshold transverse limit, when the 
transverse transverse energy of outgoing particles is close to its maximal value {x± — > 
1 at the hadronic level or, equivalently, z — »• 1 at the partonic level). The convolution 
in Eq. is turned into an ordinary product by Mellin transformation: 

a(iV,Q2) = ^a,,(iV,g2) (32) 

a,b 



a,b 

xCab\^N,^,a,{fi')j. (33) 

As discussed in Section 12.21 in the large limit each parton subprocess can 
be treated independently, specifically, each Cab is separately renormalization-group 
invariant. 

At this point, it is interesting to discuss the differences in the large A^ behavior 
of the partonic subprocesses. The cross sections for the partonic channels with two 
quarks of different flavors {ab = qq', qq', qq, qq', qq, qq') vanish at LO and are hence 
suppressed by a factor of with respect to the subprocesses with ab = qq, qg, qg. 
Moreover, in the large A^ limit this relative suppression is further enhanced by a 
factor of 0(1/A^) because the processes with two different quark flavors involve the 
off-diagonal Alatarelli-Parisi splitting functions. Therefore, these partonic channels 
do not contribute in the large A^ limit. The partonic channel ab = gg has a different 
large A^ behavior. It begins to contribute at NLO via the partonic process g + g ^ 
J + q + q, which again leads to a suppression effect of 0{1/N) with respect to the 
LO suprocesses. However, owing to the photon-gluon coupling through a fermion 
box, the partonic subprocess g + g ^ + g \s also permitted. This subprocess 
is logarithmically-enhanced by multiple soft-gluon radiation in the final state, but it 
starts to contribute only at NNLO in perturbation theory. It follows that the partonic 



channel ab = gg is suppressed by a factor with respect to the LO partonic channels 
ab = qq, qg, qg and it will enter the resummed cross secion only at NNLL order. In 
conclusion, the partonic channels that should be resummed are ab = qq, qg, qg, gg, 
where the last channel is that coupled to the gluon via a fermion box and eneters 
resummation only at NNLL. 

Furthermore, on top of Eqs. (CI [33]) the physical process Eq. ([T]) receives another 
factorized contribution, in which the final-state photon is produced by fragmentation 
of a primary parton produced in the partonic sub-process. However, the cross section 
for this process is also suppressed by a factor of in the large N limit. This is 
due to the fact that the fragmentation function carries this suppression, for the same 
reason why the anomalous dimensions jqg and jgq are suppressed. Therefore, we will 
disregard the fragmentation contribution. 

According to Eqs.( fT9ll22l) of Section [?!2l the cross section can be written in terms 
of the physical anomalous dimensions: 

a{N,Q^) = Y,Kab{N;Ql,Q^)a,k{N,Ql) (34) 

b 

J^exp [Eat{N;Ql,Q')] aab{N,Ql), (35) 



a,b 



a,b 



where 



E,,{N;Ql,Q') = I ^^,,{N,a,{e)) (36) 

+ In CabiN, 1, a,(Q2)) _ in C,,(iV, 1, ^.(Q^)). (37) 
In the large-xx limit, the order-n coefficient of the perturbative expansion of 



1 — X I 



with 



the hadronic cross section is dominated by terms proportional to 

k < 2n — l, that must be resummed to all orders. Upon Mellin transformation, these 
lead to contributions proportional to powers of In -i:. In the sequel, we will consider 
the resummation of these contributions to all logarithmic orders, and disregard all 
contributions to the cross section which are suppressed by powers of (1 — x±), i.e., 
upon Mellin transformation, by powers of jj. 

The resummation is performed in two steps as in chapter [31 First, we show that the 
origin of the large logs is essentially kinematical: we identify the configurations which 
contribute in the soft limit, we show by explicit computation that large Sudakov logs 
are produced by the phase-space for real emission with the required kinematics as logs 
of two dimensionful variables, and we show that this conclusion is unaffected by virtual 
corrections. Second, we resum the logs of these variables using the renormalization 
group. 

The l-th order correction to the leading 0{as) partonic process receives contri- 
bution from the emission of up to / + 1 massless partons with momenta ki, . . . , 



Four-momentum conservation implies: 



Pi + 'P2= Vi + ki + . . . h+i. (38) 
In the partonic center-of-mass frame, according to Eqs. (Hf5l) . we have 

(Pi +P2 = — (1 - v^coshr/^) = k^k'^il - cos%) > 0, (39) 

where 9ij is the angle between ki and kj. Hence, 

l<coshr}^<^. (40) 



Therefore, 



J2 kX,a - cos%) = ^il-z) + [(1 - zY] . (41) 

Equation fl4ip implies that in the soft limit the sum of scalar products of momenta 
ki of emitted partons Eq. (|39l) must vanish. However, contrary to the case of deep- 
inelastic scattering or Drell-Yan, not all momenta ki of the emitted partons can be soft 
as z ^ 1, because the three-momentum of the photon must be balanced. Assume thus 
that momenta fcj, z = 1, . . . , n; n < / + 1 are soft in the 2; — > 1 limit, while momenta 
ki,i> n are non-soft. For the sake of simplicity, we relabel non-soft momenta as 

k'j = kn+j] i < j < m + 1; m = I — n. (42) 

The generic kinematic configuration in the z = 1 limit is then 

fcj = l<'i<n 

m+l 

% = 0; ^fc;°=px l<z,j<m+l. (43) 

j=i 

for all n between 1 and Z, namely, the configuration where at least one momentum is 
not soft, and the remaining momenta are either coUinear to it, or soft. 

With this labelling of the momenta, the phase space can be written, using twice 
the phase space decomposition of Eq. (fT2i) in Appendix [HI as 

P2] P7) ^1; • • • ; ^n; ^1; • • • ; ^m+l) 

-,2 



J '^^^n+i{pi + P2;q,ki,---,kn) 

■ d(f)m+i{k'; k[, . . . , fc^+J d(j)2{q; p^, k'). (44) 



We shall now compute the phase space in the z 1 limit in c? = 4 — 2e dimensions. 
Consider first the two-body phase space d(f)2 in Eq. (jSj). In the rest frame of q we 



have 



87rr(l - e) 



sin " 6^^ d Iplyl (icos6'^ — P) 



where 



P 



(45) 



(46) 



Because momenta i < n are soft, up to terms suppressed by powers of 1 — z, the 
rest frame of q is the same as the center-of-mass frame of the incoming partons, in 
which 



1^7 1 = P± cosh?}^ 
COS 9^ = tanhr}^. 



Hence, 



d(f)2iq;Py,k') 



{Any (QV4)- 
87rr(l - e) ^ 



dp± dfj^ S cosh fj^ — 



2P 

7^ 



The conditions 



coshr)^ = — := > 1; /c'^ > 0, 



together with Eg. ( 14611 . restrict the integration range to 



0< k'^ <q^-VQ 



2n2 



It is now convenient to define new variables m, v 



q^ = Q^ + u{s-Q^ 



Q^[l+u{l- z)\+0{{l- z) 



k'' = v{q^ - ^/W) = Q'\uv{l-z) + 0{{l-zf) 



< M < 1 



<v <h 



in terms of which 



P 



+o[{i-zr] 



Thus, the two-body phase space Eq. (I49l) up to subleading terms is given by 
d(j)2{q;p^,k') 



87rr(l -e) 



^u{l-v){l-z) 



(47) 
(48) 



(49) 



(50) 



(51) 
(52) 



(53) 
(54) 
(55) 

(56) 



(57) 



where 



In 



2P 

7q' 



± 



/4P2 



(5^ 



We now note that the phase-space element d(f)n+i{pi + P2', Q, ki, . . . , kn) contains 
in the final state a system with large invariant mass > Q^, plus a collection of 
n soft partons: this same configuration is encountered in the case of Drell-Yan pair 
production in the limit zdy = (f I s — >■ 1, discussed in Section [X^ . Likewise, the phase 
space for the set of coUinear partons d(j)m+i{k'; k[, . . . , k'^_^_{) is the same as the phase 
space for deep-inelastic scattering (discussed in Section [3?Ti) . where the invariant mass 
of the initial state fc'^ vanishes as 1 — z (see Eq. (15^ ). We may therefore use the 
results obtained in chapter [3l where, in the case of deep-inelastic scattering, one of 
the outgoing parton momenta (fc^+i' ^^y) identified with the momentum of the 
leading-order outgoing quark p'. Hence Eq. flUUj) is obtained from the corresponding 
result in chapter [3] for deep- inelastic scattering by the replacement p' — > fc^+i: 



{Pl +P2]q, ki,...,kn) 

d(j)^+i{k';k[,...,k'^^^) = 2tt 
where iV(e) = l/(2(47r)2-2^) and 
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dZfi 





{k 



/2\m— 1- 



(59) 
(60) 



-2e 



dn'^'^H 



X f dz2z\-^\l - Z^f-''' 

Jo 

dn\ . . . dn'^ f dz'Ji--'^-^--'^^{i - z'j^' 

Jo 



(61) 



2e 



X / d4z;-(i - 4)l-2^ 



(62) 



Here, the variables defined as in chapterOfor the Drell-Yan and deep- inelastic 

scattering respectively. 

Equations fl53f54l) imply that the phase space depends on (1 — 2;)^^ through the 
two variables 



/2 ^2, 

{s - r 



k'' o^Q^{l-z) 



ocQ'{l-z) 



(63) 
(64) 



where the coefficients of proportionality are dimensionless and z- independent. By 
explicitly combining the two-body phase space Eq. (JSTj) and the phase spaces for soft 



radiation Eq. ( ]59l) and for collinear radiation Eq. (I60i) we get 

{Pl + 

Ari {^ ^\2n+m-(2n+m)e 

_ ( f^2\n+m-(n+m+l)t "/^-L \^ ^) 
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u Jo \/\-v 



(iw ^ / (if — [(^(f/^ — r),),) + biji^i — fj^] 



(65) 



In the limiting cases n = and m = we have 



#i(pi+P2;g) = 27r<5(s-g2) = _i!L_5(i_^) (ee) 

o 47r 

rf0i(fc';p') = 27i6{k") = --— -6{vy, (67) 

- z) 

the corresponding expressions for the phase space are therefore obtained by simply 
replacing 

(1 -M)-ic/fi(")(e) (5(1 -m); v"Un'^"'\e) 6{v) (68) 

in Eq. (lU^ for n = 0, m = respectively. 

The logarithmic dependence of the four-dimensional cross section on 1 — 2; is due 
to interference between powers of (1 — 2;)^^ and ^ poles in the c?-dimensional cross 
section. Hence, we must classify the dependence of the cross section on powers of 
(1 — 2;)"'^. We have established that in the phase space each real emission produces 
a factor of [Q^(l — z)"^]^ if the emission is soft and a factor of [(5^(1 — z)]^ if the 
emission is collinear. The squared amplitude can only depend on (1 — 2)"^ because 
of loop integrations. This dependence for a generic proper Feynman diagram G will 
in general appear, as discussed in chapter [3], through the coefficient (see Eq. (l22i) of 
Section 13.11) 

[DaiPE]]'''^'-', (69) 

where L and / are respectively the number of loops and internal lines in G, and 
Dg{Pe) is a linear combination of all scalar products Pe of external momenta. In 
the soft limit all scalar products which vanish as z ^ 1 are either proportional to 
Q2(l - z) or to g2(l - z)2 as shown in Eqs. fH71l53]l5iD of Section O • Equation (El) 
then implies that each loop integration can carry at most a factor of [Q^(l — z)"^] ^ 
or [Q^{l-z)r- 

This then proves that the perturbative expansion of the bare coefficient function, 
for each sub-process which involves partons a, b, takes the form 



CW(^,Q^ao,e) = ««o(Q')"^ ^(z, Q^ e) (70) 

(71) 



1=0 



where the factor l/r(l/2) was introduced for later convenience and terms C^^^, with 
k + k' < I a.t order a[ are present in general because of loops. The coefficients Cl^l, 
have poles in e = up to order 21. To understand this, we have to count the indepen- 
dent variables for the prompt photon process. We have 2 incoming particles and / + 2 
outgoing partons (a leading-order parton, the photon and / extra emissions). There- 
fore, imposing the on-shell conditions and the constraints due to Poincar invariance, 
we get 

4(/ + 4) - (/ + 4) - 10 = 3/ + 2, (72) 

independent variables. Now, we need to understand which are these independent 
variables: the general expression of the phase space in the threshold limit Eq.( |65ll is 
written in terms of 3/ -|- 3 variables which are 

where n+m = I and each solid angle depends on two parameters. Clearly, one of them 
must be a function of some of the others because of Eq. [721 In fact, from Eq. fl^Sj) . 
we know that rj-y depends on u, v, 4p^, s. Thus, the 3/ -|- 2 independent variables on 
which depends the square modulus amplitude can be chosen as 

Now, each of the / integrations over a solid angle can produce a pole 1/e from the 
collinear region. Furthermore, each of the I integrations over a dimensionless variable 
u,v , Z2, ■ ■ ■ , Zn, Z2, . . . , z'^ can produce a pole 1/e from the soft region. This explains 
why the coefficients C;^°^, can have poles in e = up to order 21. 



4.4 Resummation from renormalization group im- 
provement 

The Mellin transform of Eq. flTUj) can be performed using 

r(i/2) 



Jo Vn \nJ 



with the result 



C(°)(iV,Q^ao,' 
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x[(gr^ao] + 0(1). 



(75) 



(76) 



Equation ( 1761) shows that indeed as — ^ cxd, up to corrections, the coefficient 

function depends on through the two dimensionful variables ^ and The argu- 
ment henceforth follows in the same way as in chapter 131 in this more general situation. 



The argument is based on the observation that, because of coUinear factorization, the 
physical anomalous dimension 

7(iV,a.(Q2)) = Q'-^^\nC{N,Q^/^^\as{^^')) (77) 

is renormalizat ion-group invariant and finite when expressed in terms of the renor- 
malized coupling Q;s(/i^), related to by 

«o(/i',«s(/i')) =/i'^«.(/i')^^"^n«s(/i'),e), (78) 

where Z("»)(a,(/i2), e) is a power series in Q;s(/i^). Because is manifestly indepen- 
dent of /i^, Eq. ( ITHl) implies that the dimensionless combination ((5^)~''ao(«s(/^^)5 /^^) 
can depend on Q"^ only through asiQ"^): 

(Q2)-^ao(as(/x'), /i') = asm ^(°=^(as(Q'), e). (79) 

Using Eq. ( 1791) in Eq. ( 1761) . the coefficient function and consequently the physical 
anomalous dimension are seen to be given by a power series in ^^(Q^), as{Q'^/N) and 

oo oo oo 

7(iV,a.(Q^),e) = J] J] J] 7Lp(e) «r(Q')<(QViV') «?(QViV)- (80) 

m=0 n=0 p=0 

Even though the anomalous dimension is finite as e ^ for all N, the individual 
terms in the expansion Eq. (1371) are not separately finite. However, if we separate 
A/^-dependent and A^- independent terms in Eq. ( 1571) : 

7(iV, asm, e) = 7^^H«.(Q2), e) + 7^')(iV, as(Q'), e), (81) 

we note that the two functions 

j'^^\asm,e) ^ 7^^)(«.(Q'),e)+7^')(l,a.(Q2),e) (82) 
7«(A^,«,(Q2),e) ^ ^^^\N,asm,e)-l^^\l,asm,e) (83) 

must be separately finite. In fact, 

7(iV, asm, e) = l^'\asm, e) + 1^'\N, asm, e), (84) 

is finite for all N and 7*^'-* vanishes for = 1. This implies that 'j^^\as{Q'^), e) is 
finite in e = and that j^'-^N, asiQ"^), e) is also finite in e = for all N because of 
the A^- independence of 'y^^\as{Q'^), e). 
We can rewrite conveniently 

7«(iV,«.(Q2),6)= / —g{asm,o^s{Q^/n''),as{Q^/n),e), (85) 

where 

g{asm, as{QVn^), asiQ^n), e) = n-^^^'\n, asm, e)- (86) 



is a Taylor series in its arguments whose coefficients remain finite as e — 0. In four 
dimension we liave tlius 



7(iV,«.(Q^)) = j^'\N,asm,0)+i^''\N,asm,0) + O 



j^'\N,a,m,0) + O{N') 
^ dn 

— asiQ'/n'), a,{Q'/n)) + O (iV") 



n 



17) 



wliere 



g{as{Q'), «.(Q7^)) = lim (7(as(g'), a.W'K), «s(Q7^), e) (88) 

e— ►0 

is a generic Taylor series of its arguments. 

Renormalization group invariance thus implies that the physical anomalous dimen- 
sion 7 Eq. ( 1771) depends on its three arguments Q^, Q'^/N and Q'^/N'^ only through 
as- Clearly, any function of and can be expressed as a function of q;s(<5^) and 
as{Q'^/N) or asiQ"^ /N"^)- The nontrivial statement, which endows Eq. (IHTj) with pre- 
dictive power, is that the log derivative of 7, giasiQ"^), C(s{Q^ /n^), C(s{Q'^/nj) Eq. flHUj) . 
is analytic in its three arguments. This immediately implies that when 7 is computed 
at (fixed) order a^, it is a polynomial in In of k-th order at most. 

In order to discuss the factorization properties of our result we write the function 
g as 

^(a,(g2),a,(QVn'),«.(QV^)) = 9iMQ'),asiQVn)) 

+ ^?2(«.(Q'),«s(QV^')) 

+ gsiasiQ'),asiQyn),asiQyn')), 

(89) 

where 

00 00 

gi{as{Q'),as{Q'/n)) = Y.Y.3m0v<{Q^)<{Q' /n) (90) 

m=0 p=l 

00 00 



m=0 n=l 



00 00 00 



(91) 



m=0 n=l p=l 



xa:{Q'/n')al{Q'/n). (92) 

The dependence on the resummation variables Q^, Q'^/N and Q'^/N'^ is fully factorized 
if the bare coefficient functions has the factorized structure 



This is argued to be the case in the approach of Refs. [T^ fT2] . If this happens, 
the resummed anomalous dimension is given by Eq. (IHTl) with all Qmnp = except 

dr) dn 
7(iV,a,(Q2))= / _^^(o,a,(QVn))+ / ^ g,{0^ a.iQ'/n')). (94) 
Jl n Jl n 

We recall that the coefficient function depends on the parton subprocess in which the 
incoming partons are a,b (compare Eq. ([7j)). So, the factorization Eq. (lU^ applies 
to the coefficient function corresponding to each subprocess, and the decomposition 
Eq. (19^ to the physical anomalous dimension computed from each of these coefficient 
functions. 

A weaker form of factorization is obtained assuming that in the soft limit the 
A^-dependent and A^-independent parts of the coefficient function factorize: 

CW(Ar, Q', e) = C^''^\Q\ ao, e) C^'^^^Q'/N', Q'/N, ao, e). (95) 

This condition turns out to be satisfied [H] in Drell-Yan and deep-inelastic scattering 
to order a^. It holds in QED to all orders ^3] as a consequence of the fact that each 
emission in the soft limit can be described by universal (eikonal) factors, independent 
of the underlying diagram. This eikonal structure of Sudakov radiation has been 
argued in Refs. [21 HI] to apply also to QCD. If the factorized form Eq. ( 195|) holds, 
the coefficients Qmnp Eqs- (190119 If 921) vanish for all m 7^ 0, and the physical anomalous 
dimension takes the form 

/■^ dn dn 
7(Ar,a,(Q2)) = / _ ^^(q, a,(QVn)) + / — ^^(0, a,(QVn^)) 

Jl ^ Jl 
dn 

+ —gsiO,asiQ'/n'),a.iQyn)). (96) 
Jl 

It is interesting to observe that in the approach of Refs. [T^ [T2] for processes where 
more than one colour structure contributes to the cross-section, the factorization 
Eq. flUHl) of the coefficient function is argued to take place separately for each colour 
structure. This means that in such case the exponentiation takes place for each colour 
structure independently, i.e. the resummed cross section for each parton subprocess 
is in turn expressed as a sum of factorized terms of the form of Eq. ( 1931) . This 
happens for instance in the case of heavy quark production [El El]. In prompt 
photon production different colour structures appear for the gluon-gluon subprocess 
which starts at next-to-next-to-leading order, hence their separated exponentiation 
would be relevant for next-to-next-to-leading log resummed results. 

When several colour structures contribute to a given parton subprocess, the coef- 
ficients of the perturbative expansion Eq. fl76l) for that process take the form 

clti^)=ciS)i^) + cTie), (97) 

(assuming for definiteness that a colour singlet and octet contribution are present) 
so that the coefficient function can be written as a sum C(o) = + The 



argument which leads from Eq. (1761) to the resummed result Eq. ( IHTl) then implies 
that exponentiation takes place for each colour structure independently if and only if 

71 = 9 In /d In 78 = 5 In Cf^ /9 In (gg) 

are separately finite. 

This, however, is clearly a more restrictive assumption than that under which 
we have derived the result Eq. (l87|l . namely that the full anomalous dimension 7 is 
finite. It follows that exponentiation of each colour structure must be a special case 
of our result. However, this can only be true if the coefficients gijk of the expansion 
Eq. (iQOll of the physical anomalous dimension satisfy suitable relations. In particular, 
at the leading log level, it is easy to see that exponentiation of each colour structure 
is compatible with exponentiation of their sum only if the leading order coefficients 
are the same for the given colour structures: ^^qq^^ = g^^^ and ^^q^^q = (^q^q. This is 
indeed the case for heavy quark production (where (7001 = 0). 

Note that, however, if the factorization holds for each colour structure separately 
it will not apply to the sum of colour structures. For instance, the weaker form of 
factorization Eq. flU3|) requires that 

C\l,{e)=F,^k'{t)Gi^kMt), (99) 

but 

Fl^^,{t)GU_^,{t) + F,V(e)Gf_,_,,(e) ^ F,+,,(e)Gz„,-,.(e). (100) 

Hence, our result Eq. fl87|l for the sum of colour structures is more general than the 
separate exponentiation of individual colour structures, but it leads to results which 
have weaker factorization properties. 



4.5 Comparison with previous results 

In this section, we want to make a comparison with the resummation formula for 
prompt photon production previously released. In order to do this, we need to rewrite 
the NLL result of Ref. [H] in our formalism. The physical anomalous dimension can 
be obtained performing the Q^-logarithmic derivative of the NLL resummed exponent 
in the MS scheme of Ref. [H] . We obtain for a particular partonic sub-process: 

7(iV,a,(g2)) = / dx— [g^as{Q\l-xf)+g'^al{Q\l-xf) 

Jo i — X 

+^ia,(Q2(i _ x)) + g[a^AQ'i^ - ^-))] , (101) 



where 
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(104) 
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Here, A\ is the coefficient of ln(l/A^) in the Melhn transform of the Paa Altarelh-Parisi 
sphtting function at order a*, (3q is the a1 coefficient of the (3 function (Eq. (!HHl) in 



section 11.21) and B\ is a constant to be determined from the comparison with the 
fixed-order calculation. In Eqs. (1102111051) a, h are the incoming partons (on which the 
coefficient function implicitly depends) and d is the LO outgoing parton uniquely de- 
termined by the incoming ones. For completeness, we list explicitly these coefficients: 



a,=q,q 



(1) 



^a=g 



K 



67 



B 



d=q,q 
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where Nf is the number of active flavors. Now, performing the change of variable 



1 



in the integral Eq. 



n = 

1 — X 

we obtain at NLL 
^ dn 



(109) 



n 



[fi'oioa. 



+ fi'OOltts 



n 



+^020"^ ( ^ 



9QmOis 



«!)], 



n 



where 



9oio 



fl'ooi 



-92, 



-91, 



9020 



9002 



92 



9i 



JE^O 

2n 
Ait 



92 



-91 



(110) 

(111) 
(112) 



with 'Je the Euler constant. Hence, according to Ref.[Tl], we know exactly the value 
of the coefficients goio, go2o, fl'ooi and (7002- This enables us to compute predictions of 
high-order logarithmic contributions to the physical anomalous dimension performing 
a fixed order expansion of 7. 



We shall now show that the resummation formula of Ref.[Tl] predicts the coef- 
ficient of ln^(l/A^) of the fixed order expansion of 7, while in our approach it is 
required in order to perform a NLL resummation. We need to expand Eq. flllOp to 
order and this is obtained using the change of variable 

dn^_da^{^^ a = 1,2 (113) 
n apOs 

to perform the integral and expanding the two loops running of ag (see Eq. (1431) in 
section 11.21) . We find 



7 = [-(^001 +^oio)]«s(<5^)ln^ 

+ [-(5-002 + 5-020) - (^i/&o)(5ooi + 5oio)]a^(Q^)ln-^ 



+ 



^(5001 + 25010) 



aliQ^) In^ 1 + [-(^7003 + 5o3o)]«s(Q') In ^ 



+ [(36i/2)((7ooi + 2(7oio) + &o(5oo2 + 2go2oM{Q') In^ ^ 

+ [-{bl/3){gooi + 45oio)]«.'ln=^ ^ + 0{at). (114) 

In our approach, in order to determine the NLL resummation coefficients goio, go2o, 
gooi and goo2, we must compare this expansion to a fixed order computation of the 
physical anomalous dimension, which in the general has the form 

k i ^ 

i=i j=i 

where k is the fixed-order at which it has been computed (see Chapter [7] for a gen- 
eral discussion about the determination of the resummation coefficients). Hence, we 
determine the 4 NLL resummation coefficients through the following 4 independent 
conditions: 

5001 + 5010 = -ll (116) 

5001 + 25-010 = 7-72 (117) 

Oo 

5002 + 5020 = -7? - 77(5001 + 5oio) (118) 

Oo 

X 3& 

5002 + 2^020 = 7-72 - TTT" (5001 + 25010)- (119) 

Oo ZOq 

Thus, according to our formalism, all the coefficients 7J, 7I, 7^ and 7I should be 
known. The first three coefficients are all known thanks to the explicit 0{a'^) calcu- 
lation of the prompt photon cross section [551 EEl [57]. The last one (7I), is not yet 
known from explicit 0{al) calculation, but, according to the approach of Ref. [Hj, it 



is predicted using Eqs. ffTTMTTTImD : 
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The correctness of this result could be tested by an order calculation. If it were to 
fail, the more general resummation formula with (yfo2o determined by Eq. flll9l) should 
be used, or one of the resummations which do not assume the factorization Eq.( 



Chapter 5 



Resummation of rapidity 
distributions 



In this chapter, we present a derivation of the threshold resummation formula for the 
Drell-Yan and the prompt photon production rapidity distributions. Our arguments 
are vahd for all values of rapidity and to all orders in perturbative QCD. For the case 
of the Drell-Yan process, resummation is realized in a universal way, i.e. both for the 
production of a virtual photon 7* and the production of a vector boson W^, Z^. We 
will show that for the fixed-target proton-proton Drell-Yan experiment E866/NuSea 
used in current parton fits, the NLL resummation corrections are comparable to NLO 
fixed-order corrections and are crucial to obtain agreement with the data. This means 
that the NLL resummation of rapidity distributions is necessary and turns out to give 
better results than high-fixed-order calculations. We consider first the resummation 
of the DY case and its phenomenology and then the resummation of the prompt 
photon case. 

5.1 Threshold resummation of DY rapidity distri- 
butions 

5.1.1 General kinematics of Drell-Yan rapidity distributions 

We consider the general Drell-Yan process in which the collisions of two hadrons 
{Hi and H2) produce a virtual photon 7* (or an on-shell vector boson V) and any 
collection of hadrons (X): 



In particular, we are interested in the differential cross section -^^^{x, Q'^,Y), where 
is the invariant mass of the photon or of the vector boson, x is defined as usual 
as the fraction of invariant mass that the hadrons transfer to the photon (or to the 
vector boson) and Y is the rapidity of 7* {V) in the hadronic centre-of-mass- frame: 



Hi{Pi) + H2{P2) ^ Y{V){Q)+X{K). 



(1) 
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where E and pz are the energy and the longitudinal momentum of 7* {V) respectively. 
In this frame, the four- vector Q of 'y*{V) can be written in terms of its rapidity and 
its tranverse momentum 



Q = (Q°,Q) = {^JQ^ + QlcoshY,Q^,^^/Q^ + QlsmhY), (3) 
or in terms of the scattering angle 6 



Q = Q) = {\/Q' + \Q\\ Q±, \Q\ cos{e)). (4) 

For completeness, we recall that when = (which is not our case), the rapidity Y 
defined in Eq.(l2]) reduces to the pseudorapidity t] according to Eq.(jll): 

r/ = -ln(tan(e/2)). (5) 

At the partonic level, a parton 1(2) in the hadron Hi {H2) carries a fraction of 
momentum Xi {x2)'- 

j9l = XiPi = xi — (1, 0_L, 1)), P2 = X2P2 = X2 — {l,0±,-l)). (6) 

It is clear that the hadronic center-of-mass frame does not coincide with the partonic 
one, because xi is, in general, different from X2- Furthermore, from Eq.(l2]), we see 
that the rapidity is not an invariant. Hence, in order to define the rapidity in the 
partonic center-of-mass frame (y), we have to perform a boost of Y which connects 
the two frames. This provides us a relation between the rapidity in these two frames: 

y = y-lH^)- (7) 

A X2 

In order to understand the kinematic configurations in terms of rapidity, it is 
convenient to define a new variable u, 

Q-P2 X2 

With no partons radiated as in the case of the LO, the rapidity is obviously zero. 
Beyond the LO, one or more partons can be radiated. Now, if these partons are 
radiated coUinear with the incoming parton 2, then the partonic rapidity reaches its 
maximum value and u its minimum one. Similarly the minimum value of y and the 
maximum value of u is achieved when the radiated partons are collinear with the 
incoming parton 1. To be more precise, suppose that in the first case the radiated 
partons (collinear with the parton 2) carry away a fraction of momentum equal to 
(1 — z)p2, so that by momentum conservation Q = pi + zp2. In this case, we obtain 
immediately the lower bound for u, which is z. In the second case the collinear 
radiated partons have momentum (1 — z)pi, hence Q = zpi +p2 and the upper bound 
of u is So, z can be interpreted as the fraction of invariant mass that incoming 
partons transfer to 'y*{V). In fact: 

Q'^ X 



where we have neglected the quark masses. Therefore, we have that the kinematic 
constraints of u are: 



1 . N 

z<u<-. (10) 

z 

Then, since Xi <1 and X2 < 2, the lower and upper bounds of z are: 

x<z<l. (11) 

Thanks to Eqs. flTfSj) . the first relation can be translated directly into a relation for 
the upper and lower limit of the partonic center-of-mass rapidity: 

1 , 1,1 

-\nz <y < - \n-. (12) 

Now, we need to obtain the boundaries of the hadronic center-of-mass rapidity. Sub- 
stituting Eqs. (El ED into the two conditions u > z and -u < 1 we obtain the lower 
kinematical bound for xi and X2'- 

Xi > \fxe^ = x?, X2 > \fxe~^ = X2 (13) 

and the obvious requirment that x^2) ^ 1 implies that the hadronic rapidity has a 
lower and an upper bound: 

-lnx<F<-ln-. (14) 

2 - - 2 X ^ ' 

5.1.2 The universality of resummation in Drell-Yan processes 

According to standard factorization of collinear singularities of perturbative QCD, 
the expression for the hadronic differential cross section in rapidity has the form. 



da 



dQ^dY 



J2 [ dxi I dx2Fl''{xi,fi^)Fj''{x2,fi' 



■)2 



da^j f Q / 2\ 



where y depends on Y, xi and X2 according to Eq.([7]). The sum runs over all possible 
partonic subprocesses, Pj^^^ ,Fj'^^ are respectively the parton densities of the hadron 
Hi and H2, /i is the factorization scale (chosen equal to renormalization scale for 
simplicity) and daij / {dQ'^dy) is the partonic cross section. Even if the cross section 
Eq.( |T5i) is /i^-independent, this is not the case for each parton subprocess. However, 
the /i^-dependence of each contribution is proportional to the off-diagonal anoma- 
lous dimensions (or splitting functions), which in the threshold limit, {z 1) are 
suppressed by factors of 1 — z. Therefore, each partonic subprocess can be treated 
independently and is separately renormalizat ion-group invariant. Furthermore, the 
suppression, in the threshold limit, of the off-diagonal splitting functions implies also 



that only the gluon-quark channels are suppressed. So, in order to study resumma- 
tion, we will consider only the quark- anti-quark channel, which can be related to 
the same dimensionless coefficient function C{z,Q'^/ fi'^,as{fi^),y) for both, the pro- 
duction of a virtual photon and the production of a on-shell vector boson. In fact, if 
for the production of a virtual photon, we define C{z,Q'^/^'^,as{fi'^),y) through the 
equation, 

where the prefactor X1X2 has been introduced for future convenience, we find that for 
the case of the production of a real vector boson. 



X c(^z,^,a,{fi'),yy (17) 

where is the Fermi constant. My is the mass of the produced vector boson. The 
coefficients Cqq', for the different Drell-Yan processes, are given by: 

Cqq' = Ql^qq for 7* , (18) 
I 2 t^^ TJ/± 



Cqq' 



\Vqq,\' ioiW^, (19) 

^M? + {9l?%, for^". (20) 



Here, is the square charge of the quark g, Vqqi are the CKM mixing factors for the 
quark flavors g, q' and 

14 1 
9l = 2 ~ 3 ^^^^ ^1^' = 2 up-type quark, (21) 

12 1 
gy = ~2 3 3l = "2 a down-type quark, (22) 

with 6w the Weinberg weak mixing angle. As a consequence of these facts, resumma- 
tion has to be performed only for the quark- ant i- quark channels omitting the over- 
all dimensional factors of C(2;, Q^//^^, as(/i^), y) in the different Drell-Yan processes. 
Thus, we are left with the following dimensionless cross section, which has the form: 



a{x,Q\Y) ^ r^Ff^(xi,;.^)Ff^(x„^^ 



X C(^z,^,a,(/i2),yJ , (23) 

where Fi and F2 are quark or anti-quark parton densities in the hadron Hi and H2 
respectively. This shows the universality of resummation in Drell-Yan processes in 
the sense that only the renormalization-group invariant quantity defined in Eq.( |23l) 
has to be resummed. 



5.1.3 Factorization properties and the Mellin-Fourier trans- 
form 



For the case of the rapidity-integrated cross section, resummation is usually done in 
Mellin space transforming the variable x into its conjugate variable A^, because the 
Mellin transformation turns convolution products into ordinary products. Further- 
more, the Mellin space is the natural space where to define resummation of leading, 
next-to-leading and so on logarithmic contribution, because in this space momentum 
conservation is respected as shown in [58] ■ In the case of the rapidity distribution, 
the Mellin transformation is not sufficient. In fact, rewriting Eq.( |23ll in this form 



a{x,Q^,Y) = [ dxidx2dzF^'{xi,fi'^)F^^{x2,n^) 
Jo 



Q 



2 



X C yz,—,as{iJ^ ),y j 6{x - X1X2Z), (24) 
we see that the Mellin transform with respect to x, 

a{N,Q^,Y)= [ dxx^~^a{x,Q^,Y), (25) 



does not diagonalize the triple integral in Eg. (1241) . This is due to the fact that 
the partonic center-of-mass rapidity y depends on xi and X2 through Eq.(I71). The 
ordinary product in Mellin space can be recovered performing the Mellin transform 
with respect to x of the Fourier transform of Eq. fl2^ with respect to Y. Calling the 
Fourier moments M, using Eq.(I71) the relations ( fT2|T4l) and the identity 

_„.(,),y__,„_j = 1^ iyc{^.,- 



(26) 



we find that 



a{N,Q^,M) = I dxx^~^ [ dYe'^^^a{x,Q\Y) (27) 

Jo Jin 

= F^'{N + iM/2, n^)F^''{N - iM/2, jj^) 

X c(^N,^,asifi'),My (28) 



where 



Fl'^N ±iM/2,^?) = f dxx''-^^"''/^F^\x,^i^), (29) 

Jo 

c(N,'^,as{fi'),M] = dzz""-' dye^'^y 

V J Jo Jln^ 



'In 

X C[z,^,asifi'),y]. (30) 



Eg. (1281) shows that performing the Mellin- Fourier moments of the hadronic dimen- 
sionless cross section Eg. (1231) . we recover an ordinary product of the Mellin- Fourier 
trnsform of the coefficient function and the Mellin moments of the parton densities 
translated outside the real axis by ±zM/2. Because the coefficient function is sym- 
metric in y, we can rewrite Eg.( !30l) in this way: 



q2 \ /-Inl/v^ 



C [N,^,as{^l^),Mj = 2j dzz""-' J dycos{My) 



X c(^z,^,asifi'),yy (31) 

From this last eguation and Eg.(l29l). we see that the dependence on M, the Fourier 
conjugate of the rapidity y, originates from the parton densities, that depend on 
N =1= iM/2, and from the factor of cos{My) in the integrand of Eg. (131 p . 

5.1.4 The all-order resummation formula and its NLL imple- 
mentation 



In this section, we show that the resummed expression of Eg. flzSj) is obtained by 

simply replacing the coefficient function C (^N, ^,as(/i^),Mj with its integral over 

y, resummed to the desired logarithmic accuracy. This is eguivalent to saying that 
the factor of cos{My) in Eg. (l3ip is irrelevant in the large-N limit. Indeed, one can 
expand cos (My) in powers of y, 



2„,2 



cos{My) = 1 - ^ + 0(MV)- (32) 

and observe that the first term of this expansion leads to a convergent integral (the 
rapidity- integrated cross section), while the subseguent terms are suppressed by pow- 
ers of (1 — ^), since the upper integration bound in Eg. (|3Tl) is 

lni= = i(l-.) + 0((l-z)2). (33) 



Hence, up to terms suppressed by factors 1/A^, Eg. (l30l) is egual to the Mellin transform 
of the rapidity-integrated Drell-Yan coefficient function that we call Ci{N, Q^/ /i^, as(/i^)). 
This completes our proof. We get 



a''''{N,Q'',M) = F^'{N + iM/2, i/)F^'{N -iM/2, i/'' 

This theorethical result is very important: it shows that, near threshold, the 
Mellin-Fourier transform of the coefficient function does not depend on the Fourier 
moments and that this is valid to all orders of QCD perturbation theory. Furthermore 
this result remains valid for all values of hadronic center-of-mass rapidity, because we 



X Cr(N,%,as{fi')). (34) 



have introduced a suitable integral transform over rapidity. The resummed rapidity- 
integrated Drell-Yan coefficient function to NLL order has been studied in Section 
13.41 It is given by 

Cr (n, as{^^'')] = [exp{lniV^i(A, 2) + g^iX, 2)}]^.^^. (35) 



where A = boasi^f) InN and where the resummation functions 5'i(A,2) and g2{X,2) 



are given in Eqs. fll01lll02p of Section \TM with the resummation coefficients in the MS 
scheme given in Eq.(l46l) of Section [2.3. 1[ 

Now, we want to arrive to a NLO and NLL expression of the rapidity-dependent 
dimensional cross section. This is achieved ffistly taking the Mellin and Fourier 
inverse transforms of (t^'^^{N, Q^, M) Eq. flM|) in order to turn back to the variables x 
and Y: 

r°° rIM rC+ioo <nT 

a-{x, Q\ Y) = / ^e-^*^^ / |^ x-^^a^iN, Q\ M). (36) 

In principle the contour in the complex A^-space of the inverse Mellin transform in 
Eq. flHUl) has to be chosen in such a way that the intersection of C with the real axis lies 
to the right of all the singularities of the integrand. In practice, this is not possible, 
because the resummed coefficient function Eqs.f ll011)102p of section [33] has a branch 
cut on the real positive axis for 



N > Nl = e-^o-^c?^), (37) 

which corresponds to the Landau singularity of as{Q'^ /N"-) (see Eq. (l4T]) in Appendix 
IL2p . This is due to the fact that if the A^-space expression is expanded in powers of a^, 
and the Mellin inversion is performed order by order, a divergent series is obtained. 
The "Minimal Prescription" proposed in [SB] gives a well defined formula to obtain the 
resummed result in x-space to which the divergent series is asymptotic and is simply 
obtained choosing C = Cmp in such a way that all the poles of the integrand are to 
the left, except the Landau pole Eq.(l37l). Recently, another method has been prposed 
in Ref.[nn]- Here, we will adopt the "Minimal Prescription " formula, deforming the 
contour in order to improve numerical convergence and to avoid the singularities of 
the parton densities of Eq.( |MI) which are transated out of the real axis by ±2M/2. 
Hence, we perform the A^-integral in Eq.( |36l) over a curve P given by: 



P = Ti + Pa + Ps (38) 
M 

ri(t) = CMP~i— + t{l + i), tG(-cx),0) (39) 
M 

T2{s) = Cmp + is—, sG(-1,1) (40) 
M 

T-i{t) = CMP + ^Y-t(l-0, tG(0,+oo) (41) 

The double inverse tranform of Eq. ( l36l) over the curve P then becomes: 

1 dm 

a''Hx,Q^X) = - cod-Y\nm)a''\x,Q\-\nm), (42) 

vr ./n m 



where we have done the change of variableM = — Inm. The factor a'^'^^{x, Q^, M) of 
the integrand in Eq. (l42l) is given by 

r/s r 

- .-CMP-lns+m/2+l)^res^^^^ + 111 S - i{M/2 + 1), Q\ M) 



(43) 



1 



TT 



-3? 



X 



x(l - z) + ■^x-^^'^-''^'/^a''\CMP + M) 



where we have done another change of variables {t = —Ins). Eqs. fl42|43l) are the 
expressions that we use to evaluate numerically the resummed adimensional cross 
section in the variables x and Y Eq. flHUj) . Furthermore, we need to know the analytic 



continuations to all the complex plane of the parton densities at the scale in 
Eq.l lMI) . Here, we need to evolve up a partonic fit taken at a certain scale solving the 
DGLAP evolution equations in Mellin space. The solution of the evolution equations 
is given in Section 11.51 

Finally, we want to obtain a NLO determination of the cross section improved with 
NLL resummation. In order to do this, we must keep the resummed dimensionaless 
part of the cross section Eq. fl42l) . multiply it by the correct dimensional prefactors 
Eqs.f l57l[T9|20l) and parton densities, add to the resummed part the full NLO cross 
section and subtract the double-counted logarithmic enhanced contributions. Thus, 
we have 

da da^^^ da'^' 



dQ^dY dQ'^dY ^ dQ^dY 



' da""^' ' 


- as 


' d 


/ da''^' Y 




_dQW_ 


das 







(44) 



The first term is the full NLO cross section given in [T71 EHl ED E2]- We report the 
complete expression in Appendix O 

The third and the fourth terms in Eq. (1441) are obtained in the same way as the 
second one, but with the substitutions 



CY 



N, ^,a<,(/i^; 
^, %,as(/i^; 



(45) 



a,(/i2)2Ai <^ In^iV + lniV 



27i 



respectively. The terms that appear in the second in the second line of Eqs. (l45l) are 
exactly the 0{as) logarithmic enhanced contributions in the MS scheme. 

We note that the resummed cross section Eq. (l44l) is relevant even when the variable 
X is not large. In fact, the cross section can get the dominant contributions from the 
integral in Eq. fl23p for values of z Eq.(^ that are near the threshold even when x is 
not close to one, because of the strong suppression of parton densities Fi{xi, /i^) when 
Xi are large. 



5.1.5 NLL impact of resummation at E866 experiment 

To show the importance of this resummation, we have calculated the Drell-Yan rapid- 
ity distribution for proton-proton collisions at the Fermilab fixed-target experiment 



Figure 5.1: Y-dependence of (P a / [dQ"^ dY) in units of pb/GeV^. The curves are, from top to 
bottom, the NLO result (red band), the LO+LL resummation (blue band) and the LO (black 
band). The bands are obtained varying the factorization scale between /i^ — 2Q^ and /i^ — 1/2Q^. 



E866/NuSea [3H]- The center-of-mass energy has been fixed at \fS = 38.76 GeV and 
the invariant mass of the virtual photon 7* has been chosen to be = 64 GeV^ in 
analogy with [19]. Clearly the contribution of the virtual can be neglected, because 
its mass is much bigger than Q^. In this case x = 0.04260 and the upper and lower 
bound of the hadronic rapidity Y Eq.( |T4l) are given by ±1.57795. We have evolved up 
the MRST 2001 parton distributions (taken at /i^ = 1 GeV^) in order to compare to 
Ref. |19]. where the NNLO calculation is performed. However, results obtained using 
more modern parton sets should not be very different. The LO parton set is given in 
p] with af^{mz) = 0.130 and the NLO set is given in [6l] with a^^°{mz) = 0.119. 
The evolution of parton densities at the scale yU^ has been performed in the variable 
flavor number scheme. The quarks has been considered massless and, at the scale 
of the transition of the flavor number {Nf —>■ Nf + 1), the new flavor is generated 
dynamically. The resummation formula Eq. flMl) together with Eqs. fll01H35l) has been 
used with the number of flavors Nf = 4. 

In figure Em we plot the rapidity-dependence of the cross section at LO, NLO and 
LO improved with LL resummation. The effect of LL resummation is small compared 
to the effect of the full NLO correction. We see that, at leading order, the impact of 
the resummation is negligible in comparison to the NLO fixed-order correction. This 
means that, at leading order, resummation is not necessary. 

The LO, the NLO and its NLL improvement cross sections are shown in figure 
15.21 The effect of the NLL resummation in the central rapidity region is almost as 
large as the NLO correction, but it reduces the cross section instead of enhancing it 
for not large values of rapidity. The origin of this suppression will be discussed in 
the next Section. Going from the LO result to the NLO with NLL resummation, we 
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Figure 5.2: Y-dependence of d?a/{dQ'^dY) in units of pb/GcV^. The curves are, from top to 
bottom, the NLO result (red band), the NLO+NLL rcsmnmation (green band) and the LO (black 
band). The bands are obtained as in figure [5T] 



note a reduction of the dependence on the factorization scale i.e. a reduction of the 
theoretical error. 

Now, we want to establish if the leading logarithmic terms that are included in 
the resummed exptonent represent a good approximation to the exact fixed order 
computation. Only if this is the case, we can believe that our resummation is reliable 
in perturbative QCD. In order to do this, we compare the full NLO DY rapidity 
cross section with the one obtained including only the large- leading terms of the 
coefficient function. For simplicity, we choose the factorization scale yU^ to the scale 
of the process . The leading large- coefficient function is given by: 



C^^^'^ (iV, a,(g2)) = 1 + a,(Q2)2Ai In^ AT + 27^^ In iV + 7I - 2 + — , (46) 



where we have added the constant terms at 0{as) which are not resummed. For an 
explicit derivation of these constant terms see for example Ref. [17] Section 3. We 
plot the result in Figure 15. 3[ We see that the leading terms Eq. fH6l) represent a good 
approximation to the exact NLO computation, because they account for more than 
90% of the full NLO computation for all relevant rapidities. 

In Figure 15.41 we plot only the 0{as) correction. Here we see that the 0{as) 
contribution of Eq.( H6l) represent a good approximation to the exact 0{as) NLO 
contribution, because it accounts for more than 80% of the full 0{as) NLO correction 
for all relevant rapidities. 

In figure [^31 we plot the experimental data of Ref. [SH] converted to the Y variable 
together with our NLO and NLL resummed predictions. 



The data in Ref. [38j are tabulated in invariant Drell-Yan pair mass yQ^ and 

Feynman xp bins. To convert the data to the hadronic rapidity F, we have used the 





Figure 5.3: (Pa/{dQ'^dY) in units of pb/GcV^ for the full NLO computation (upper red line) 
and for the leading terms of Eq. (|lB|) (lower black line). It has been calculated for one value of the 
factorization scale fj.^ ~ Q^. 




Figure 5.4: d?a /{dQ'^dY) in units of pb/GeV^ for 0{as) correction of the fuU NLO computation 
(upper red line) and for the leading 0{as) term of Eq. lH^l) (lower black line). 



Figure 5.5: Y dependence of (P a / [dQ"^ dY) in units of pb/GeV^. The curves are, from top to 
bottom, the NLO resuh (red band) and the NLO+NLL resummation (green band) together with 
the E866/NuSea data. The bands are obtained as in figure [5Tl 



definition of tlie Feynman xp wliicli is 



wliere we liave used Eq.(j3]). Solving Eg. ( 1471) in Y we liave 



(47) 



F = In 



H 



Vs 



(4J 



Witli tliis equation and witli tlie aid of tlie Q± distribution, which is also given 
in Ref.[3H], we have converted the data from xp to Y. Furthermore, for each xp 
bin, we have done the weighted average of three \/Q^ bins (7.2 < ^/Q^ < 7.7; 
7.7 < < 8 2 and 8.2 < ^/Q^ < 8.7 with the energies in GeV). 

The agreement with data is good and a great improvement for not large rapidity 
is obtained with respect to the NLO calculation. We note also that the NLL resum- 
mation gives better result than the NNLO calculation performed in [19]. The NNLO 
prediction has a worse agreement with data than the NLO one for not large values 
of rapidity. This result suggests that, for the case of rapidity distributions, NLL re- 
summation is more important than high-fixed-order calculation and that it can be so 
even at higher center-of-mass energies. 



5.2 The origin of suppression 

In this section, we shall show that the suppression of the cross section of the NLL 
correction with the parameter choices of the experiment E866 is due to the shift in 



the complex plane of the dominant contribution of the resummed exponent. We shall 
do it using a simplified toy-model. 

Consider the collision of only two quarks with parton density 

F(x) = (l-x)l (49) 

Its Mellin transform is given by: 

Furthermore, we take the double-log approximation (DLA) which is obtained per- 
forming the limit A — * in the resummed exponent Eq. flH^j) . Thus, in this simple 
model, the Mellin- Fourier transform of the NLL resummed cross section Eq. (1441) can 
be written in the following form: 

(r{N, M) = a^^{N, M) + \F{N + iM /2)\^ ^a^^^{N), (51) 
where a^^{N,M) are the exact NLO Mellin- Fourier moments and where 

A(T°L^(Ar) = [e"^2Ai \n^N ^^2Ai In^ A^j . (52) 

If there is a suppression, this means that the quantity 

aiN, M) - a^^iN, M) = r^^-- ^jj-^^ (53) 

^ ' ^ ' ^ [iV2 + M^] [(iV + l)2 + Mi] [(iv + 2)2 + Mi]' ^ ^ 



should produce a negative contribution in performing the inverse Mellin and Fourier 
transform. It is given by the integral 



2vr 7c-.oo 27r^ [N^ + Mi] [(N + 1)2 + Mi] [(jy + 2)2 + Mi] 

The integrand function of this expression has not only a cut on the negative real axis 
(as it happens in the inclusive case), but also poles that are shifted in the complex 
plane: 

iM , , 

-n±— ; n = 0,1,2. (55) 

Because of the factor in the inverse Mellin integral in Eq. (1541) . its dominant 
contribution comes from the poles with n = in Eq. fl55l) . The contribution of the 
pole at +iM/2 is given by 

^^.^-.M(y+i„ ^) 4A^°"^(iV = ^ ^^^^ 



oo 



27r iM{iM + l){iM + 2)' 



where 



Aa^^^(A^ = iM/2) = exp 



as2Ai In^ ^ — ^ + ivr In ^ — ^ 

' 2 4 2 

. , , 2 1^1 TT^ , IMI 

-l-as2Ai W^- — ^ + Z7rln^ — ^ 

'24 2 



+ (57) 



The important thing to notice of this contribution is the fact that the imaginary 
pole has produced an oscillating prefactor in front of the resummed exponent which, 
together with the oscillating factor of the Fourier inverse integral in Eg. (1561) at zero 
hadronic rapidity Y, is given by 



exp 



0^27: Ai In -iM\ny/x 



{51 



We note that, with the inclusion of the contribution of the other pole at —iM/2, the 
real and immaginary part of Eq. ( 157|) contribute to the integral of Eq. (154|) . The real 
and immaginary part of Eq. fl57l) are given by: 



3fJ[Aa^^^(^M/2)] = e"^^"H'" — ~-Vcos(a327rAiln^) (59) 



as2Ai ( In^ 



„2 m 



\M\ TT 



»[Aa°^^(«M/2)] = e^^'^n'-^ ^--;sin(a,27rAiln^ ) (60) 

|M| 



|M| 



—as2T[Ai In ■ 



2 



Now, to roughly estimate the effect of the oscillating factor in Eq. (!3U|) . we use the 
value of M = Mq where the phase of Eq. ( !58l) is stationary and is given by: 

as2'n'Ai , , 

Substituting this in Eqs. ( 159)1601) . we find a suppression of about 35% for the parameter 
choice of E866 experiment and a suppression of about 10% for the W boson production 
at RHIC with a center-of-mass energy of = 500GeV which has more or less the 
same value of x. We should now recall that this is a rough estimation and that there 
is also the contribution of the cut on the negative real axis which usually produces an 
enhancement. However, comparing this estimation with the result for the W boson 
production at RHIC (see e.g. figure 1 and 2 in reference [l7j), where resummation 
produces an enhancement of about 4% we can believe that the ignored contribution 
in this section produce an enhancement of about 15%, thus giving a suppressio of 
about 20% at E866 experiment. 

5.3 Resummation of prompt photon rapidity dis- 
tribution 

5.3.1 General kinematics of prompt photon rapidity distri- 
bution 

Here, we consider the rapidity distribution of the prompt photon process discussed 
in chapter m 

i7i(Pi) + H2{P2) ^ 7(P7) + ^- (62) 



Specifically, we are interested in the differential cross section p\ ^^^^y (^-L ? pI ? '^7) ; 
where as in section HIT] of chapter |l]p_L is the transverse momentum of the photon, r]^ 
is its hadronic center-of-mass pseudorapidity and 

x^ = —. (63) 

The pseudorapidity of the direct real photon in the partonic center-of-mass frame 7)^ 
is related to //-y through Eq.(l6]) in section WAk 

1 1 ^1 

^7 = ^7-oln— . (64) 

Furthermore, as in chapter HI we use the following parametrizations of the photon 
and of the incoming partons' momenta 

P7 = (p± coshi7^,pj_,p_L sinhr)^), (65) 

Vi = XiPi = ^(1,0^,1), (66) 

P2 = X2P2 = ^(l,0x,-l). (67) 

The transverse energy that the partons can transfer to the outgoing partons must be 
less than the partonic center-of-mass energy y/s = y/xiX2S. This means that 

z cosh^ < 1, (68) 

where we have defined the parton scaling variable 



X1X2 



(69) 



with, as in chapter IU = 4p^. Eg. (1681) implies that the upper and lower boundaries 
for the partonic center-of-mass pseudorapidity are given by 

V- <fi^<fi+, (70) 

where 

77±=ln(^± 1 =±ln(^ + (71) 



z \ z \ \ z V z 



Using Eq. ( l64l) , we can rewrite the transverse energy condition Eq. ( 168|) as a condition 
for the lower bound of X2- 



X2 > ^ V- — = 4- 72 



Now, the requirement that 0:2 < 1 implies the lower bound for Xi. 



The upper and lower bounds of the hadronic center-of-mass pseudorapidity can be 
found with the obvious condition that x^^) < 1- In this way, we find, 

V- <V^<V+, (74) 

where 

V± = In f-L ± J—-l] = ±ln + J—-l] ■ (T5) 



5.3.2 Mellin-Fourier transform and all-order resummation 

The expression with the factorization of colhnear singularities of this cross section in 
perturbative QCD is 

pI- — —{x±,pl,r]^) = / dxi dx2XiF^^{xi, fi^)x2F^''{x2, fi'^) 
dp±dr]^ ^ J^o J^o 

xCab (^z, j^, as(/i^), 17^^ 5{x± - ZX1X2), (76) 

where F^i(xi,/i^), F^^{x2,n'^) s-'^s the distribution functions of partons a, 6 in the 
colliding hadrons and //^ is the factorization scale equal to the renormalization scale. 
fj^ is a function of t]^, Xi and X2 as defined by Eq. fl64l) . The coefficient function 
Cab{z, as(/i^), fj-y) is defined in terms of the partonic cross section for the process 
where partons a, b are incoming as 



Q ^ 2^ ^ \ _ 3 da, 



C.,^.,^,a.(/.^),r),J=pl^. (77) 

To allow the Mellin transform to deconvolute Eq.( |76l) . we first perform the Fourier 
transform with respect to f]^, thus obtaining 

a{N,Q\M) = / dx^xl~^ / d7]^pl-——{x^,pl,7]^) (78) 
Jo Jr^_ dp^dr]^ 



F'^'iN + 1 + tM/2, ii'')F^\N + 1 - fi' 



a,b 

")2 



where 



xCab{N,^,as{i^'),M]. (79) 



F^^{N + l±iM/2,n^) = f dxx''^'''''^Fl'^{x,^?), (80) 



Cab [ N, ^, a,(/i2), = 2 ^ dzz"""' j^^ d% cos{M% 



X C[z,^,as{^l^),%]. (81) 



A resummed expression of Eg. (1831) in the threshold limit for the transverse energy 
[z —>■ 1 or equivalently — > oo) is obtained in the same way as we have done at the 
beginning of section [3.1.41 since the upper integration bound of fj^ in Eq. flHTj) is 




r)+ = lnl ^ + ) =,/T^-(l-z) + 0{{l-zf/'). (82) 



Thus, up to terms suppressed by factors the resummed expression of Eq.(l83l) is: 
a'''{N,Q^,M) = J]FfHAr + l+zM/2,/i2)F,^2^iV+l-zM/2,/i2) 

a,b 

xCr:,(^N,^,a,{f,')y (83) 

where CJ^| is the resummed pseudorapidity-integrated coefficient function for the 
prompt photon production for the subprocess which involves the initial partons a, b. 
These resummed coefficient function has been studied in chapter H] and in Refs. p31 
l65] . This result is analogous to that of the Drell-Yan rapidity distributions case, 
in the sense that the resummed formula is obtained through a translation of the 
parton densities' moments by 1 ± iM/2 and the pseudorapidity-integrated coefficient 
functions. 



Chapter 6 



Renormalization group 
resummation of tranverse 
distributions 



We prove the all-order exponentiation of soft logarithmic corrections at small trans- 
verse momentum to the distribution of Drell-Yan process. We apply the renormaliza- 
tion group approach developed in the context of integrated cross sections. We show 
that all large logs in the soft limit can be expressed in terms of a single dimensional 
variable, and wc use the renormalization group to rcsum them. The resummed result 
that we obtain is, beyond the next-to-leading log accuracy, more general and less pre- 
dictive than those previously released. The origin of this could be due to factorization 
properties of the cross section. The understanding of this point is a work in progress. 

6.1 Drell-Yan distribution at small transverse mo- 
mentum 

We consider the Drell-Yan process 



and, in particular, the differential cross section ^^y if^' gl, 3:2), where g_L is the 

transverse momentum with respect to colliding axis of the hadrons H\ and i72, is 
the virtuality of photon and X\^X2 are useful dimensionless variables, that, in terms 
of the hadronic center-of-mass squared energy S — {P\ -|- P-i)^ and the photon center- 
of-mass rapidity y, are given by 



The relation between these two variables and the fraction of energy carried by the 
virtual photon is 



//i(Pi) + 7/2(^2) ^ 7*(Q) + 



(1) 





(3) 



2 
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According to standard factorization of perturbative QCD, the expression for the dif- 
ferential cross section is 



da 



dqldY 



{Q^,q'i,Xi,X2) = dzi dz2fi{zi,n'^)f2{z2,n'^ 



dq\dy ' 

where /x^), /2(^2, /^^) are the parton distribution functions of the colliding quark 
and anti-quark in the hadrons Hi and H2 respectively. The arbitrary scale fi"^ is the 
factorization scale, which, for simplicity, is chosen to be equal to the renormalization 
scale. The condition that the invariant mass of the emitted particles cannot be 

2 

negative, imposes that {zi — Xi){z2 — X2) > ^ and, taking the small q'j_ limit, we 
obtain that 

= Xi, Z™" = X2. (5) 

The partonic center-of-mass squared energy s and rapidity y are related to the 
hadronic ones by a scaling and a boost along the collision axis with respect to the 
longitudinal momentum fraction Zi, Z2 of the incoming partons: 

s = ziZ2S\ y = Y-\\M—. (6) 

I Z2 

We define analogous variables to that of Eqs. (O]) at the partonic level 



Z\ \ S Z2 \ S 



with inverse relations 



y = o In 



{XxlZx){x2lz2y 2 0:2/^2' 

Now, thanks to these equations, we can define a dimensionless differential cross section 
and coefficient function 



ly(gi/Q^Xl,X2) = -^:T|^(Q',gLxi,X2), (9) 

X\X2 dq^dY 

(1^' Q"^' Zl Z2' " ) {xi/ zi){x2/ Z2) dq\dy 

(10) 

in such a way that Eq.(jl]), together with the conditions Eqs.(l5]), takes the useful form 
of a convolution product 



TT7 / Q]_ X2 I 2\\ tn,2 2 2 t 2\\ 



W{ql/Q\xi,X2) 



-_ f ^ f ^f^(^zi,f/)f2{z2,fl'^) 
Jxi Z\ ^2 

X W{(fl^^, ql/Q\ xi/zi, X2/Z2, (11) 

which is valid only for small q\. 



6.2 The role of standard factorization 



It is known that this expression (or equivalently Eq.(jl])) is originated by the factor- 
ization of colhnear divergences in the impact parameter (6) which is conjugate upon 
Fourier transformation to the transverse momentum {q±): 

W{Q%^,x,,X2) = J d^q±e''^'W{ql/Q^,Xi,X2). (12) 

In d = A — 2e dimensions this factorization has the form 

#(QV/x^QV,Xl,X2,a.(/x')) = 
' — —Ziz,, a,{^\ e)Z{z,, a,{^% 6)#W(Q^ h\ ^, ^, ao, e). (13) 

Note that the universal function Z that extracts the colhnear divergences from the 
bare coefficient function doesn't depend on the Fourier conjugate (6) of the transverse 
momentum Ref . [1 1] . In Fourier space Eg. (ITT]) , becomes 

-,27,2 \ '^^1 '^^2 r , 2\r/ 2\ 



iy(QV,Xi,X2) = / ^ / ^/i(^,,^^)/2(z2,^' 
Jx2 ^1 JX2 ^2 

X W^(QV/X^ Q%\ X2/^2, «s(/i')) (14) 

Furthermore, Eq. fllip tells us that the differential cross section is a convolutional 
product which is diagonalized by a double Mellin transform. Thus, performing the 
double Mellin and Fourier transform, the coefficient function takes the simple factor- 
ized form: 

W{Q%\N^,N2) = /l(iVl,/i2)/2(iV2,/x')Ty(gV/x^gV,iVl,iV2,a,(/i2)). (15) 

Our goal is to resum the large logarithms \n.Q%^ to all logarithmic orders. These logs 
are present to all orders in the contributions to this differential cross section. They 
come from the kinematical region of soft and colhnear emissions. However, we know 
from Eq.( |T3i) that colhnear divergences that arises in the limit g_L — are absorbed 
in parton distribution function evolution. Consequently, we will resum only the large 
logarithms \nQ%'^ that come from soft contributions. 

We define the usual physical anomalous dimension: 
Q''-^W{Q%\ iVi, N2) = iiw){Q^h\ iVi, iV2, as{Q^))W{Q%\ N,, N2). (16) 

It is clear that ■j(w){Q^b'^, ^i, N2,C(s{Q^)) is a renormalization group invariant and 
we will show that it is also independent of A^i and N2 when choosing the arbitrary 
scale /i^ equal to 1/6^ and taking into account only soft contributions. Thus, in the 
soft limit and with the convenient choice /x^ = 1/6^, we can write 

ifwriQ'b', N,, N2, a.m) = 7(1, Q'b', a.m) (17) 



So, the resummed expression for the cross section Eq.( fTT]) in Fourier space, in which 
the colhnear contributions to the large InQ^b^ are separated from the soft ones, has 
the general form: 

W-{Q'b\x„x,) = ['^ r^Mz„l/b')Mz„l/b')K-{b\Ql,Q') 

xW^-'^qy, Qlb\ xi/zi, «s(l/&')), (18) 

where 

K-{Q%\QlQ^) = exp | f ^r-{%,fi'b',a,m'))\ , (19) 

The scale Qq, must be larger than the lower limit of the perturbative analysis (Qq > 
Aq(jjj). Hence, in order to absorb the possible large correction of the type InQo^^ 
we will always choose Ql = 1/b^. Accordingly, the condition Qq > Aq^^ becomes 
6^ < 1/Aq(J2j and the resummed exponent of Eq. (fT9i) is related to the resummed 
physical anomalous dimension Y'^'^ through the logarithmic derivative: 

7-^(1, QV, a.(Q2)) =g2-^ r ^r-{%,fi'b^a,m')). (20) 

It is now clear that resummation of collinear emissions is realized by the parton 
distribution evolution thanks to the fact that the factorization scale yU^ is arbitrary. 
Resummation of soft gluon emissions can be achieved by the resummation of the ex- 
ponent that appears in this expression. This is the subject of the next section. 



6.3 The q\ singularities of soft gluon contributions 

We now proceed through the calculation of the resummed exponent using kinematics 
analysis and renormalization group improvement. According to the Appendix [B], the 
phase space measure in d = 4 — 2e dimensions for n extra emissions of the partonic 
Drell-Yan subprocess can be written for n = and n > 1 respectively as: 

dqldy = gl^(l- ^1)^(1 (21) 



dqj^dy Jq 27r 



e = M'- Ml = ^[(l-ei)(l-6)+gi(l-ei-6)], (22) 



where N{e) = l/(2(47r)2-2^), = Xi/zi, q\ = ql/Q^ and d(]"-^(e) stands for the 
integration of n — 1 dimensionless variables {zi, i = 1, . . . , n — 1). and ^2 are related 
to the partonic center-of-mass rapidity (y) and energy (s) by the relations: 



ei6 ' 



6 



(23) 



The phase space measure d(j)n{k; fci, . . . , is the same as the phase space measure 
of the DIS process with an incoming momentum with a nonzero invariant mass {k'^ = 
M ^) and n outgoing massless particles. This phase space has been analyzed in Section 
13.11 and is given by, 



d(f)i{k; ki) 
d(f)n{k] ki,..., kn) 



N{e) 



n = 



2tc 



2n 



n-l 



(M2)"-2-("-i)^(il]"-i(e), 



(24) 

n > 1. (25) 



According to this, we can rewrite Eqs. fl21l) .( !22l) in this form: 



d(pi{pi +P2\q) 
dq\dy 
#2(Pi +P2;q,ki) 



dq,dy 



Ml 



— 5(1-^1)5(1 -6)5(gl) 
{qlrN{e)6{Ml) 

^[(i-^i)(i-6) + g!(i-ei-6)]. 



(26) 
(27) 
(28) 



for n = 0, 1 respectively, and 

(pi +P2; q,ki,..., kn) 
dqidy 



{ql 



2-n 



(M2)"-2-(«-i)^rff]"-i^^^ ^29) 



for n>2. The dependence of the phase space on q]_ comes entirely from the factors: 



{qlr5{Ml), n 

{Mir-\qir%Mi) 



2\-(n-l)e-l 



n>l. 



(30) 
(31) 



The phase space measure must be combined with the square modulus of the 
amplitude, in order to determine the logarithmic singularities in = which are 
regularized in d = 4 — 2e dimensions. Studying the behavior of the invariants that 
can be constructed with the external momenta, we can establish in which kinematical 
region the square modulus of the amplitude can be singular in -^0. From the 
study of the DIS-like emissions (see Section 13. 2p we know that 



k^ 



,^n—l ' ' ' 



1)^/^(1 



1 < i < n 



(1 - 2„_i) 



A:?. 



(32) 

(33) 
(34) 



This means that all the invariants that can appear in the function Dg{i3,Pe) in 
Eq. (l22l) of Section 13.11 can be expressed in terms of the following ones: 



Pi ■ h 



P2 ■ ki 



p\=pI = = 0, 



Pi ■P2 



6^' 
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[(i-eO(i-6) + gl(i-6-6)] 



1/2 



(35) 
(36) 
(37) 
(38) 



In the case n = 1, the single emission squared amplitude at tree level has a.l/q\ 
singularity for g_L — >^ 0. This can be easily seen by an explicit 0{as) computation 
(see for example Refs.[66 t Wl \ 168]). Hence, in the general case, we expect that in the 
q± ^ Q limit the squared amplitude has the following behavior 

l^n+ir - ^(Mo^)"HMo^)""^?„,„,(ei,6), (39) 

where and k are integer or half-integer numbers (see Eqs.( !36|371) ). However, as 
discussed in Section I3.2[ here we will assume that only the integer powers of Mq 



contribute. Then, we know that phase space contributes the factors of Eqs.( l30]l3T|) 
and, hence, Eq. 0391) implies that a generic contribution to the coefficient function 
W^^^ has the folowing structure: 



n = 1, 
n > 1, 



(40) 
(41) 



where for the moment we do not care about the overall dimensional factor. Now, we 
are interested in taking into account only the 1/^i singularity, because more singular 
terms are forbidden and less singular ones are suppressed. As proven in Appendix [Dl 
in the limit — > 0, 



66 



(1-6)+ (1-e 



lngi5(l-^05(l-6) 



+ 0{qll (42) 



and, for 7] = —{n — n'2 — l)e (with e < n — — 1 > 0), 

[(i-ei)(i-6) + gi(i-6-6)r^ 



(43) 



;i-6)''-'(i-6 



-5(1-^1)5(1- 6) + 0(gi 



Therefore, ni 



-n + 1, n'n < n — 1. Furthermore, we note that the terms not 



proportional to 5(1 — 6)^(1 — 6) ^-^e divergent (in the q± —>■ limit) due to collinear 
emissions. These divergences are absorbed and resummed by the Altarelli-Parisi 
evolution of the parton distributions /i(zi, /x^) and f2{z2, /i^) when they are evaluated 



at the scale /x^ = 1/6^ in Fourier space (see Eq.(|T8ll). As a first conclusion, we obtain 
that the contributions that must be resummed in the limit — ^* are those which are 
proportional to 5(1 — ^i)5(l — ^2) and, thus, belong to the kinematical region of only 
the soft extra emissions. Hence, we obtain that the soft part of the coefficient function 
that must be resummed has, after the inclusion of loops, the following general form: 



W{Q', ql, 6, 6, ao, e) = 1] a^W^iQ^ ql, ao, e)6{l - 6)5(1 - 6), (44) 

n=0 



with 



W^n(g^gi,«o,6) = (Q^)^"1C7i;)(e)5(gi) + E^n?(e)(^l)-^-'=^ 

k=2 

n 



k=l 



where the factor of (Q^) "'^ has been introduced for dimensional reasons. We, now, 
perform the double Mellin transform and the Fourier transform using the fact that: 



d^q^e^'-^^Siql) = tt, (46) 
I d'q^e''-^-{ql)-'-'^ = 7^F,{e){h^f\ (47) 
j d^q^/'-'^^{ql)-^-'''\nql = -nFk{e){bY' \nb^ - 7r^^{bY% (48) 



According to this. Eg. (1441) has, after Mellin and Fourier transform, this structure: 

00 n 

#(°)(Q^6^ao,e) = ^^2(e)[(Q')"^«o]"-1( W)~^«o]'= + (50) 

n=0 k=0 

00 n 

+ In Q V 5^ ^ ^1°) (6) [(Q^)-ao]--1(l/6^)-^ao]'= 



ra=l k=l 



6.4 The resummed exponent in renormalization 
group approach 



At this point, we calculate the resummed exponent that appears in Eq.( fT9l) : 



l/fe2 /i 



T^'\Q',fi',b',ao,e) = In 



\W(^){l/b\b\ao,e) 



oo n— 1 



n=l k=0 



(oo n 
n=l k=l 



(51) 

n— fc 



x[(l/62)-ao]' 



where the last term has not been expanded because we must take into account that 
in the bare coefficient function Eq. (jH]) there is only one explicit logarithm. From the 
explicit calculation to order 0{ao), we find that 



<(.) 



2tt {AtiY 

^r(i-e) 

Att {Any 

Yrri-e 



(52) 
(53) 



Now, we want to rewrite Eq. fl5ip in a renormalized form. To do this, we use, as 
expalined in Chapter [3], the fact that (Q^)~^a;o and (l/6^)~^a;o are renormalization 
group invariant. Consequently, we may write: 



{Q'r^ao = «,(Q2)Z("^)(a,(Q2),e), 
(l/62)-ao = a.(l/62)Z("=)(a,(l/62),e) 



(54) 
(55) 



where Z^"^^(as(/i^), e) has multiple poles at e = and fj? is the renormalization scale 
which for simplicity has been chosen equal to the factorization scale. Furthermore we 
note that the universal functions Z^^\Ni, as{fi^), e)Z^^\N2, asifJ^'^), e) that extract 
the coUinear poles from the coefficient function simplify in the first line of Eq. flFIj) . 
Thus the renormalized expression of Eq. fl3Tl) has the form: 



-p/V -2j2 /-2\ \ 
1/62 /i^ /i^ 



oo n— 1 



EE^-(e)«r'^(/i^)«r(i/&^) 

m=l n=0 

(oo m 
l + lnQV5^^^1(e) 
m=l n=l 



(56) 



^2 = 1/62 



xa,(Q^)'"-"<(l/6^) . 



The resummed exponent is clearly pole-free and so we can exploit the cancellation of 
the poles that could be present in the coefficients E^^{e) and i?^„(e). Furthermore 
we want to perform a comparison with previously released resummation formulae 
given in [HI EZ]- In order to do these two things we rewrite Eq. fl^Ul) in terms of the 
renormalized physical anomalous dimension r(^, Q;s(/i^), e) and therefore, we 
calculate, according to Eq. (l20l) . the logarithmic derivative of Eq. (j56ll : 



7(l,/iV,a,(/i^),e)= (57) 

rj OO n—1 

^ ^ ' m=l n=0 

1 + in/i^&^ E^=i Er=i EUehsm^r-a^iim 

1 + ln/^^&^ E^=i EZi EUe)a.{fi'r--c.^{l/b^) 



where 



P^''\asifi^)) = -tasif) + (58) 

and [3{as{ji^)) = —Plasi^fi^) + 0{a1) is the usual four-dimensional /3-function. Now, 
in order to isolate the terms which contain the explicit logs from the rest we add and 
subtract the term 

oo m 

EE^-n(e)«.(/i^)'^-"<(l/6^). (59) 

m=l n=l 

After this, we re-expand the various terms in powers of Q!s(/i^) and agil/h'^), but not 
in powers of ln/i^6^. The result that we find in this way has the following structure: 



/-^^ dn' ( infi'b' YZ=. EZi rL(^)«r"(/^^)<(i/&^) 



To show the cancellation of divergences we rewrite the integrand separating off the 
6^-independent terms as in Chapter [3] and Chapter |H 

7(l,/iV,a,(/i2),e) = f(^)(a,(/i2),e)+f«(a,(/i2),a,(l/fe2),e) 

+f'«(ln/iV,«,(/i2),«,(l/62),e), 

where 

oo 

fW(a.(/i2),6) = Er-W«r(/i') (61) 

m=l 

oo oo 

f«(a.(/i^),a.(l/62),e) = E E r^+-(e)«r(/^')<(l/&') (62) 



m=0 n=l 



(63) 

We know that 

a,(l/6^) = /K(/x2),ln/iV), (64) 
so, in principle, it is possible to invert this relation in order to obtain: 

ln/iV = (7(as(/u'),«s(W)). (65) 
The function g has the following perturbative expression: 

1 



gia.i,% a.Hm ^ - - ) | 1 + > J /.«;(,^) | . (66) 

Substituting this expression in Eg. (1631) and re-expanding in powers of asifi"^) and 
asil/b"^), we obtain that: 

f'«(ln/i26^ a,{fi'), a,{l/b'), e) = f ^(^.(/x^), as{l/b'), e) - f «(a,(/.2), ^.(/x^), e),(67) 
where 

oo oo 

f«(a.(^^),a.(/i^),e) = f ^+„„(e)ar(/i^)<(/i^). (68) 

m=l n=l 

We choose as counterterm, 

Z(^\a,m'),e) = f«(a,(/x2),a,(/i2),e). (69) 

With this choice we obtain: 

7(l,/iV,«,(/i2),e) = f(^)(«.(/S'),e) + r^'n«s(/i'),a.(/i'),e) + 

+f«(a,(^2), «,(l/62), e) - f«(«,(^2), a,(^2), e) + 
+f«(a,(/.2),a,(l/62),6) -f«(a,(/.2),a,(/i2),e). (70) 

The first line is a power series with coefficients which are pole-free for each 6^, because 
the second and the third lines vanish when 6^ = Hence, the sum of the second 

and the third line must be finite at e = 0, but it is not necessarily analytic in ^^(/x^). 
To find its perturbative expression in powers of ^^(a^^) we rewrite the last two lines 
of Eq.(I70]) as 



l/b2 



There could be other residual cancellations of e = poles between these two terms, 
but their sum must be finite at e = and analytic in a.^ifi'^) and a<((l/6^). Thus, 



we get a perturbative expression of the physical anomalous dimension with finite 
coefficient: 



7(l,/xV,a.(/i')) = -/ '^A{a,{^^^))-B{as{p?))- 



2 



J^C{asm'),a,{f,')), (72) 



l/b2 



where 



A{a,{fi')) = 5^A„<(/z^) (73) 

n=l 

oo 

B{a,m')) = Y^B^aTifi') (74) 



m=l 

oo oo 



C(a.(/i2),a,(/i2)) = 5^5^C„„ar(/i')<(/i') (75) 

(76) 



m=l n=l 



After an integration by parts of the first term we obtain an expression for the all-orders 
resummed exponent: 



ln^A(a,(/x2)) + fi(«,(/x2)){77) 

Ji/b^ 

To obtain a LL resummation we need only a coefficient {Ai) and to obtain a 
NLL resummation we need four coefficients Ai, A2, Bi,Cii. However it has been 
demonstrated by explicit calculations [671 EH] that 

Cu = 0. (78) 

Therefore for the resummation at the NLL level we need only three coefficients 
Ai, A2, Bi. Our general formula reduces to that of Ref. [H] when 

C^n = 0. (79) 

This restriction could be a consequence of the factorization of soft emissions from the 
hard part of the coefficient function, but this remains unproven. 

The result for the anomalous dimension in Eg. (1721) can be rewritten in our for- 
malism performing the change of variable 



We get 



-^^'^ dn' 



7(l,r& ,«.(r)) = - / — G(«,(/i'),«,(/i7n')) + G(«.(/x')), (8i: 

\i n' 



where 



m=0 

OD 



12) 
i3) 



m=l 



The case of the resummation formula of Ref . [6] is obtained when Gmn is non- vanishing 
only when m = 0. In this case, we have 

' dv' 

7(1, asifi')) = - I — G{a,{p?/n')) + G{a,{-^')). (84) 



6.5 Logs of q\ vs. logs of 6^ to all logarithmic or- 
ders 



Large logarithms of q±_ appear in the perturbative coefficients in the form of plus 
distributions. We define 



InP(gi)^ 



in such a way that 



InP(gi) 



Let us consider the Fourier transforms 



lnP(gj 



2vr / q±dq_LJo{q_ib) 
Jo 



\nP{qj 



where we have used the definition of the 0-order Bessel function Jq- 



2tt 



iz cos d 



Jo{z) = — dee 

Jo 

We now exploit the definition of the plus distribution: 

Ip{Q%') = 2n / du[Mq±b) - l]—^ + 27r / dq^Mq±b)- ^ 



(1± 



^5) 



(86) 



^7) 



^9) 



Writing In^g^ as the p^^ a-derivative of at a = 0, we get 



2n 
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TT- 



QP 

daP 

QP 



2a- 1 



dq±q^ 



2a- 1 



daP 



0%' 
4 



T(a) 



rfl-a) 



1 

a 



(90) 



where the last equahty follows from the identity 

.ir(l/2 + z//2 + /i/2) 



dxx^Jy{ax) = 2^ a ^ 



r(l/2 + p/2 - /i/2) 
a > 0; -Rez/ - 1 < Reji < 1/2. 



From Eq.( l90i) . we read off the generating function G{a) of Ip 



(91) 
(92) 



G{a) 



TT 



a 



Til -a) 



(93) 



in the sense that 



2i,2^ 



dP 
daP 



G{a) 



(94) 



a=0 



Now, the generating function of logarithms of Q^lP' is {Q'^lP') " in the sense that 

' dP 



Lp ^ ln^(l/(Q^6^)) 



daP 



(95) 



a=0 



Inverting Eq. (!93|l , we find the relation between the generating function of Lp and the 
generating function of Ip, which is 



—S{a)[aG{a) + vr], 

TT 



where 



S{a) 



1 r(l-a) 
4^r(l + «)■ 



(96) 



(97) 



Performing the Taylor expansion of the r.h.s. of Eq. (1961) around a = and using 
Eq.(l9ll), we obtain: 



^Q%-)-^ = - E ^ E (""l^^'-l^^^^^HO) 
TT ^-^ ml ^-^ \ I I 



m=0 i=0 



(98) 



where S^^\0) is the j-th derivative of S{a) evaluated at a = 0. Now, using Eg. (195 



we the relation between Lp and Jpi 



(99) 



Thanks to the first equality in Eg. (1871) and to the fact that 

1 d'' 



InP-'' ql 



p{p-l)---{p-k + l) d\ir ql 



\ifqi, 



(100) 



we arrive at a relation to all logarithmic orders between the logs of 6^ and the logs of 



i^ g(^-^)(o) /• Ax^.,-,., 



:[{k-i)\ J Q 



d^ 



qi dhi^ q]_ 



lYi^ql 



(101) 



This relation allows us to derive the relation between a generic function of \n{l/Q%'^) 
and a function of Ing^. Indeed, given a function 



/ 1 \ °° 1 



(102) 



Eg. (|TOTll implies: 
1 



h In 



1 - ^(fe-i)(0) rd\^ 



TT ^ 



d^ 



qidhi^ql 



■h{\nq]_] 



(103) 



The r.h.s. of Eg. dlOSp can be viewed as the Fourier transform of a function (more 
properly a distribution) h{\nqY): 



h In 



Q262 
hihi q\ 



d'q± 



e"^^-^h{\nql), 



1 - 5(fc-i)(0) 

71 {k-iy. 



d^ 



qi dhi'^ql 



Mlnql] 



(104) 



6.6 Resummation in g^-space 

In this section, we investigate the conseguences of our general result Eg. (11041) for 
the resummation at the NLL level of logarithmic accuracy. According to eg (ll9l) and 
the discussion below and according to Eg.( !77|) . we have that our resummation factor 
formula in Fourier space is: 



K'''iQ%l l/h\ Q^) = exp {E''\Q%\ l/6^ Q^) } , (105) 



where 

E-{Q'bM/b^Q') = r ^T-{%r^'b\a,m')), (106) 
and where at NLL level 

r^lL(C,/^V,a.(/z2)) = -ln%[A^a.{^^') + A2a',{^')]-B^asif^') 
fx /i 

-— a«(/i)ln — — , (107) 

Po asiix-^) 

where we have used the definition of the /5-function: 

/i'^a.(/u2) = Pas = -l3oal - P^al + 0(af) (108) 



and where we have used the change of variable 

(i/x^ dag 

ix^ (3{as) 



(109) 



to compute the integral that appears in the last term of Eg. 0771) . Now, thanks to 
Eq. fll04p . we can rewrite the resummed exponent in 6-space Eq. fll06p in terms of a 
resummed exponent defined in g_|_-space. Thus, up to NNLL terms, we obtain: 



r-^(gi,gi,g^) 



;iio) 



where 



-a,(g^)ln ^ (111) 



(3o asiqj 

and where the relation of the constant coefficients of this last equation and the of 
Eq. ffTOTl) is 

ii = ^ (112) 
A. ^ -7^ + ^/3oln^) (113) 

TT vr 4 
= (115) 

TT 

Here je is the usual Euler gamma. Now, we want to define a resummation factor 
in g^-space. Looking at Eq.( JT8i) . we note that large InQ^fo^ of coUinear nature are 



resummed by the parton distribution function. So, in order to define a resummation in 
(jx-space, we must take them into account. For simphcity, we consider the resummed 
part of non-singlet cross section, because the non-singlet parton distribution functions, 
which are defined as 



(AT, = /„(iV, ^^^) - f,(N, /i^) a,b^g, 



(116) 



evolve independently. In particular, in Mellin moments they satisfy the following 
evolution equations: 



d 



;ii7) 



Hence, the non-singlet parton distribution functions evaluated at /x^ = 1/6^ are related 
to the ones evaluated at /i^ = by. 



/:,(iV,l/62)=exp 



1/62 II 



118) 



Thus, the resummed part of the cross section with the non-singlet parton distribution 
functions evaluated at /x^ = becomes 



exp 



(119) 



exp 



l/fe2 /^^ 



r-(§,/iV,a.(/.2)) -X^7'(iV„«.(/^')) 



The general relation between a function of InQ^^^ and its Fourier anti-transform 
Eq. fll04p . immediately enables us to define a resummed exponent of the non-singlet 
part of the cross section in g^-space, which is: 



T^res(;-2 2 1 \- 5'^^ (0) 



X 



exp 



1 

g(^-^)(o) 

r {k-iy- 



r,2 



r,2 ' 



11 



/i 



X 



d d 



k-l 



exp 



djj 



r,2 



dqldln'^-^ql 

S^''\0) d^ d 



,0^ 



/iVgi,a,(/i2)) - ^y(iV,-,a,(/i2)) 



lim — 



X 



9{ql -r/)exp 



2 /i2 



r^^^(^,/iV?i,«s(/i')) - 5^7'(iv„«s(/x')) 



(120) 



where the last equation defines imphcity the g_L-space resummation exponent: 



^^^^(gi, gl, Q') = exp {E^^'iql, ql, Q')} 



(121) 



All the previously released expressions for this exponent given in [69l ITOl 171] are par- 
ticular cases of this general expression. They differ essentially in the criteria according 
to which the sub leading terms are kept. 

We want to calculate the NLL result in gj_-space. Thus, keeping only the terms 
up to NNLL in Eq.([l20]) we obtain 



K^LLiql ql Q') = -i\ 0{qi - r/) exp 
TT qj_ 



'^A''^ f -rres /Q'^ -2 / 2 /-2\\ 



j—l / J j,_n ■ J 



k=0 



In order to compare this result at NLL to that of [70j, we define a new variable h: 

h = 2\nqlAias{ql). (123) 



In terms of this variable and using Eq. fl971) the series that appears in Eq. fll22p can be 
computed: 



fc=0 



k\ 



[-\MlA,aMl)T = S{-h/2) = 2^ ^|| + |^/^| . 



(124) 



In conclusion, we obtain that the NLL resummation factor becomes 



K^NLiql ql Ql = -ii 0{ql - ^) exp 



'^A'"^ f T^res /Q^ -2 / 2 / -2\\ 



2 /i^ Y /i^ 



5^7'(iV„«,(/i2)) 



r(i + V2) ) 
r(i-V2)J' 



(125) 



which gives the same result given in Ref. [70] in the case that the coefficient Cu that 
appears in Eq. (11071) is equal to zero and that the arbitrary constants Ci and C2 also 
defined in [TOj are equal to one. It is clear that the last two terms of the exponential 
of our result let the non-singlet parton distribution densities, which enter the q]_ 
derivative, evolve from the scale to the scale q\. 

We conclude the chapter noting that also in this case the resummed results using 
the renormalization group approach are less predictive than results obtained with the 
approach of Ref.fT], as it is shown in Ref. [6]. Furthermore the conditions that reduce 
our results to those of Ref. |6]| in terms of factorization properties is still an interesting 
open question. 



Chapter 7 



Predictive power of the 
resummation formulae 



We have shown that the renormahzation group resummed expressions are less predic- 
tive than those obtained with other approaches discussed in Sec J2.3[ In this Chapter 
we shall compare the various approaches quantitatively. We will show that all the re- 
summation coefficients can be determined by a fixed order computation. This can be 
useful, because the determination of the resummation coefficients from a fixed order 
computation represents a possible way to check the correctness of the resummation 
formulae with strong factorization properties. 

In particular, in this chapter, we will show how the resummation coefficients Qmnp 
(for the prompt photon case), gmn (for the DIS and DY cases) and Gmn,Gm (for 
the DY transverse momentum distribution) can be determined. For the rapidity 
distributions of DY and DIS they are the same of the all-inclusive cases (see Chapter 
|5]) The resummation coefficients are determined by comparing the expansion of the 
resummed anomalous dimension 7 in powers of as(Q^) with a fixed-order calculation, 
which in general has the form: 



where jpolNjas) is the physical anomalous dimension for each individual partonic 
subprocess for the promt photon case, for the qq channel in the DY case and for 
the q channel in the DIS case. For the case of the small transverse momentum DY 
distribution it ha the form: 



The number fcmin is the minimum order at which the anomalous dimension must be 
calculated in order to determine its N^~^LL resummation. 

For prompt photon production, the number of coefficients Nk that must be deter- 
mined at each logarithmic order, and the minimum fixed order which is necessary in 




(1) 



(2) 
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Prompt photon 




Eq.flMl) seclOl 


Eq.fl96l) secJO 


Eq.flHTD secIO 




2k 


fc(A:+3) 
2 


fc(fc+l)(fc+5) 
6 


h ■ 


k + l 


2k 


3A; - 1 



Table 7.1: Number of coefficients iV^, and minimum order of the required perturbative 
calculation fcmin for inclusive prompt photon N^~^LL resummation. 





DIS/DY 




Eq.rtHT]) sec JO 


Eq.dHUD sec JO 




k 


fc(fc+i) 

2 


h ■ 


k 


2k - 1 



Table 7.2: Number of coefficients and minimum order of the required perturbative 
calculation fcmin for inclusive DIS and DY N^~^LL resummation. 





DY transverse distribution 




Eq.(l84D secJf).4l 


Eq.flHll) secJf).4l 




2k - 1 


fc^+3A;-2 
2 


h ■ 


k 


2A; - 1 



Table 7.3: Number of coefficients A^^^ and minimum order of the required perturbative 
calculation k^m for small transverse momentum DY N^~^LL resummation. 



order to determine them are summarized in Table 17.11 according to whether the coeffi- 
cient function is fully factorized [Eq.( j9^ sec 14. 41 ]. or has factorized A^-dependent and 
iV-independent terms [Eq.(l96l) sec J4.4j . or not factorized at all [Eq. (IHTll se d4.4j . In the 
approach of Refs. ITi] the coefficient function is fully factorized, and furthermore 
some resummation coefficients are related to universal coefficients of Altarelli-Parisi 
splitting functions, so that femin = k. For prompt-photon production, available re- 
sults do not allow to test factorization, and test relation of resummation coefficients 
to Altarelli-Parisi coefficients only to lowest 0{as)- 

The results for DIS and Drell-Yan, according to whether the coefficient function 
has factorized A^-dependent and A^-independent terms as in Refs. [HH El [1] [Eq. ( JSTl) 
sec J3.3j or no factorization properties as in [(Eq JHOil sec J3.3j . are reported in table 17.21 
Current fixed-order results support factorization for Drell-Yan and DIS only to the 
lowest nontrivial order 0{a1). 

In table 17.31 we report also the results for the small transverse DY resummation. 
We list Nk and /cmin for the approach of Ref. [6] [Eq.lJH^ sec i6.4j and for the renormal- 
ization group approach [Eq. (IHTI) sec J6.4j . If the two cases are related by factorization 
properties of the cross section is not yet understood even if probable. 

In the following, we present all the proofs of these results. 



7.1 Prompt photon production in the strongest 
factorization case 

This is the case of Eq.f l94p in section l4?4l In this case there are Nk = 2k non- vanishing 
coefficients gooi and goio i = 1,2, ... ,k. The resummed expression of the anomalous 
dimension at N^~^LL is given by: 



-N T 

7(iV,a.(Q^))= / ^[Y^g,,,a\{Q-'/n' 



n 



^ [t.9oo.al{Qyn)] . (3) 



We consider ffist the second integral in Eq.(l3]). Noting that: 

dn das{Q'^/n) 



n 



where 



6n = 



/3o 



bi 



47r (47r)2' 

with (3q and (3i given in Eq. flHHj) in section II. 2^ we can rewrite it in the form: 



da^ 



Ek i 
i=i QooiO^s 



(4) 

(5) 
(6) 

(7) 



Now, we expand up to order a^"^ each term that compares in the integrand of this 
last expression and collect all the coefficients that correspond to the same power of 
as- Doing this, we have that the integral ([7]) can be rewritten in the following form: 



da^ 



I Ja,{Q'^) 



9ooi 
a. 



+ {blgooi + 9002) + (blgooi + blgoo2 + fi-oos)". H h 



+ (&fe_ifi'ooi H \- bl_\gook^i + gook)a 



k-2 



where k > 1 and the numbers 6^ are build up with the coefficients of the (3 function. 
Now, we perform the integral over a^. We get: 

k 



ij^^ooi In (^^^^) + ib\gooi + goo2)[as{Q' /N) - a^iQ^)] + 

+ ^(^'2^001 + &2^002 + ^003)[«'(QV^) ' + " " " + 



+^(&Li^7ooi + ■ ■ ■ + 6^:1^700.-1 + gook)[a^s-\QVN) - «^'(g')] ) (9) 



To perform the first integral of Eq. ([3]) , we it is sufficient to note tfiat in this case 



and proceed as before. The result that we obtain is: 



(10) 



dn 
n 



^goioaliQ^/^ 

i=l 

= 1^010 In ( "^^^(gl^'^ ) + (^1^010 + 9o2o)[as{QVN') - a.(g^)] + 

+^(^2^7010 + blgo2o + go3o)[al{QVN^) - a'(Q')] + ■ ■ ■ + 

+^(&Li^?oio + ■ ■ ■ + btlgok^io + gokoM'^Q'/N') - a^l(Q2)]| (11) 

At this point, we take the first term of Eq.(l9]) together with the first term of Eq.l fTTl) 
in order to isolate the first contributions of Eq.(l3]). We have: 



1 

2^0 



^010 In , . + 2^001 In , , 



(12) 



From this contribution, we want to extract the firs two LL terms. Hence, using the 
one loop running of ag^Q"^ /N°-),a = 1,2 (EqJMl of section [TT^ . we obtain for this 
contribution: 

a,{Q^)ln^[-{gooi+gow)] + a^(g') In^ l[6o/2((7ooi + 2(7oio)] 

+ 0{al\n{N)) + 0{al+'\n\N)), (13) 

where i > 0,1 < j < i + 3. Now, we take the second term of Eq.Q together with 
the second term of Eq. ffTT]) in order to keep the second the second contributions of 
Eq.® and we have: 

7^ [{blgoio + 9o2o){as{QVN') - as{Q')) + 2{blgooi + goo2){as{QVN) - a.(Q'))] . 

ZOq 

(14) 

From this contribution, we want to extract the first two NLL terms. In order to 
do this, we use the one loop running of as{Q'^/N°'),a = 1,2 and observe that the 
coefficients of (7001 and goio are modified by the last two terms of Eq. (JT3l) . We get: 



ln(l/A^)[-(c}5(oio + dlgooi + 5-020 + 5-002)] + 
a^iQ^) \n\l/N)[bo{2clgoio + dlgooi + 25020 + 5002)] 
+ Oialln{N)) + C>(af+Mn^(Ar)), (15) 

where c]^, c/j^, cj, d\ are coefficients (that are of no concern to us) and ^ > 0, 1 < j < 
i + 3. This procedure can be repeated for all the other contributions. So, we take the 



k-th term of Eq.(l9]) together with the k-th term of Eq.l JTT]) in order to keep the k-th 
contributions of Eq.(l3]). For the general k-th contribution, we get: 



1 

2b'o 



k-1 



(16) 



From this contribution, we want to extract the first two N^~^LL terms. In order to 
do this, again, we use the one loop running of UsiQ"^ /N""), a = 1, 2 and note that the 
coefficients 5^001 and g^iQ with z = l,2,...,/c — 1 are modified by the previous terms. 
For this generic term, we get: 



a] In — [-(4_i5(oio H h c^_}fi'ofc-io + 4-i5'ooi H h dl_\gook-i + 5'ofco + 5'oofc)] + 

In^ i- [60 A;/2(24_ 1^010 + ■ ■ ■ + 2c^:}^ofc-io + 4-i^ooi + • • • + dt\9mk-i 
+2^ofco + ^?oofc)] + In(iV)) + ln^(iV)), (17) 



where i > 0, 1 < j < i + 3. 

To summarize, Eqs. fll3lll5|17p tell us that from the expression of the physical 
anomalous dimension Eq. ([3]) , we can extract the following linear combinations of the 
coefficients ^foio, ■ • ■ , goko, fl'ooi, • • • , Qook- 



h — —{gooi + fi'oio) 

h = -^{9001 + '^goio) 

h = -{c\goio + d\gooi + go2o + goo2) 

h = bo{2clgoio + d\gooi + 2go2o + goo2) 

1 /I k— 1 jl jk— 1 \ 

hk-l = —[Ck-l90W + ■ ■ ■ + C^„i5'0fc-10 + "fc-lfi'OOl + \- "fc_i5'00A:-l + fi'OfcO + 900k) 

hk = ~;^(2Cfe_i(7oio + ■ ■ ■ + 2c^_]^(7ofc-io + '^fc-i5'ooi + ■ ■ ■ + c?fc_i5'oofc-i + 2(7ofco + fi'oofc)) 



with k > 2 and k, i = 1,. . . ,2k the known terms. These are 2k independent lin- 
ear combinations that determine the 2k coefficients (7010, • • • ,9oko,9ooi, ■ ■ ■ ,9ook of a 
N^~^LL resummation comparing them with the correspondent terms of Eq.(l2]) up to 
order aj"*"^. This is a direct consequence of the fact that the two vectors (1, 1) and 
(1,2) are independent. This shows that, in order to obtain a N^~^LL resummation 
in the case of the strongest factorization (Eq. fl94l) in section , we need to know a 
N^^^LO fixed order calculation of the physical anomalous dimension. Hence, in this 



7.2 Prompt photon production in the weaker fac- 
torization case 



This is the case of Eg. fl96p in section I4.4[ In this case in order to perform a LL 
resummation, we need two coefficients {goio, gooi)', to perform a NLL resummation 
three more coefficients are added ((70205 5^002, fi'oii); in general to perform a N^~^LL 
resummation k + 1 coefficients are added {g^ij + j = k) to those of the N^~^LL 
resummation. Hence, in order to perform a N^~^LL resummation, we need to deter- 
mine 

iV. = E(p + l) = ^^, (18) 

coefficients. We want to determine fcmin in Eq.([2]) so that all the k{k + 3)/2 are 
fixed by the same number of independent conditions obtained from the fixed order 
expansion of the resummed physical anomalous dimension. We note, first of all, that 
this happens if we can extract 2 independent conditions from the LL contributions, 
3 from the NLL one and k + 1 from the N^~^LL. The N^~^LL expression of the 
physical anomalous dimension in this case is given by: 



/N 1 k s-1 
^ E E 90^s^^<iQV^')<"iQ'M■ (19) 
s=2 1=1 



The first two LL contributions have been already computed and are given in Eq. f[T3|) . 

Now, we extract the first 3 NLL contributions of Eg. (1191) . Recalling the deriva- 
tion of Eq.([T3j), we obtain that these contributions are contained in the following 
expression: 

7^ [{b\goio + go2o)MQVN') - «.(g')) + 2(61(7001 + goo2)MQVN) - a,{Q'))] 
+ / —goiias{Q^/n^)ois{Q^/n). (20) 

Since 

c^.(QVp")= , , T^^l,, , +Q(af-ln^l), ^>0, l<J<^ 



l + a/5oa.(Q2)lni ^ ^ p 



1 



p 

00 



r = E(-)'«'^o«i(Q')in^'^, (21) 



l + a6oa,(Q2)ini ^ ^ ^ p 



and keeping in mind that corrections to the NLL come from Eq.( fT3l) . we have that 



Eg. (1201) become: 



oo 

— [2^(4^.(7010 + gm,) + 2(4^?ooi + ^?oo2)]6'o«^'(Q') In^' ]^ + (22) 



The Cauchy product of the two series in Eq. (l22l) is given by 

oo ^ CO ^ 

{Y.^-)%ai{Q') In^- -)E(-)^2^6X(Q') In^ -) = 



i=0 



and 



Now, because 



= 2(2^' - 1). 



i=l 



n n 



(23) 



(24) 



(25) 



we can perform the integration. We get 

3 



2(2^ — 1) 

2\cljgoio + go2o) + 2(^^5-001 + ^'002) H ^ 5-011 



n 



(26) 



where Cy\ are certain coefficients we don not need to worry about. The first three 
NLL contributions are given by j = 1, 2, 3. 

The last step is to extract the first k + 1 N^~^LL contributions that come from 
Eq. ffT^ . Recalhng how Eq.( ffTBl) ) was computed, we have that the desired contribu- 
tions are contained in the following expression: 



1 
^0 



^(&Li5oio + ■ ■ ■ + &t}5ofc-io + 9,M-\QVN^) - «s"'(Q')) 



+ 



+ ;^(&Li5ooi + ■ ■ ■ + hl-_\9o,k-i + 9mk){a]-\Q'lN) - a'-\Q')) 



+ 



(27) 



We use the following relations 



1, 



<(Q'/P'') = A fn2^^ iv + 0(«f Mn^ -), ^ > 0, 1 < j < ^ + 1, 



(l + a6o«,(g2)inl 



P' m=0 



r + m — 1 
m 



a'"6™a™(g2)ln'" 



1 

P 



{2i 



where (^) are the usual binomial coefficients. With this, we can compute the integral 
in Eg. (1271) performing the Cauchy product of the two series expansion of as{Q'^/n'^) 
and of as{Q^/n) and performing explicitly the integral using the change of variable 
Eq.®. We get: 

/jV , k—l k—1 oo ^ 

i=l i=l m=0 

(29) 

where 
and where 

n \ r(n + 1) 



-ij r(o)r(n + 2) ^^^^ 

is equal to 1 for n = —1 and otherwise. Thefore, keeping in mind the calculation 
of Eq. fll6l) and taking the first k + 1 N^^^LL contributions of Eq. fl27p . we have 

C (\m+l /y^ _|_ ^ _ ]^\ 

2^ V^^l )[2™(4-i™^7oio + --- + c^:L^ofc-io + ^7ofco) + 

m=0 V "T / 



fc-1 t-l 

+ (4-lmfi'001 H \- d\z\jgOOk-l + fi'OOfc) + XI XI 90U-ifum + 

t=2 i=l 

+ '^9o^k-^Ci^''~'^\b^a':^-{Q') ln-+^ 1 (32) 

i=l ^ 

where c^_im; '^fc-im; fum^^ certain coefficients we do not have to worry about. At 
this point, we can make some simplifications. In fact, since 

rio,k) _t}r^(^ + k~l\ _ (-ir+i /m + A:-l\ 

^"^ " m + 1 V ^-1 )~ k-1 V rn + 1 J ^^^^ 

and (see Appendix [E]) 

(^(fc,o) _ 2"C(,°'^\ (34) 
we can write Eq. (l32l) in the following form: 

X I [^m'°^('^fc-lm.fi'010 -I ^ C^_lrnfi'Ofe-lo) + 

m=0 

k-1 t-l 

+C^'''\dl_^j^gOOl H h C?fc_ijfi'OOfc-l) + C'm'''^ ^^^^OOit-ifum + 

+ i29o^k-^C^::'-Ab^a':+-{Q') ln-+i 1, (35) 

i— n J 



What this result tells us, is that passing from the N^~'^LL to the N^~^LL resum- 
mation, k + 1 new resummation coefficients are added. In Eq.( |35ll . we have + 1 
conditions for this coefficients (one for each m) to be set equal to the correspond- 
ing fixed order contribution of Eq.([2]) • We shall now show that this conditions are 
independent. This is equivalent to showing that the k + 1 linear combinations 

k 

Y,9^^k~^C';i;^-'^ m = 0,l,...,fc (36) 

j=0 

with 

Ci^^) ^ + = 2' ^) ("^ -]^_[ - (37) 

1=0 ^ / \ J / 

are independent. Moreover, this is equivalent to showing that for each k the columns 
of the (A; + 1) X (A; + 1) matrix A^- = Cm^~^^ are independent vectors. To show this, 
we need to use two identities proved in Appendix [El 

^ 2"'C^y (38) 
C^^^ = 2Cl^^-''> - Cf'^h z,j>l, (39) 

(k) 

which allows us to to compute the columns of the matrix explicitly for all k. To 
show their independence, we use induction on k, i.e. we demonstrate the independence 
of the columns for k = 1 and then we assume that the property is valid for A; — 1 to 
prove that it remains valid for k. In the case A; = 1, we have a 2 x 2 matrix: 

(C'r^ C'£'°^) = (l 2-), (41) 



and now it is clear the two columns are independent, because 1 and 2"^ are independent 
functions of m. Now, using the induction hypothesis 

fe-i 

J2 a.C'(;''=-^-*) = ^ a, = 0, ^ = 0, . . . , A: - 1, (42) 

i=0 



we want to show that it is sufficient to prove that: 



k 

J] /^.C'^i''^-^) = ^ A = 0, z = 0,...,fc. 

1=0 



(43) 



Using the relations ( l38|39l) . we have: 

k k-l 

Y,^M''-'^ = (/3o + 2-/?fc)Ci°'^) + J]/5,C'^''=-^-) 

i=0 3=1 

k-l k-l 

j=i j=i 

fc-l k-2 

j=l j'=0 
k-2 

+2/?,_iC'i^-i'0) = (44) 

Now, from Eg. (1401) . we know that Cm''^ is a degree-(A; — 1) polynomial in m. Fur- 
thermore, from Eq. (!39ll . we know that the vectors Cm^ ^ with j = 0, . . . ,k — l are 



at most polynomials of degree k — 2 in m. Consequently, Eq.(jlll) can be satisfied if 
and only if 

(/?o + 2"^(3,)C^^''^ = 0, (45) 

k-2 

-PiCi^''-'^ + - Pj+i)Cl^'-'~'^ + 2/?.-iC'i^-^'°) = 0. (46) 

From the first, it follows that: 

Po = Pk = 0, (47) 
while from the second, thanks to the induction hypothesis Eq.f l42p . we have that 



A = Pk-i = 0, (48) 

and that 

Pi = Ip2 = \p3 = --- = l^^Pk-1. (49) 

In conclusion from Eqs. fl47ll48l|^9l) it follows that 

/3, = 0, j=0,...,fc. (50) 

This completes the proof that the columns of the squared (/c + 1) x (A; + 1) matrices 
A^^- = Cm^~''^ are independent for all k. This shows that, in order to obtain a N^^^LL 
resummation in the case of the weaker factorization (Eq. (!96l) in section I^Til) . we need 
to know a N'^^LO fixed order calculation of the physical anomalous dimension. Hence, 
in this case /cmin = 2k. 



7.3 Prompt photon in the general case 



Let us now consider the most general case, in which the coefficient function does 
not satisfy any factorization property. This is the case of Eq. fl87p in section I4.4[ In 
this case, in order to perform a LL resummation we need 2 coefficients (fifooi? fi'oio); 
to perform a NLL resummation 5 coefficients are added ((7002, S'loi, 5^020? fi'iio? fi'oii) and 
to perform a N^~^LL resummation k{k + 3)/2 coefficients are added {gmnp with 
m + n + p = k but without n = p = 0). Thus, in order to perform a N^~^LL 
resummation, we need to determine 

^ p(P + 3) fc(fc+l)(fc + 5) ^^^^ 
p=i ^ ^ 

coefficients. The N^-^LL expression of the physical anomalous dimension is given, 
in this case, by 

J k s—l s—i 

7(iV,a.(g^))= / -J2T.T.9^^^-^-^<iQ>iiQV^")^r''iQVn). (52) 

/I ^ 

""J^ s=l i=0 j=0 

Now, we proceed in the same way as we have done in section [72] and we find that the 
k{k + 3)/2 new coefficients that are added passing from the N^~'^LL to the N^~^LL 
contributions appear only in the following combinations: 



fe— 1 k—i 00 ^ 

E E 9^^>^-^-^ E b-a.m'^"' In"-^' ^. (53) 

i=0 j=0 m=0 

Each term with fixed m in the expansion Eq. (!66l) provides a new condition on these 
coefficients. However, these conditions are not linearly independent for all choices of 
m. Indeed, let us define the matrix Cm'' ' "'^ = D^'^}. where the lines are labelled 
by the index m and the columns by the multi-index This matrix gives the linear 

combination of the coefficients gmnp in Eq.(|66ll to be determined and it turns out to 
be of rank 

rj«',.,) = 2k < M^±^. (54) 



D^'^}. is a M X matrix , whose columns are the M-component vectors 



We shall now prove this statement: 
, is a M X ^ ■ 

DOJ) = C^.'^-^-J); 0<i<k-l; 0<j<k-i; < m < M. (55) 
We use induction on k. For k = 1, D^^^ is a 2 x 2 matrix with columns 

/I \m+l 

^« = W'^\ = (1, 2-), (56) 

that are linearly independent; the rank of D^^^ is 2. Let us check explicitly also the 
case k = 2. In this case 

n(2) _ (fi{0,l) M1,0) M0,2) ^(1,1) (-y(2,0)\ 

^m — 5 '^m ' '^m ' '^m ' '^m ) ■ W J 



The first two columns are the same as in the case k = 1: they span a 2-dimensional 
subspace. The last three columns are independent as a consequence of Eg. (1431) with 
k = 1. Furthermore, Cm'^"* and Cm^^ = 2'"Cm'^'* are independent of all other columns, 
because they are the only ones that are proportional to a degree-1 polynomial in m. 
Finally, Ct^^ is a linear combination of the first two columns, as a consequence of 
Eqs. flHTllHUD with i = j = 1. Thus, the rank of D^^) is 2 + 2 = 4. 

We now assume that Z}^^"^) has rank 2 (A; — 1), and we write the columns of D^''^ 

as 

Dl^^ = {C^ri^'-'-'-^\ C^^:'-'^) (58) 
0<i<k-2, < j < k -1 -i 0<l <k. (59) 

By the induction hypothesis, only 2{k—l) of the columns Cm''^^^ * are independent. 
The columns Cm'^~^^ are all independent as a consequence of Eq.(|in]); among them, 
those with 1 < / < A; — 1 can be expressed as linear combinations of Cm''^"^ * "'^ by 
Eq. Only Cm^^ and Cm'°'* are independent of all other columns because they are 
proportional to a degree- (fc — l) polynomial in m, while all others are at most of degree 
(fc — 2). Hence, only two independent vectors are added to the 2{k — l)-dimensional 
subspace spanned by Cm'' ^ * and the rank of D'^''^ is 

2(A; - 1) + 2 = 2A;. (60) 

It follows that each individual terms in the sum over m in Eq. ( l66i) depends only on 
2k independent linear combinations of the coefficients Qijk-i-j, 0,<i<k — 1, 0< 
j < k — i. 

This means that the N'^~^LL order resummed result depends only on 2k indepen- 
dent linear combinations of the k{k+3)/2 new coefficients that are added passing from 
the N^~'^LL resummation to the N^~^LL one and that the remaining coefficients are 
arbitrary. Because a term with fixed m in Eq. fl66|) is of order aj"'"™', this implies that a 
computation of the anomalous dimension up to fixed order fcmin = 3A; — 1 is sufficient 
for the N'^^^LL resummation, because m = 0, 1, . . . , — 1. Note that when going 
from N'^~^LL to N'^LL, at this higher order, in general some new linear combinations 
of the k(k + 3)/2 coefficients, that we have added from the N^~'^LL to the N'^'^LL, 
will appear through terms depending on Hence, some of the combinations of co- 
efficients that were left undetermined in the N'^~^LL resummation will now become 
determined. However, this does not affect the value k^i^ of the fixed-order accuracy 
needed to push the resummed accuracy at one extra order. In conclusion, even in 
the absence of any factorization, despite the fact that now the number of coefficients 
which must be determined grows cubically according to Eq. ( ISTl) . the required order 
in as of the computation which determines them grows only linearly. 



7.4 DIS and DY in all cases 



This is the case of Eq. ( |80ll in section 13.31 Here we discuss the case without assuming 
any factorization property, because the factorized case will be recovered as a particular 



case. So, in the general case, at LL we need to determine 1 coefficient [goi)', at 
NLL 2 coefficients are added {go2,gii) and at N^~^LL, k coefficients {gik-i with 
i = 0, 1, . . . , /c — 1) are added to the N^~'^LL ones . Thus, in order to perform a 
N^~^LL resummation, we need to determine 

p=i 

coefficients. The N^~^LL expression of the physical anomalous dimension is given by 
j{N,a,{Q')) = a ^ ^ ^ ^_,a:(Q2)ar (Q'K), (62) 

-'1 ^ s=l i=0 

where a = 1 for DIS and a = 2 for DY. Now, we proceed in the same way as we have 
done in the previous sections and we find that the k new coefficients that are added 
passing from the N^~'^LL to the N^~^LL contributions appear only in the following 
combinations: 

fc— 1 oo ^ 

Y^gik-i <^"'^'Cl^''-'^ h^a,{Q^f^^ In'^+i -. (63) 

j=0 m=0 

Again, each term with fixed m in this expansion provides a new condition on these 
coefficients. These conditions are all linearly independent for any choice of m, because 
the k^k matrix = C^'^"* with m, i = 0, 1, . . . , fc — 1 has all independent columns. 
This is a direct consequence of the fact that each column is a polynomial with differ- 
ent degree in m. This implies that, in order to determine a N^~^LL resummation, 
a computation of the physical anomalous dimension up to order /cmm = 2/c — 1 is 
sufficient. In the more restrictive case Eg. 0811) in Section the only non- vanishing 
coefficients are go^k- This means that going from the N^~'^LL to the N^~^LL we need 
to take only one combination of the expansion Eq.( l63l) . Hence, in this more restrictive 
case, Nk = k and fcmin = k. 



7.5 DY small transverse momentum distribution 

This is the case of Eq. (18 ip in section 16.41 obtained with the renormalization group 
approach. The result of the approach of Ref . j6] (reported in Eq. ( IHT]) of section 16.41) 
will be recovered as a particular case. In the most general case, at LL we need to 
determine 1 coefficient (Goi); at NLL 3 more coefficients are added (Gn, 6*02, Gi) and 
at N''~^LL, k + 1 coefficients {Gik-i with i = 0, 1, . . . , fc — 1 and Gk-i) are added. 
Therefore, in order to perform a N'^~^LL resummation, we need to determine 

E. . k ~\~ 3/l 2 . . 

{p+l)-l = ^ (64) 



coefficients. The A^'^ ^LL expression of the physical anomalous dimension is given by 



V E E G.s-.aim')armVn') + Y,G^aim'). (65) 

—1 i=0 1=1 



Now, we proceed as before and we find that the k new coefficients that are added 
passing from the N^~'^LL to the N'^'^LL contributions appear only in the following 
combinations: 



fe— 1 oo / 1 \ 

1=0 m=0 ^ 



in')- (66) 



As before, each term with fixed m in this expansion provides a new independent 
condition on this coefficients. This implies that, in order to determine a N^~^LL 
resummation, a computation of the physical anomalous dimension up to order /cmin = 
2/c — 1 is sufficient. In the more restrictive case of Eg. (15^ sec l6.4p . the only non- 
vanishing coefficients are Go,fc and Gk-i- This means that going from the N^~'^LL 
to the N^~^LL only 2 coefficients are added but the LL where we have only one 
coefficient. Hence, in this more restrictive case, Nk = 2k — 1 and /cmin = k. 



Chapter 8 
Conclusions 



In this thesis, we have studied the renormahzation group approach to resummation for 
all inclusive deep-inelastic and Drell-Yan processes. The advantage of this approach 
is that it does not rely on factorization of the physical cross section, and in fact it 
simply follows from general kinematic properties of the phase space. Then we have 
analyzed some of its generalizations. 

In particular, we have presented a generalization to prompt photon production 
of the approach to Sudakov resummation which has been described in Chapter [3] 
for deep-inelastic scattering and Drell-Yan production. It is interesting to see that 
also with the more intricate two-scale kinematics that characterizes prompt photon 
production in the soft limit, it remains true that resummation simply follows from 
general kinematic properties of the phase space. Also, this approach does not re- 
quire a separate treatment of individual colour structures when more than one colour 
structure contributes to fixed order results. 

The resummation formulae obtained here turn out to be less predictive than pre- 
vious results: a higher fixed-order computation is required in order to determine the 
resummed result. This depends on the fact that here neither specific factorization 
properties of the cross section in the soft limit is assumed, nor that soft emission 
satisfies eikonal-like relations which allow one to determine some of the resummation 
coefficients in terms of universal properties of collinear radiation. Currently, fixed- 
order results are only available up to for prompt photon production. An order 

computation is required to check nontrivial properties of the structure of resum- 
mation: for example, factorization, whose effects only appear at the next-to-leading 
log level, can only be tested at 0{al). The greater flexibility of the approach pre- 
sented here would turn out to be necessary if the prediction obtained using the more 
restrictive resummation were to fail at order a^. 

We have also proved a resumation formula for the Drell-Yan rapidity distributions 
to all logarithmic accuracy and valid for all values of rapidity. Isolating a universal 
dimensionless coefficient function, which is exactly that ones of the Drell-Yan rapidity- 
integrated, we have shown a general procedure to obtain resummed results to NLL 
for the rapidity distributions of a virtual photon 7* or of a real vector boson W^, Z^. 
Furthermore, we have outlined a general method to calculate numerical predictions 
and analyzed the impact of resummation for the fixed-target experiment E866/NuSea. 
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This shows that NLL resummation has an important effects on predictions of differ- 
ential rapidity cross sections giving an agreement with data that is better than NNLO 
full calculations. We have found a suppression of the cross section for not large values 
of hadronic rapidity instead of enhancing it. This suppression arises due to the shift 
in the complex plane of the dominant contribution of resummed exponent. These 
leaves open questions for future studies about possible suppression of the rapidity 
integrated cross sections at small x. 

The study of the renormalization group resummation applied to the case of small 
transverse momentum distribution of Drell-Yan pairs has opened further interestingly 
aspects about the relation between factorization properties of the cross section and 
the final structure of the resummed results which has been not yet well understood. 

Furthermore, because of its generality, renormalization group resummation lends 
itself naturally to some important future applications. They are the factorization of 
resummation of rare meson decay processes (like B —>■ Xs'j) and the resummation of 
generalized parton densities in deeply virtual proton Compton scattering. 



Appendix A 



Relations between logarithms of N 
and logarithms of {1 — z) 



In this Appendix we want to find the general relations between the logarithms of 
and the logarithms of {1 — z). Let's consider a generic logarithmical enhanced term 
in z space 

"ln^'(l - z) 



and take its Mellin transform: 



1-z 



1 ^^-1 - 1 



dz- 



1 - z 



\nP{l-z). 



(2) 



To find a general relation between the terms of the type of Eq. ([T]) and the logs of A^, 
we first notice that all integrals Ip Eq.([2]) can be obtained from a generating function 



dP 



T{N + r]) 7]' 

Now, using the Stirging expansion of the F function at large A^ 



(3) 
(4) 



we get 



r(Ar + 1) = V2;^e^"^^-^ + 



(5) 



G(r/) 



7] 







(6) 



At this point, we notice that A^ ^ is just the generating function of the log of A^: 

1 dP 



L„ = lnP^ = ^e"'"^/^ 



r]=0- 



(7) 
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Hence, Eq. ([6]) can be viewed as a relation between the generating function for Ip and 
for Lp. In particular, Taylor-expanding r(l + rf) in Eq.(l6]) leads to leading, next-to- 
leading,. . . In relations: 



G{ii) 



1 

7] 



oo 



k=0 



k\ ' 



k\ dln'1/NJo 
Now, if we put this last equation in Eq.(l3]), we obtain: 



^^e(^-i)in(i-z)_ 



-E 



rW(i) rf'^ 



p + 



-I p+i 

f 1 ^ V 



^ lnP(l -z) ^ / 1 



p+i^r('=)ri)(ini 



A:=0 



where in the last equality we have used the identity 
d" 



dln^l/N Jo 



N lnP(l - z) 



dz- 



1- z 



p + 1 \ k 



1 \ p+i-fc 

In — 

N 



(9) 
(10) 



This last expression is equal to zero when k > p + 1. This result can be expressed in 
terms of derivatives with respect to ln(l — z) with the identity 



d'' 



dz- 



dz d''\nP{l-z) 

dln^ 1/N Jo l-z ~ Jo 1-z d\n\l-z) 

thus obtaining for Ip Eq.© the following all- logarithmic-order relation: 

"^-^ dz dMnP(l-z) 



(12) 



^rW(i) 

fc=0 



k\ 

r(p+i)(i) 
P+i 







1- z dln^il- z) 
' 1 



(13) 



The inverse result, expressing Lp in terms of Jp, can be analogously found inverting 
the relation between the generating functions Eq.(l6]): 



riG{j]) + 1 

r(i + r7) ■ 



(14) 



Proceeding as before, we get 



(15) 



^ (k-1 



dz ■ 



Id" ln"fl-2) 



+A(")(1) + 
where A^''\r]) is the fcth derivative of 



{k-l)\ Jo l-z (iln^(l-z) 

1 ^ 
N 



(16) 



r(r^)' 



(17) 



Because of the p- independence of the coefficients of the expansion Eq. (|T^ . we 
can determine exphcitly the Melhn transform of a generic logarithmical enhanced 
function 



^(ln(l - zr) 



= 

+ p=0 



lnP(l - zY 
1 - z 



1 - z 

Indeed, its Melhn transform up to non-logarithmical terms is given by 

"1 ^N-l 



/ dz— ^(ln(l-^)") = 5^^,a% 

p=0 ^ k=Q ^ ^ ^ ^ 



p+l-fc 



(19) 



-E 

fc=0 



rW(i) /-i-* dz 



k\ 



I- zdln^il- z) 



^ ^ dz 



l-z 



g{\n{l - z) 



-g{\n{l - zf) 

""'dn r 
— ^(In-), 

n n 



(20) 



where in the last equality we have done the change of variable n = {1 — z) and 
where 

^rW(l)a'= d" 
^ In/C = - V M In/C . 



k=0 



(21) 



The inverse relation can be analogously derive. Namely, we can cast the integral of 
any function 

oo 

g{ln{l-zr) = Y,9p'^^''i^-^r (22) 
Mellin transform, up to non-logarithmic terms: 



— ^ In- 

n n 



1 ^7V-l _ I 

dz 

1 — z 



g{Hl-zr), 



(23) 



where 

fc=0 

For completeness, we recall that r'(l) = — 7^; = —0.5772 and that in some cases it 
can be useful to perform the change of variable n = n " so that we have 

1 Z"^" dn 1 dn' 1 

i / ^^(lni)= / ^^(ln4^) (25) 



Appendix B 
Phase space 



The n-hodj phase space can be expressed in terms of (n — m)-body and (m + l)-body 
phase spaces. To prove this statement, let us consider the definition of the phase 
space in d = 4 — 2e dimensions for a generic process with incoming momentum P and 
n on-shell particles in the final state with outgoing momenta Pi,P2, ■ ■ ■ ,Pn- 

iMP-.Pu . . . .PJ ^ . . . - p. - ■ ■ - PJ- (1) 

The momentum-conservation delta function can be rewritten as 

6d{P -Pi Pn) = J d''Q6d{P -Q-Pl Pm) 

xSd{Q-Pm+l Pn)- (2) 

The integration measure of this equation can be rewritten in this way: 
where 

q2 = (Q0)2-|g|2 (4) 

= (Pm+1 H VPnf (5) 

= {P-Pl Praf- (6) 

We shall now find the minimum and the maximum value of the variable Q^. If we 
use Eq.Q in the center-of-mass frame of the momentum Pm+i + ■ ■ ■ + Pn, we have 

The minimum value of is achieved when all the energies p^+i) • • • ^Pn equal to 
their invariant masses (mj = \/pf)- Hence 



Qmin — i\/ Pm+1 + " " " + yP- 
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because is a Lorentz ivariant. From Eq.(l6]), we have that in the center-of-mass 
frame of the momentum P — pi — • ■ ■ — p^, 

Q' = {P'-v\ V'mf- (9) 

Now, the minimum value of + ' ' ' + Pm is achieved when all these terms reduces to 
their invariant masses {mi = -\/pf). In this case 

Pi = ■■■ =Pm = (10) 
and this implies that P = and that P*^ = -\/P^. Therefore, we obtain 

QL. = (v^-^/p? (11) 

again as a consequence of the Lorentz invariance of Q^. We obtain immediately the 
general phase space decomposition formula using Eqs. ( fT|2|3l) together: 

d(j)n{P]Pl,-..,Pn) = / ——d(l)n,+i{P]Q,Pi,...,Pm) 

^ min 

xd(j)n-m{Q]Pm+l, ■ ■ ■ ,Pn), (12) 

where Q^in Q^ax given by Eq. (IHl) and Eq. ( ITTi) respectively. Using recursively 
Eq.( |T2i) with m = 1, it is possible to rewrite a n-body phase space in terms of n 
two-body phase spaces. We shall now compute the two-body phase space in the 
center-of-mass frame. We have 

In the center-of-mass frame P = and thus P" = \^P^. If we can neglect the invariant 
mass of p, we have that in this frame \Q\ = \p\ = p^ and that 



/P2 V 2VP2 



5[p'- -—^ . (14) 



/P2 V 2VP2 
Now, 

d'^-^p = ipi'^-^pidn,^^ = {py-^dp'^dn..^., (15) 

where dfl^-i is the solid angle in d—1 dimensions. Substituting Eqs. fll4fl5l) in Eq. (|T^ 
and integrating over p^, we obtain 

dMP;Q,p) = ^ [ ^ ^ ^^^-1 

= Nie)iP'r'U-j^) dn,.,, (16) 



where 

^(^) = 

Finally, we want to calculate dQ^-i- To do this, we use its definition 

dQd-i = d9d-ism''-^9d-idnd-2, (18) 

which can be applied recursively i times till d — 2 — i > with < 6d-i < vr. The 
normalization of the solid angles can be obtained performing the gaussian integral in 
spherical coordinates thus giving 

Qd-i = / dQd^i = — — -— , (19) 

J r((d- l)/2)' ^ ^ 

where F is the usual gamma function 

POO 

F(a) = / dte-'r-^ (20) 
Jo 

In our case (ci = 4 — 2e) we can use Eq.( IT8l) two times and we odtain 

dQa-i = dVL^^2e = sin^'^' 9 d9 sm'^' dcp dQi_2e- (21) 
In many cases it is useful to rewrite this equation in terms of other variables variables 

1 + cos 9 1 + cos (h , , 

1/1 = — ^ — , y2 = — 0<y, <1. (22) 

Doing this change of variables and recalling that the two-body squared amplitude 
cannot depend on more than two angles, we obtain 

= YiY/^ I' - ^i)]"' "^y^^y^^^ - y2r'^'~'. (23) 

We recall for completeness that the integrals in Eqs. (11211231) are indicated only to 
remind the integration range. They can be performed without the matching with 
the square amplitude of the corrispondig process only for the determination of their 
normalization. We can check explicitly that Eq. fl^Hl) is correct performing its integral. 
This is easily done using the definition of the B function 

B{z,w)^ ['dtt-\l-tr-' = ^^^f^, (24) 
Jo T{z + w) 

and the Legendre duplication formula 

T{2z) = {27t)-^/^2'''-'/^T{z)T{z + 1/2). (25) 

We find 

27r3/2-e 

which is in agreement with Eq. fll9l) . 



Appendix C 



Full NLO expression for DY 
rapidity distributions 



We report here the complete expression of the NLO DY distributions given 
[601 [611 [62] with the factorization scale equal to the renormalization scale: 



NLO 



da 
dQHY 



q,q i ^ 



X 



27r 



/TT ^ 
2 



27r 



where 



3^ 



for 7* 



5(g^ - M^) for Z° and , 
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and where 



with 



and 



Cgf{xi,X2,Y) = xiX26{xi - a;")5(a;2 -x^), x^2) = 
C^q (^Xi,X2, Y, = xiX2Cf!^S{xi - x°i)S{x2 - xl) 

+2Li2K) + 2U2{x'i) + \n\l - + ln2(l - + 2 In 

1 — x" 



x21n 



^0 
X2 



1 — X 



6{xi — x^ 
x\ ^x\^ -^x\ (\a.{l - x%lx2) 

X2 



^0 



„0 2 I 2 
•^2 ^ -^2 



X In — + ^ 



x\ 



X2 — X2 



X2 X2 X\{X2 — X2) 
+ 



T^0 2 I ™2 
X2 ^ -^2 



X2 



1 



X 



■In 



2x0(1 



X 



Ob n 



X\{X2 + 0:2) 
,0 ^0^ 



+ (1^2) 



{X2 - 

I G^{xi,X2,x\,xl) . ^OX , n ^ 

[(Xl - X\){X2 - X%)U ^ " 



1 — X 



1 - X?)(5(X2 - 4) [3 + 2 In + 2 In 

(^02 I ™2 1 \ N N 

p. 



r5(x2-xg 



x\ 



xl' + {x,-x\f) 

' ^ G''{x^,X2,Y) 
. . (2^2 - 3^2) + 

X2, Y) + 1^^| ^(^2-X°) ^^02 ^ _ ^0)2) 
11 1 ,Xi 



X 



CW(xi,X2,r) = CW(a:2,a:i,-F), 

ln(l - t) 



Li2(x, 

G^(a:i,a:2,F) = 

i/^(a;i,a;2,F) = 

G'^(xi,a:2,y) = 



^ - dt- 



Xlxl{xi+Xl)(x2 + X^) 

AxiX2{xiX2 + aj^a^O) 
,Tia;2(xiX2 + X2x'()'^ ' 
2a;2(a;?^a;2^ + {xiX2 — XiX2Y){xiX2 + x^Xa) 
xfx2(a;iX2 + a;2a;5)(a;2 + X2) ' 



2xTXn(xTXn + a;iX9)fa;ia;?Xo + x?a;nfa;ia;o + 2xtXo)) 



Appendix D 



Proof of some identities of chapter 



D.l Proof of Eq. (HI) 



We compute in the limit q'j_ ^ the distribution S{Mq) where Mq is given by Eq.( l28|l . 
Firstly we have 



= - - 6) + - 6 - 6))- (1) 

Then taking a generic test function f{C,i,C,2), 

J d^i J - - 6) + - ^1 - 6)) 
i + V 1-6 + 91 

dUi^U 1) / d^2 , / , f 1 - 6 - ^^^.2 



^^^^ 



i 1 - 6 + 91 



1-6 + 91 9i 



+ /^e /^V ; /^.r^ +/(l.l)ln^. (2) 
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Now, taking the limit ^ and using Eq.([T]), we obtain the following identity in 
the distribution sense 



71 + 71 n lng^5(l - - ^2) 



(3) 



which is exactly Eq. (l42i) . 



D.2 Proof of Eq. (143]) 



Let us compute for a generic test function /(^i,^2) the distribution 



(4) 



in the small-g^. limit with 77 = a|e|, a > and e 7^ 0. The integration range is fixed 
by the requirement 



6^: 



■[(i-ei)(i-6) + gl(i-ei-6)]>o, 



which gives 



or equivalently 



0<?,<1 o<&<ii-M^=& 



o<6<i o<e,<M#^ = 6- 



We observe also that 



6gi 



6gi 



Next, we decompose 



(5) 



(6) 



(7) 



where 



4 



(9) 



2 , (10) 



The ^2 integral in T?2 is immediately performed: 

5 I 



1 - 6 - 



Therefore, 



Tl 



■n-i 



1 

V 



(1-6 



Proceeding as above, we regularize the integral: 



where 



rp21 

hi - 


(1 








V 


hi - 


(1 








V 



(1-6 



i-6 + gi 



[/(6,i)-/(i,i)] 



/(i,i) / 



(1-6 



/(i,i) 



i-6 + gi 
(i + gir(i + gl)''-(gi)'' 



The two distributions T?2 and T?2 are now well defined as — > 0. Similarly, 



y I ^ y^ ^ ^ y^ ' 



where 



c?6(i-6 + gl) 



«2 



rf6 1 - 6 



6gi 



1-6 + gi 



{[/(6,6)-/(6,i)]-[/(i,6)-/(ii)]} 



^112 



^1 



c?6(i-6 + gIr' / ^6 l-6- 



[/(l,6)-/(l,l)] 



?2 



6gi 



i-6 + gi 



r,-l 



7,-1 



(11) 
(12) 

(13) 
(14) 



(15) 



(16) 



(17) 



(18) 

The term Tif can be computed by changing the order of integration and exploiting 
Eq.(l8]) and using eq ffTTl) . thus obtaining 



T'i = + ) 



2^r, rl 



di2 



(1-6 



■[/(i,6)-/(i,i)]. 



Collecting Eqs.( 



V Jo ^" 1 - 6 + ?! 

taking the limit — * and using the identity 

(i-er' = (i-or' + ^^(i-o, 



we obtain the following identity in sense of distributions: 

[(i-ei)(i-6) + gi(i-6-6)r' 



(19) 

(20) 
(21) 



-5(l-ei)(5(l-6) + 0(gl) 



which is exactly Eq. 



Appendix E 



Proof of combinatoric properties of 
chapter [71 



Let us consider the coefficients defined by 



-y2r""^-M "'-'""^"M. (1) 



We shall now prove that 



^(*,o) _ 2"'Cj^''^ (2) 



or, equivalently, that 



C^i'") = 2-(7(?'') (4) 
C^'^) = 2C^^~''>-Cfh z,j>l,, (5) 



where 



1=0 ^ / V J / 

Eq.(jll) follows immediately. Indeed: 

= + (7) 



m + i — 1 
i- 1 

We turn now to the proof of Eq.Q. In the case i = j = 1 it follows again 
immediately. In fact we have: 

C^y^ = 1, C^''^ = C^''^ = E 2' = 2"^+' - 1- (9) 

1=0 
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In the cases i — — 2ei — 2,j — litis quite easy: 

^(1,1) = 2^+^-1, C-^^^'O) = 2"^(m + 1) (10) 

^{1,2) ^ 2(2"^+^ - 1) - (m + 1), C^"'^) = 2"'+\m + 1) - (2"^+^ - 1), (11) 

where we have use the equality: 

V 12^ = lim y (e"^"2y ^ 2™+i(m + 1) - 2(2"*+^ - 1). (12) 

1=0 1=0 

In the other cases, i.e. for i,j > 2, it follows in a straightforward way: 

='5s--g2.'(''-:-)(™-;:/-)+g2.('t:-)("-;:- 



= c'^:'^^ (13) 

where V — \ and where in the second, the third and the last line we have used the 
following identities: 

r-2 j = ( .-1 )-\ r-1 )• 

_ r(r-l)_^^ j^^^ 



r - ly r(r)r(o) 

We note that these last two properties are vahd only for r > 2 and this is the reason 
why the cases i=j = lei = l,j = 2 and i = 2,_7 = 1 have been treated separately. 
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